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

TAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCG  >  NC_000913/4287589‑4287883
                                                                                                                                                    |                                                                                                                                                   
tAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGcaaca                                                                                                                                                   <  1:875671/150‑1 (MQ=255)
tAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGcaaca                                                                                                                                                   >  2:481231/1‑150 (MQ=255)
    tatAAATTGGTCTTTTACATGCACATACAATTGATTAAAGACAACATTCACAGTGTTGTTATTTGTAACACATAGGGGCGAGCAATGTCATGACACTGTAGGTGCGGTTACTGTCGTG‑AAAAACACTAAAAACCGCCATTGCAACAATgg                                                                                                                                               <  1:209807/148‑1 (MQ=255)
     ctAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGgc                                                                                                                                              >  2:411289/1‑150 (MQ=255)
      tAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGgcg                                                                                                                                             >  2:1012136/1‑150 (MQ=255)
      tAAAGTGGTATTTTACATGCACTTACAATAGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGTTGCAATGgcg                                                                                                                                             <  2:847415/150‑1 (MQ=255)
        aaGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGTTGCAATGgcggt                                                                                                                                           <  2:326014/150‑3 (MQ=255)
         aGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAAACAATAAAAACCGCCATTGTTGCAATGgcggt                                                                                                                                           <  2:547218/150‑3 (MQ=255)
              tATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCgg                                                                                                                                     <  1:109347/150‑1 (MQ=255)
               aTTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGa                                                                                                                                    >  2:1060674/1‑150 (MQ=255)
                 tttACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGAtt                                                                                                                                  >  1:284919/1‑149 (MQ=255)
                 tttACATGCAATTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGGTACTGTCGTG‑AAAAACAATAAAAACCGCCATTTCAACAATGGCGCAATTCGGATg                                                                                                                                  <  2:505281/150‑1 (MQ=255)
                  ttACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCACCATTGCAACAATGGCGCAATTCGGATGa                                                                                                                                 >  1:456268/1‑150 (MQ=255)
                     cATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTc                                                                                                                                       >  2:74220/1‑141 (MQ=255)
                     cATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTc                                                                                                                                       <  1:74220/141‑1 (MQ=255)
                        gCACTTACAATTGATTAAAGACAACATTCACAGTGTGGGTATTTGTTACACATAGGGGCGAGCAATGTCATGCCAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACTGCCATTGCAACAATGGCGCAATTCGGATGAAGcccc                                                                                                                           <  2:851352/150‑1 (MQ=255)
                         cACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGTGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGAAATTCGGATGAAGCCCCt                                                                                                                          >  1:573841/1‑150 (MQ=255)
                            ttACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATg                                                                                                                       <  1:993377/150‑1 (MQ=255)
                               cAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGTGCGAGCAAAGTCATGACACTGTAGGTGAGGTTACTGTCTTT‑AAAAACAATAAAAACCTCCATTGCAACAATGGCGCAATTCCGATGAAGCCCATATGACa                                                                                                                    >  1:208970/1‑150 (MQ=255)
                                aaTTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCGATGACaa                                                                                                                   >  2:767797/1‑150 (MQ=255)
                                 aTTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACaaa                                                                                                                  >  1:820495/1‑149 (MQ=255)
                                 aTTGAATAAAGACAACATTCACAGTGTGCTAATTTGTTACACATACGGGCGAGCAAAGTCATGACAGTGTAGGTGCTGTTACTGTCGTG‑AAAAACAAAAAAAACCGCCATAGCAACAATGGCGCAATTCGTATGAATCACCTATTAAAAg                                                                                                                  >  1:810716/1‑150 (MQ=255)
                                  ttGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAgg                                                                                                                 <  1:863452/150‑1 (MQ=255)
                                    gATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGAGGTTACTTTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAATGAt                                                                                                               >  1:58094/1‑150 (MQ=255)
                                     aTTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATa                                                                                                              >  2:605612/1‑150 (MQ=255)
                                      ttAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaa                                                                                                             >  2:712908/1‑150 (MQ=255)
                                       tAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaa                                                                                                            >  1:582692/1‑150 (MQ=255)
                                       tAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGGCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaa                                                                                                            <  1:320393/150‑1 (MQ=255)
                                        aaaGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCATTTCGGATGAAGCCCGTATGACAAGGTTaaaa                                                                                                           >  2:184851/1‑150 (MQ=255)
                                           gACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATa                                                                                                        >  1:496530/1‑150 (MQ=255)
                                              aacaTTCACAGTGTGGTTATTTGTTACACATAGTGGCGATCAATGTCATGACATTGTAGGTTCGGTTACTGTAGTG‑AAAAACAATAAAAACCTCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAac                                                                                                     <  1:246798/150‑1 (MQ=255)
                                                        gtgtGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTa                                                                                           >  1:852385/1‑150 (MQ=255)
                                                        gtgtGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTa                                                                                           >  2:137941/1‑150 (MQ=255)
                                                        gtgtGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTa                                                                                           >  1:322150/1‑150 (MQ=255)
                                                          gtgGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATc                                                                                         <  1:1116427/150‑1 (MQ=255)
                                                            ggTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGAGTACGTGCGGTTACTGTCGTG‑AAAAACAAAAAAAACCGCCATTGCAACAATGGCGCTATTCGGATGAATCCCCTATGACAAGGAAAAAAATAAACGCACGCCGTATAtt                                                                                       >  2:1172713/1‑150 (MQ=255)
                                                              ttATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCa                                                                                     <  2:463174/150‑1 (MQ=255)
                                                                 tttGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCt                                                                                  <  1:229819/150‑1 (MQ=255)
                                                                  ttGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCtt                                                                                 <  1:816656/150‑1 (MQ=255)
                                                                     ttACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAtt                                                                              <  1:712908/150‑1 (MQ=255)
                                                                     ttACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAAATCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCGTCCGCTTAtt                                                                              >  2:578653/1‑150 (MQ=255)
                                                                      tACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTg                                                                             >  1:573135/1‑150 (MQ=255)
                                                                           aTAGGGGCGAGCAATGTCATGACAGTTTTGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACAGCCATTGCAACAATGTCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTcc                                                                        <  2:922833/150‑1 (MQ=255)
                                                                               gggCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtttt                                                                    <  1:605612/150‑1 (MQ=255)
                                                                               gggCGAGCAATGTCATGACAGTGTAGGTGCGGATACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAATCCCCTTTGACAATGATAAAAATAAACGCACTCCGTATCTTCATCTTATTGATTcgattt                                                                    >  2:71511/1‑145 (MQ=255)
                                                                                ggCGAGCAAGGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGTTTAAAATAAACGCACGCCGTTTCTTCAGCTTATTGTTTCGttttt                                                                   >  1:625857/1‑150 (MQ=255)
                                                                                 gCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtttttt                                                                  <  1:1146452/150‑1 (MQ=255)
                                                                                   gAGCAATGTCATGACAGTGTAGGTGCGGTTTCTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGGGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttct                                                                <  1:1133121/150‑1 (MQ=255)
                                                                                     gCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttt                                                              <  1:994379/150‑1 (MQ=255)
                                                                                       aaTGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATTAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt                                                            >  2:689553/1‑150 (MQ=255)
                                                                                       aaTGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAACCGCCCGCCGTATCTTCAGCTTATTGATTCCTtgtttcttttt                                                            >  2:170604/1‑150 (MQ=255)
                                                                                       aaTGTCATGACAGTGTAGGTGAGGTTACTGTCGTG‑AAAAAAAATAAAAACCGCCATTGCAACAATGGAGCAATTCGGATGAAGCCCCTAAGACAAGGATAAAAATAAACGCACGCAGGATCTTCAGCTTATTGATTCCTtgtttcttttt                                                            >  2:166900/1‑150 (MQ=255)
                                                                                          gTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCAc                                                         >  2:351646/1‑150 (MQ=255)
                                                                                          gTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCAc                                                         >  1:733090/1‑150 (MQ=255)
                                                                                          gTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCAc                                                         >  1:386860/1‑150 (MQ=255)
                                                                                            cATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACtt                                                       >  2:553026/1‑150 (MQ=255)
                                                                                              tGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGCAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGGATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCt                                                     >  1:348745/1‑150 (MQ=255)
                                                                                              tGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCCGCTTATTGATTCCTTTTTTCTTTTTCCATTCt                                                     >  2:124374/1‑150 (MQ=255)
                                                                                              tGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAATGAGTCCTTTTTTCTTTTTCACTTCt                                                     >  1:892309/1‑150 (MQ=255)
                                                                                                    tgtAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCATTTTCACTTCTGTTCAc                                               >  2:150773/1‑150 (MQ=255)
                                                                                                       tggTGCGGTTACTGTAGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCt                                            <  2:573841/149‑1 (MQ=255)
                                                                                                       aGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATGCCTTTTTTATTTTTAACTTCTGTTCACGCt                                            >  2:243603/1‑150 (MQ=255)
                                                                                                        ggTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTg                                           >  1:710170/1‑150 (MQ=255)
                                                                                                          tGCTGTTACTTTCGTG‑AAAAACAATATAATCTGCCTTTGATTCATTGGCGCAATTCGGATGAAGCCCCTATGACAAGGATTAAAATATAAGCACGCCGTATCTTCAGGTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGaa                                         <  2:810716/150‑1 (MQ=255)
                                                                                                                  cTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgc                                 <  2:436679/150‑1 (MQ=255)
                                                                                                                  cTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgc                                 <  2:1032415/150‑1 (MQ=255)
                                                                                                                  cTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgc                                 >  2:1239184/1‑150 (MQ=255)
                                                                                                                  cTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgc                                 <  1:1315567/150‑1 (MQ=255)
                                                                                                                    gTCGTG‑AAAAACAATTAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAACCGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctg                               <  1:184192/150‑1 (MQ=255)
                                                                                                                    gTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctg                               <  1:711286/150‑1 (MQ=255)
                                                                                                                     gcGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctgc                              <  2:519689/149‑1 (MQ=255)
                                                                                                                       gTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTc                            >  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)
                                                                                                                                                    |                                                                                                                                                   
TAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCG  >  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: ‑