Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,287,736 A→G 7.3% intergenic (‑365/‑28) acs ← / → nrfA acetyl‑CoA synthetase/nitrite reductase, formate‑dependent, cytochrome

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,287,7360AG 7.3% 237.0 / 8.4 84intergenic (‑365/‑28)acs/nrfAacetyl‑CoA synthetase/nitrite reductase, formate‑dependent, cytochrome
Reads supporting (aligned to +/- strand):  ref base A (47/30);  new base G (3/3);  total (50/34)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.03e-01

TAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCG  >  NC_000913/4287589‑4287883
                                                                                                                                                    |                                                                                                                                                   
tAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGcaaca                                                                                                                                                   <  1:875671/150‑1 (MQ=255)
tAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGcaaca                                                                                                                                                   >  2:481231/1‑150 (MQ=255)
    tatAAATTGGTCTTTTACATGCACATACAATTGATTAAAGACAACATTCACAGTGTTGTTATTTGTAACACATAGGGGCGAGCAATGTCATGACACTGTAGGTGCGGTTACTGTCGTGAAAAACACTAAAAACCGCCATTGCAACAATgg                                                                                                                                               <  1:209807/148‑1 (MQ=255)
     ctAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGgc                                                                                                                                              >  2:411289/1‑150 (MQ=255)
      tAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGgcg                                                                                                                                             >  2:1012136/1‑150 (MQ=255)
      tAAAGTGGTATTTTACATGCACTTACAATAGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGTTGCAATGgcg                                                                                                                                             <  2:847415/150‑1 (MQ=255)
        aaGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGTTGCAATGgcggt                                                                                                                                           <  2:326014/150‑3 (MQ=255)
         aGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAAACAATAAAAACCGCCATTGTTGCAATGgcggt                                                                                                                                           <  2:547218/150‑3 (MQ=255)
              tATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCgg                                                                                                                                     <  1:109347/150‑1 (MQ=255)
               aTTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGa                                                                                                                                    >  2:1060674/1‑150 (MQ=255)
                 tttACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGAtt                                                                                                                                  >  1:284919/1‑149 (MQ=255)
                 tttACATGCAATTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGGTACTGTCGTGAAAAACAATAAAAACCGCCATTTCAACAATGGCGCAATTCGGATg                                                                                                                                  <  2:505281/150‑1 (MQ=255)
                  ttACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCACCATTGCAACAATGGCGCAATTCGGATGa                                                                                                                                 >  1:456268/1‑150 (MQ=255)
                     cATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTc                                                                                                                                       >  2:74220/1‑141 (MQ=255)
                     cATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTc                                                                                                                                       <  1:74220/141‑1 (MQ=255)
                        gCACTTACAATTGATTAAAGACAACATTCACAGTGTGGGTATTTGTTACACATAGGGGCGAGCAATGTCATGCCAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACTGCCATTGCAACAATGGCGCAATTCGGATGAAGcccc                                                                                                                           <  2:851352/150‑1 (MQ=255)
                         cACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGTGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGAAATTCGGATGAAGCCCCt                                                                                                                          >  1:573841/1‑150 (MQ=255)
                            ttACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATg                                                                                                                       <  1:993377/150‑1 (MQ=255)
                               cAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGTGCGAGCAAAGTCATGACACTGTAGGTGAGGTTACTGTCTTTAAAAACAATAAAAACCTCCATTGCAACAATGGCGCAATTCCGATGAAGCCCATATGACa                                                                                                                    >  1:208970/1‑150 (MQ=255)
                                aaTTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCGATGACaa                                                                                                                   >  2:767797/1‑150 (MQ=255)
                                 aTTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACaaa                                                                                                                  >  1:820495/1‑149 (MQ=255)
                                 aTTGAATAAAGACAACATTCACAGTGTGCTAATTTGTTACACATACGGGCGAGCAAAGTCATGACAGTGTAGGTGCTGTTACTGTCGTGAAAAACAAAAAAAACCGCCATAGCAACAATGGCGCAATTCGTATGAATCACCTATTAAAAg                                                                                                                  >  1:810716/1‑150 (MQ=255)
                                  ttGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAgg                                                                                                                 <  1:863452/150‑1 (MQ=255)
                                    gATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGAGGTTACTTTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAATGAt                                                                                                               >  1:58094/1‑150 (MQ=255)
                                     aTTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATa                                                                                                              >  2:605612/1‑150 (MQ=255)
                                      ttAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaa                                                                                                             >  2:712908/1‑150 (MQ=255)
                                       tAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaa                                                                                                            >  1:582692/1‑150 (MQ=255)
                                       tAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGGCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaa                                                                                                            <  1:320393/150‑1 (MQ=255)
                                        aaaGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCATTTCGGATGAAGCCCGTATGACAAGGTTaaaa                                                                                                           >  2:184851/1‑150 (MQ=255)
                                           gACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATa                                                                                                        >  1:496530/1‑150 (MQ=255)
                                              aacaTTCACAGTGTGGTTATTTGTTACACATAGTGGCGATCAATGTCATGACATTGTAGGTTCGGTTACTGTAGTGAAAAACAATAAAAACCTCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAac                                                                                                     <  1:246798/150‑1 (MQ=255)
                                                        gtgtGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTa                                                                                           >  1:852385/1‑150 (MQ=255)
                                                        gtgtGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTa                                                                                           >  2:137941/1‑150 (MQ=255)
                                                        gtgtGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTa                                                                                           >  1:322150/1‑150 (MQ=255)
                                                          gtgGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATc                                                                                         <  1:1116427/150‑1 (MQ=255)
                                                            ggTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGAGTACGTGCGGTTACTGTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCTATTCGGATGAATCCCCTATGACAAGGAAAAAAATAAACGCACGCCGTATAtt                                                                                       >  2:1172713/1‑150 (MQ=255)
                                                              ttATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCa                                                                                     <  2:463174/150‑1 (MQ=255)
                                                                 tttGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCt                                                                                  <  1:229819/150‑1 (MQ=255)
                                                                  ttGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCtt                                                                                 <  1:816656/150‑1 (MQ=255)
                                                                     ttACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAtt                                                                              <  1:712908/150‑1 (MQ=255)
                                                                     ttACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAAATCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCGTCCGCTTAtt                                                                              >  2:578653/1‑150 (MQ=255)
                                                                      tACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTg                                                                             >  1:573135/1‑150 (MQ=255)
                                                                           aTAGGGGCGAGCAATGTCATGACAGTTTTGGTGCGGTTACTGTCGTGAAAAACAATAAAAACAGCCATTGCAACAATGTCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTcc                                                                        <  2:922833/150‑1 (MQ=255)
                                                                               gggCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtttt                                                                    <  1:605612/150‑1 (MQ=255)
                                                                               gggCGAGCAATGTCATGACAGTGTAGGTGCGGATACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAATCCCCTTTGACAATGATAAAAATAAACGCACTCCGTATCTTCATCTTATTGATTcgattt                                                                    >  2:71511/1‑145 (MQ=255)
                                                                                ggCGAGCAAGGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGTTTAAAATAAACGCACGCCGTTTCTTCAGCTTATTGTTTCGttttt                                                                   >  1:625857/1‑150 (MQ=255)
                                                                                 gCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtttttt                                                                  <  1:1146452/150‑1 (MQ=255)
                                                                                   gAGCAATGTCATGACAGTGTAGGTGCGGTTTCTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGGGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttct                                                                <  1:1133121/150‑1 (MQ=255)
                                                                                     gCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttt                                                              <  1:994379/150‑1 (MQ=255)
                                                                                       aaTGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATTAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt                                                            >  2:689553/1‑150 (MQ=255)
                                                                                       aaTGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAACCGCCCGCCGTATCTTCAGCTTATTGATTCCTtgtttcttttt                                                            >  2:170604/1‑150 (MQ=255)
                                                                                       aaTGTCATGACAGTGTAGGTGAGGTTACTGTCGTGAAAAAAAATAAAAACCGCCATTGCAACAATGGAGCAATTCGGATGAAGCCCCTAAGACAAGGATAAAAATAAACGCACGCAGGATCTTCAGCTTATTGATTCCTtgtttcttttt                                                            >  2:166900/1‑150 (MQ=255)
                                                                                          gTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCAc                                                         >  2:351646/1‑150 (MQ=255)
                                                                                          gTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCAc                                                         >  1:733090/1‑150 (MQ=255)
                                                                                          gTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCAc                                                         >  1:386860/1‑150 (MQ=255)
                                                                                            cATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACtt                                                       >  2:553026/1‑150 (MQ=255)
                                                                                              tGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGCAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGGATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCt                                                     >  1:348745/1‑150 (MQ=255)
                                                                                              tGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCCGCTTATTGATTCCTTTTTTCTTTTTCCATTCt                                                     >  2:124374/1‑150 (MQ=255)
                                                                                              tGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAATGAGTCCTTTTTTCTTTTTCACTTCt                                                     >  1:892309/1‑150 (MQ=255)
                                                                                                    tgtAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCATTTTCACTTCTGTTCAc                                               >  2:150773/1‑150 (MQ=255)
                                                                                                       tggTGCGGTTACTGTAGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCt                                            <  2:573841/149‑1 (MQ=255)
                                                                                                       aGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATGCCTTTTTTATTTTTAACTTCTGTTCACGCt                                            >  2:243603/1‑150 (MQ=255)
                                                                                                        ggTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTg                                           >  1:710170/1‑150 (MQ=255)
                                                                                                          tGCTGTTACTTTCGTGAAAAACAATATAATCTGCCTTTGATTCATTGGCGCAATTCGGATGAAGCCCCTATGACAAGGATTAAAATATAAGCACGCCGTATCTTCAGGTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGaa                                         <  2:810716/150‑1 (MQ=255)
                                                                                                                  cTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgc                                 <  2:436679/150‑1 (MQ=255)
                                                                                                                  cTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgc                                 <  2:1032415/150‑1 (MQ=255)
                                                                                                                  cTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgc                                 >  2:1239184/1‑150 (MQ=255)
                                                                                                                  cTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgc                                 <  1:1315567/150‑1 (MQ=255)
                                                                                                                    gTCGTGAAAAACAATTAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAACCGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctg                               <  1:184192/150‑1 (MQ=255)
                                                                                                                    gTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctg                               <  1:711286/150‑1 (MQ=255)
                                                                                                                     gcGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctgc                              <  2:519689/149‑1 (MQ=255)
                                                                                                                       gTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTc                            >  1:268487/1‑150 (MQ=255)
                                                                                                                             aaaCAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGCTTCCTTTTTTCTTTTTCACTTCTGATCACGCGGAACAAACGGCTGCGCCCGCa                       >  1:243559/1‑150 (MQ=255)
                                                                                                                             aaaCAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTATTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCa                       >  2:115819/1‑150 (MQ=255)
                                                                                                                              aaCAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCaa                      >  1:996274/1‑150 (MQ=255)
                                                                                                                                      aaaCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTa              <  2:58094/150‑1 (MQ=255)
                                                                                                                                       ttgcGCCATTGTTGCAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAACCGGCTGCTCCCGCAAAACCTGTaa             >  2:952061/4‑150 (MQ=255)
                                                                                                                                       aaCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTaa             >  1:824038/1‑150 (MQ=255)
                                                                                                                                        tgcGCCATTGTTGCAATGGCGCAATTCGGATGACGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTATTTTTCACTTCTGTTCACGCTGCACAAACGGCTGCTCCCGCAAAACCTGTAAc            >  1:635260/3‑150 (MQ=255)
                                                                                                                                        tgcGCCATTGTTGCAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTAACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAAc            >  1:977409/3‑150 (MQ=255)
                                                                                                                                         ccGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTACCt           >  1:1165076/1‑150 (MQ=255)
                                                                                                                                            ccATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACtgtg        <  1:1060674/150‑1 (MQ=255)
                                                                                                                                             cATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTgg       >  1:296452/1‑150 (MQ=255)
                                                                                                                                              aTTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGa      <  2:852385/150‑1 (MQ=255)
                                                                                                                                                tGCAACAATGGCGCAATTCGGCTGAAGCCCCTATGACAAGGATAAAAATAAACTCACTACTTTTCTTAAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAg    <  2:369064/150‑1 (MQ=255)
                                                                                                                                                tGCAACAATGGCGCAATTCGGATGAAGCCCCTATGCCAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCATTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAg    >  2:828074/1‑150 (MQ=255)
                                                                                                                                                  caacaaTGGCGCAATTCGGATGAAGCCCCTATGCCAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCg  >  1:442673/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCG  >  NC_000913/4287589‑4287883

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: