Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,256 G→A 7.2% intergenic (+265/+142) 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,2560GA 7.2% 236.4 / 6.3 83intergenic (+265/+142)fdnI/yddMformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (45/32);  new base A (4/2);  total (49/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

CCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCG  >  NC_000913/1552109‑1552405
                                                                                                                                                   |                                                                                                                                                       
ccTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcg                                                                                                                                                       <  1:735454/150‑1 (MQ=255)
            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:1371678/1‑150 (MQ=40)
                                                                        tcatcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAg                                                                             <  2:538509/150‑1 (MQ=40)
                                                                        tcatcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAg                                                                             >  1:1493119/1‑150 (MQ=40)
                                                                         catcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc                                                                            >  2:708604/1‑150 (MQ=40)
                                                                         catcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc                                                                            >  1:184126/1‑150 (MQ=40)
                                                                          ttcaGTTTCAAATGGCGCTGTAAAAGGCGGCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCg                                                                           <  2:1661148/149‑1 (MQ=25)
                                                                          atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCg                                                                           >  2:920781/1‑150 (MQ=39)
                                                                          atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCg                                                                           >  2:559492/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                                                                   >  2:499453/1‑150 (MQ=32)
                                                                                  cAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAt                                                                   >  1:134692/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:1639463/1‑150 (MQ=32)
                                                                                             gTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGtt                                                        >  1:828752/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  >  2:490370/1‑150 (MQ=255)
                                                                                                                                                   gcggcgTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCg  <  1:1080075/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                       
CCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCG  >  NC_000913/1552109‑1552405

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: