Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,823,288 C→T G161D (GGC→GAC)  recA ← DNA recombination/repair protein RecA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,823,2880CT96.4% 176.5 / NA 55G161D (GGC→GAC) recADNA recombination/repair protein RecA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (24/29);  minor base G (0/2);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

CGCTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTG  >  NC_000913/2823046‑2823522
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
cGCTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACATAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATCCGGATCTGGTTGATGAAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTc                                                                                                                                                                                                                                     <  1:94989/250‑1 (MQ=255)
 gCTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGc                                                                                                                                                                                                                                   <  1:273369/251‑1 (MQ=255)
   ccccccccccGTTTTCCCCCTCTTTCACCGCGCCGATTCGACGGATGTCGAGACAAACAGAGGCGTAGACTTTCAGCTCGTTACCTCCGTGAGTGGTTTCCGGGTTACCTAACATCACACCAATTTTCATCCGGATCTGGTTGAGGAAAATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCTGCTGGCTCATCATACGTGCCGCAAGCCCCAGGTGAGAGTCGGCGATTTCGCCt                                                                                                                                                                                                                                 <  2:406886/242‑1 (MQ=255)
          cacGTTTTCGCCCTCTTCCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTGCGCCTTCGAttt                                                                                                                                                                                                                          <  2:239700/251‑1 (MQ=255)
            ctttttCGCCCTCTTTCCCCGCGCCGATTCGACGGATGTCGCGGAGAACTGAGTAGTAAATTTTCAGCGCGTTCCCTCCTGTAGTGGTTTCCGGTTTACCGACCATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCA‑GCAGCGTTTTGGCCTGCTGCAGGGTACCCTCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGAttt                                                                                                                                                                                                                          <  2:215675/247‑1 (MQ=255)
               tttCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCAGACGGATCTGGTTGATGAAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGc                                                                                                                                                                                                                      <  2:136774/250‑1 (MQ=255)
                   gCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACCGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCg                                                                                                                                                                                                                 >  1:368621/1‑251 (MQ=255)
                            cACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGACGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTg                                                                                                                                                                                                        >  1:121013/1‑251 (MQ=255)
                                     gATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTGCATACGGATCTGGTTGATGAAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACgg                                                                                                                                                                                               <  1:68240/251‑1 (MQ=255)
                                         cgacgGCTGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTc                                                                                                                                                                                           >  1:284274/1‑251 (MQ=255)
                                         cgacgGATGTCGAGACGAACAGAGGCGTATAATTTCAGCGCGTTACCTCCGGTAGTGGTTTCCGGTTTACCGAACATCACACCAATTTGCATACGGATCTGGTTGATGAAGATCA‑GCAGCGTTTTGGACTGCTTCAGGTTGCCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGGCGATTGCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTc                                                                                                                                                                                           <  1:36280/251‑1 (MQ=255)
                                            cgGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGAGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCCTACGGATCTGGTTGATGAAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAac                                                                                                                                                                                        >  1:386353/1‑251 (MQ=255)
                                             ggATGTCGAGCCGAACGAAGGCGTAGAATTACAGCGCGTTACCACCGGTAGTGGTTTCCGGGTGACCGAACATCACCCCAATTTTCATACGGATGTGGTTGATGAAGATCA‑GCAGCGTGGTGGACTTCTGCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGGCGATGTGAGAGACGTCGATGTTGCTTTCGATTTTCGCTTTCGGCGTCAGTGCCGCCACGGCGTCAacg                                                                                                                                                                                       <  2:39203/251‑1 (MQ=255)
                                                                             agggcgTTCCCACCGGTAGTGTTTTCCGGGTTACCGAACACCCCACCAATTTTCATACGGATCTGGTTGAAGAAGACCA‑GCGGCGTTTTGGACTGCGTCAGGTTCCCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGAGTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACTTCTACTGCGCCAGAACGCGCCAg                                                                                                                                                       <  2:89496/248‑1 (MQ=255)
                                                                             aGCGCGTTACCACCGGTAGTGGTTTCCGGTTTACCGAACATCACCCCAATTTTCATACGGATCTGGTTGATGAAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAg                                                                                                                                                       <  1:145842/251‑1 (MQ=255)
                                                                               cgcgTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAggg                                                                                                                                                     >  1:177544/1‑251 (MQ=255)
                                                                                  gTTACCCCCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGt                                                                                                                                                  <  1:140322/251‑1 (MQ=255)
                                                                                       caccGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAg                                                                                                                                             <  1:21396/251‑1 (MQ=255)
                                                                                        accGGCAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTcac                                                                                                                                               >  2:242717/1‑248 (MQ=255)
                                                                                              tgTGGTTTCCGGTTTACCGACAATCATCCGAATTTTCATACGGATCGGTTTTATGAATA‑CACACAGCGTGTTGGACTGCTTCAGGTTCCCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACTCGCCAGGGCGTCACCGAttgcc                                                                                                                                       <  2:230699/249‑4 (MQ=255)
                                                                                               gTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTcc                                                                                                                                       >  2:343856/1‑249 (MQ=255)
                                                                                                    ttCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCt                                                                                                                                >  1:373969/1‑251 (MQ=255)
                                                                                                        gggatACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGAATTCCAGTGCCTGCTc                                                                                                                            >  2:68465/5‑251 (MQ=255)
                                                                                                                                ttCATACGGGTCTGGTTGAGGAAGGTCG‑TCAGCGTGTTGGACTGCTTCAGGTTCCCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGGGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGcaca                                                                                                     <  2:18436/250‑1 (MQ=255)
                                                                                                                                  cATACGGATCTGGTTGATGAAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACTACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAACAcagc                                                                                                   <  1:121622/250‑1 (MQ=255)
                                                                                                                                    tACGGATCTGGTTGATGAAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAgg                                                                                                >  1:20003/1‑251 (MQ=255)
                                                                                                                                       ggATCTGGTTGATGAAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTg                                                                                             >  1:262910/1‑251 (MQ=255)
                                                                                                                                        ggagctgggggatgaagagatgcagtGTGTTGGTCTGCTTCAGGTTTCCCGCCAGCTTACGCATCGCCTGGCTCATCCTACGGGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCGTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGCCGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTg                                                                                             <  2:285284/230‑1 (MQ=255)
                                                                                                                                                tgatgaAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATc                                                                                     <  1:382014/250‑1 (MQ=255)
                                                                                                                                                    gaAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACg                                                                                 <  1:321687/250‑1 (MQ=255)
                                                                                                                                                      aGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGCTTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGCTATCGACGccc                                                                              >  2:307662/1‑251 (MQ=255)
                                                                                                                                                               agcgtgaaggccTGCTTCAGGTTACCAGCAAGCTTACGACTCGCCTGGCTAATCATACGTGCCGCACGGACCATGTGAGAGACGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGCTAACGTCTACTGCGCCCGAACGCGCCAGGGCGTAACAGACTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGCTCTCGCCGCCCCGTTTACg                                                                      >  2:61690/12‑251 (MQ=255)
                                                                                                                                                                 cGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCCAGGCCCCTGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCCGTGCCGCCCAGGAGTCAACCACGACAACGGCTACTGCGCCAGAACGCGCCCGGGCGTCACCGCTTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGGTATCGACGCCCCGGTTACGTg                                                                    >  2:215582/1‑251 (MQ=255)
                                                                                                                                                                    gtTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGt                                                                 >  2:128075/1‑251 (MQ=255)
                                                                                                                                                                    gtTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGGCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCgg                                                                 >  1:274944/1‑250 (MQ=255)
                                                                                                                                                                       ggACTGCTGCAGGTGACCCGCCAGCTTACGCAGCGCCTTGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCTTTCTATTTCCGCTTTCGGCGTCAGTGCCTCCACGGAGTCAACGACGATATCTTCTACTGCGCCAGAACGCGCCATGGCGTCACAGTTTTCCAGTTCCTGCTCGCCGGTGTCCGGTTGGTAGCTCAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGa                                                              <  1:181975/251‑1 (MQ=255)
                                                                                                                                                                                cAGGTTGCCCGCCAGCTTACGCATCGCCGGGCTCATCATACGTTCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGGCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCa                                                     <  2:344829/251‑1 (MQ=255)
                                                                                                                                                                                 aGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGAGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTcc                                                      >  1:231525/1‑249 (MQ=255)
                                                                                                                                                                                      aCCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGcgcg                                                <  1:354114/250‑1 (MQ=255)
                                                                                                                                                                                      aCCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCgt                                               >  1:130870/1‑251 (MQ=255)
                                                                                                                                                                                              gCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCa                                       <  1:61813/251‑1 (MQ=255)
                                                                                                                                                                                                      aTCGCCTGGCTCATCAAACGTGCCGCAAGGCCCATGTGAGAGTCGTCGCTTTCGCCTTCGATTTCCGCTTTCGGCGTCCGTGCCGCCACGGAGTCAACGACGCTAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATaa                                >  2:187585/1‑250 (MQ=255)
                                                                                                                                                                                                         gCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGc                            >  1:344883/1‑251 (MQ=255)
                                                                                                                                                                                                            tGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGCTAAACGcac                          >  2:359166/1‑250 (MQ=255)
                                                                                                                                                                                                                 catcatACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTTCAGCGCGTGTTCAGCATCGATAAACGCACAGGttt                    <  1:53452/251‑1 (MQ=255)
                                                                                                                                                                                                                 catcatACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGttt                    <  1:287760/251‑1 (MQ=255)
                                                                                                                                                                                                                 catcatACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGtt                     <  1:209627/250‑1 (MQ=255)
                                                                                                                                                                                                                       ccGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGGGCCGCCACGGAGTCAACGACGATAACGGCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTGCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTGTACCtt              <  2:368621/250‑1 (MQ=255)
                                                                                                                                                                                                                       aCGTGCCTCAAGGCCCATGTGAGGGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTTCCGCCACGGAGTCAACGACTATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCTGTGCCTGCTCGCCGTTGGCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTGACGTGCGTAGATTGGGTCCGGCGCGTGTTCAGCATCGATAAACGCACAGGCTTTACCt               <  2:376759/250‑1 (MQ=255)
                                                                                                                                                                                                                       aCGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTgg                                                                                                            >  2:305183/1‑157 (MQ=255)
                                                                                                                                                                                                                       aCGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTgg                                                                                                            <  1:305183/157‑1 (MQ=255)
                                                                                                                                                                                                                           gCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCAcg          <  1:343856/251‑1 (MQ=255)
                                                                                                                                                                                                                              gCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCAcgctg       >  2:164493/1‑251 (MQ=255)
                                                                                                                                                                                                                               cAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTAATGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCAcgctgc      >  1:206460/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCAcgctgcgctg  <  1:215582/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
CGCTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCA‑GCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTG  >  NC_000913/2823046‑2823522

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑