Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,287,734 C→T 7.2% intergenic (‑363/‑30) 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,7340CT 7.2% 245.8 / 6.4 84intergenic (‑363/‑30)acs/nrfAacetyl‑CoA synthetase/nitrite reductase, formate‑dependent, cytochrome
Reads supporting (aligned to +/- strand):  ref base C (47/30);  new base T (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 = 6.86e-01

GACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCG  >  NC_000913/4287586‑4287883
                                                                                                                                                     |                                                                                                                                                     
gACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGca                                                                                                                                                      >  1:463174/1‑150 (MQ=255)
   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)
                                                                                                                                                     |                                                                                                                                                     
GACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTG‑AAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCG  >  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: ‑