Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,270 T→G 8.6% intergenic (+279/+128) fdnI → / ← yddM formate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,552,2700TG 8.6% 221.6 / 11.3 90intergenic (+279/+128)fdnI/yddMformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (47/33);  new base G (6/2);  total (55/35)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

ATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGC  >  NC_000913/1552121‑1552417
                                                                                                                                                     |                                                                                                                                                     
aTAGAATGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGAATCAGGCAATGTTGCTTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                                                                                                       <  2:1323303/150‑1 (MQ=255)
 tAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGAGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                                      <  1:380749/150‑1 (MQ=255)
      tgctttGAACGTTTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGCTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                                 <  1:112581/145‑1 (MQ=255)
        cAGTGAACGTGTAGCCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                               <  1:1592737/150‑1 (MQ=255)
         ctcgAACTTGTAGTCCTGATACGCGTAGCCCATCATTAAATGTTTATTTTTTATTAAGTTTCAAATGGCGCTGTAAAAGGCGTCTTTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATACGGCGTTCACGCCGCATCCGGCa                                                                                                                                              <  2:203596/147‑1 (MQ=9)
         cgTGAACGAGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGGTGCGTTTGTCATCACTTTCAATTGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                                                                                              <  1:1315799/149‑1 (MQ=255)
          ttgcaCGTGTAGTCATTATAAGCGTAGCGCATCAGGCAATGTTTCGTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCAt                                                                                                                                             <  2:1250807/146‑1 (MQ=18)
           tGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATg                                                                                                                                            <  1:1001128/150‑1 (MQ=25)
            gAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCGTATCCGGCATGa                                                                                                                                           >  2:360196/1‑150 (MQ=255)
               cGTGTAGGCCTGATAAGCGGAGCGCATCAGGCAATGTTTCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGa‑‑aaa                                                                                                                                      <  2:892823/150‑1 (MQ=255)
                  gtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCATGAACGCCTTATCCGGCATGA‑‑AAAccc                                                                                                                                   <  1:1302945/150‑1 (MQ=255)
                   tagccatGATATCCGTAGCTCATCATTAAATCTTTCGTTTCCCTTCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCt                                                                                                                                  <  1:191458/144‑1 (MQ=14)
                    aGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGAATACAACGTAGGCTTAATTCATTCCGGAAGCCGCGTGAAACCCTTATACCGAATGA‑‑AAACCCta                                                                                                                                 >  2:812918/1‑149 (MQ=21)
                     gtccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTc                                                                                                                                <  2:1169502/148‑1 (MQ=255)
                     ggCCAGATAAGCCTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATTTCGATGTAAAACGCGTAATTTTCATATTTTAGAAATCGTAGGATTTGTTCATGCAGGATGCGGCGTTAACGCCTTATCCGGCATGA‑‑AAACCCTTc                                                                                                                                <  2:324650/150‑1 (MQ=18)
                      gCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCa                                                                                                                               <  1:1521955/150‑1 (MQ=255)
                      gCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCa                                                                                                                               >  1:728282/1‑150 (MQ=255)
                       ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCaa                                                                                                                              >  1:1245197/1‑150 (MQ=255)
                          taTAACCGTAGCGCATCAGTCAATGTTTCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATc                                                                                                                           <  1:1512423/149‑1 (MQ=255)
                            tAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAACGTAGGCTTTGTTCATTCCGGATGAGGCGTGAAAGACATATACGGAATGA‑‑AAACCCTTaaaaacaa                                                                                                                         >  1:1192012/1‑146 (MQ=18)
                             aaGCGTAGGGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGAGTCATTTTGATATTGTATACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTAAAATCCaa                                                                                                                        >  1:939690/1‑150 (MQ=18)
                              aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAAt                                                                                                                       <  2:532518/150‑1 (MQ=255)
                              aGCGTAGCGCATCAGGAAATGTTGCGTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAAt                                                                                                                       <  2:240752/150‑1 (MQ=255)
                               gCGTAGCGCATCATGCAATGTTTCGTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATa                                                                                                                      <  2:873335/150‑1 (MQ=255)
                                  tAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGAt                                                                                                                   <  1:815112/150‑1 (MQ=255)
                                     cgcttgaGGAAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGc                                                                                                                <  2:506873/144‑1 (MQ=255)
                                     cgcATCTGGCAATGTTTCGTTTTTTTTCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGc                                                                                                                <  2:1709882/150‑1 (MQ=14)
                                           aGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTATAAGGAGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGTCGTGAAAGACTTTTCCGGCATCA‑‑AAAAAATACAAATCCAATAGATTCCAGAGaa                                                                                                          >  1:26748/1‑150 (MQ=17)
                                           aGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCAAATCCGGAATGA‑‑AAACACATCAAATCCAATAGATTGCAGTGaa                                                                                                          >  2:704112/1‑150 (MQ=17)
                                           aGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGAGGCGTGAAAGCCTTATCAGGCATGA‑‑AAACCCTTCAAATACAATAGATTGCAGTGaa                                                                                                          >  1:1189759/1‑150 (MQ=17)
                                            ggCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGAGACATTTTAATATTGTAGAAAAAGTAGGCTTTGTACATGCCGGATGCCGCGAGAACCACTAAACAGCCATCA‑‑AAACCCTTAAAAAACAATAGATTTCAGTGAAc                                                                                                         >  2:354556/1‑150 (MQ=255)
                                                            tcatcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAg                                                                                         <  2:538509/150‑1 (MQ=40)
                                                            tcatcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAg                                                                                         >  2:1371678/1‑150 (MQ=40)
                                                            tcatcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAg                                                                                         >  1:1493119/1‑150 (MQ=40)
                                                             catcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc                                                                                        >  1:184126/1‑150 (MQ=40)
                                                             catcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc                                                                                        >  2:708604/1‑150 (MQ=40)
                                                              ttcaGTTTCAAATGGCGCTGTAAAAGGCGGCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCg                                                                                       <  2:1661148/149‑1 (MQ=25)
                                                              atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCg                                                                                       >  2:559492/1‑150 (MQ=39)
                                                              atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCg                                                                                       >  2:920781/1‑150 (MQ=39)
                                                                  gTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc                                                                                   >  2:742241/1‑150 (MQ=35)
                                                                  gTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc                                                                                   >  1:1213114/1‑150 (MQ=35)
                                                                  gTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc                                                                                   >  1:667670/1‑150 (MQ=35)
                                                                    ttCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcgc                                                                                 >  2:1080075/1‑150 (MQ=34)
                                                                    ttCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcgc                                                                                 >  1:362626/1‑150 (MQ=34)
                                                                      cAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAt                                                                               >  1:134692/1‑150 (MQ=32)
                                                                      cAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAt                                                                               >  2:499453/1‑150 (MQ=32)
                                                                       aaaTGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCaaa                                                                                                                             <  2:145470/103‑1 (MQ=32)
                                                                        aaTGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑TAACCCTTCAAATCCAATAGATTGAAGTGAAAGTGTAGGCCTGATAAGCGTAGCGCAtaa                                                                             >  1:1325480/1‑148 (MQ=12)
                                                                            gcgcTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                         >  1:378925/1‑150 (MQ=32)
                                                                                 gTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGtt                                                                    >  1:828752/1‑150 (MQ=32)
                                                                                 gTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGtt                                                                    >  1:1639463/1‑150 (MQ=32)
                                                                                       ggCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGttt                                                              >  1:1278006/1‑150 (MQ=32)
                                                                                        gCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTg                                                             >  1:619824/1‑150 (MQ=32)
                                                                                           tCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGtca                                                          >  2:1529212/1‑150 (MQ=11)
                                                                                                    tatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTcaa                                                 >  2:146836/1‑148 (MQ=11)
                                                                                                     atTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTaa                                                >  1:678739/1‑150 (MQ=32)
                                                                                                     atTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGAGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTaa                                                >  2:1184583/1‑150 (MQ=17)
                                                                                                        gTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATg                                             >  2:1672389/1‑150 (MQ=34)
                                                                                                        gTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATg                                             <  1:146836/150‑1 (MQ=34)
                                                                                                        gTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCATCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATg                                             >  1:916781/1‑150 (MQ=21)
                                                                                                          aGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGgc                                           >  2:437929/1‑150 (MQ=34)
                                                                                                           gACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGgcg                                          <  1:1248704/150‑1 (MQ=34)
                                                                                                                cGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTta                                     >  2:1487660/1‑150 (MQ=37)
                                                                                                                cGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTGta                                     >  1:359912/1‑150 (MQ=17)
                                                                                                                    ggCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCCGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATaaa                                 >  2:116041/1‑150 (MQ=255)
                                                                                                                    ggCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATaaa                                 >  1:716824/1‑150 (MQ=40)
                                                                                                                      cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGt                               >  2:1327445/1‑150 (MQ=255)
                                                                                                                        ttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGc                             <  2:766959/150‑1 (MQ=255)
                                                                                                                           ttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAt                          >  1:1212521/1‑150 (MQ=255)
                                                                                                                             cATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttt                        <  2:1493119/150‑1 (MQ=255)
                                                                                                                             cATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttt                        >  2:1098384/1‑150 (MQ=255)
                                                                                                                              aTGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtt                         >  2:615037/1‑150 (MQ=14)
                                                                                                                              aTGCCGGATAAGGCGTTCACGCCGCATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGAGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttt                       >  2:284916/1‑150 (MQ=14)
                                                                                                                              aTGCCGGATAAGGCGTTCACGCCGCATCCGGCATGA‑‑AAACCCTTCAAATCAAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttt                       >  1:573210/1‑150 (MQ=14)
                                                                                                                               ttccGGATGCGGCGTTAACACCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttttt                      <  1:959944/148‑1 (MQ=255)
                                                                                                                                gCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttttt                     >  2:916814/1‑150 (MQ=255)
                                                                                                                                gCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAtttaga                     <  2:1278006/150‑4 (MQ=255)
                                                                                                                                 ccGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttttttt                    <  1:1633050/150‑1 (MQ=255)
                                                                                                                                 ccGGATGCGGCATGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttttttt                    >  2:1484388/1‑150 (MQ=255)
                                                                                                                                 ccGGATAAGGCGTTCACGCCGCATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTTTAAAAAGAGACAtttttaa                    >  2:959944/1‑148 (MQ=11)
                                                                                                                                  ctgATGCGGCGTGAACGCCTTATCCTGCATGA‑‑AAACCCTTCAAATCGAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTa                   <  1:82822/148‑1 (MQ=255)
                                                                                                                                  cGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTa                   <  2:1409962/150‑1 (MQ=255)
                                                                                                                                    gATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAtt                 <  2:440604/150‑1 (MQ=255)
                                                                                                                                       gcggcgTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCg              <  1:1080075/150‑1 (MQ=255)
                                                                                                                                       gcggcgTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCg              >  2:490370/1‑150 (MQ=255)
                                                                                                                                          gcgTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTaa           >  1:332112/1‑150 (MQ=255)
                                                                                                                                           cgTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAAc          >  2:805200/1‑150 (MQ=255)
                                                                                                                                           cgTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAAc          <  1:920781/150‑1 (MQ=255)
                                                                                                                                              gAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCac       >  2:1155046/1‑149 (MQ=255)
                                                                                                                                               caCGCCGCATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGa      >  1:1082809/2‑150 (MQ=255)
                                                                                                                                                  gCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACgg   <  2:201381/150‑1 (MQ=255)
                                                                                                                                                   ccTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGc  >  2:1035300/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGC  >  NC_000913/1552121‑1552417

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