Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,287,735 A→T 7.9% intergenic (‑364/‑29) 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,7350AT 7.9% 216.6 / 7.7 84intergenic (‑364/‑29)acs/nrfAacetyl‑CoA synthetase/nitrite reductase, formate‑dependent, cytochrome
Reads supporting (aligned to +/- strand):  ref base A (47/30);  new base T (3/4);  total (50/34)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

GACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCG  >  NC_000913/4287586‑4287883
                                                                                                                                                      |                                                                                                                                                    
gACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGca                                                                                                                                                      >  1:463174/1‑150 (MQ=255)
   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)
                                                                                                                                                      |                                                                                                                                                    
GACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCG  >  NC_000913/4287586‑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: