Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,917,896 C→A 6.1% intergenic (+83/+149) barA → / ← gudD sensory histidine kinase BarA/D‑glucarate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,917,8960CA 6.1% 321.2 / 6.4 96intergenic (+83/+149)barA/gudDsensory histidine kinase BarA/D‑glucarate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (44/46);  new base A (4/2);  total (48/48)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.11e-01

Only 100 of 102 total aligned reads displayed.

AGCTGCTGGAACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCT  >  NC_000913/2917749‑2918045
                                                                                                                                                   |                                                                                                                                                     
aGCTGCTGGAACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATgcgg                                                                                                                                                     <  2:238438/150‑1 (MQ=255)
  gtgctgGTTCTGTTGGTCGAGATGGATAATGTCGAGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATgcggcg                                                                                                                                                   <  2:857652/149‑1 (MQ=255)
  ctgctgGAACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATgcggcg                                                                                                                                                   <  1:1575970/150‑1 (MQ=255)
   tgctgGAACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGt                                                                                                                                                  >  1:283791/1‑150 (MQ=255)
    gctgGAACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGtt                                                                                                                                                 >  1:1627971/1‑150 (MQ=255)
     ctgGAACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTa                                                                                                                                                >  1:1652436/1‑150 (MQ=255)
      tgGAACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTaa                                                                                                                                               <  1:1446870/150‑1 (MQ=255)
        gAACTGGTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACg                                                                                                                                             <  1:563910/150‑1 (MQ=255)
          aCTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                                                           <  1:849493/150‑1 (MQ=255)
             gTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                                                        <  1:896340/150‑1 (MQ=255)
               tggaatagaTGGATAATGTCGCTAGCGAAGACAGCAAAATTATCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGCTTGGTGACAAAGTGACTCTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                                                      >  2:1429781/7‑150 (MQ=38)
                 gACGAGATGGTTAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                                                    >  1:1881011/1‑150 (MQ=255)
                   cGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                                                  >  1:663178/1‑150 (MQ=255)
                   cGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                                                  >  1:1862517/1‑150 (MQ=255)
                   cGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                                                  <  2:983462/150‑1 (MQ=255)
                    gagaTGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                                                 <  2:618094/150‑1 (MQ=255)
                       aTGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGGGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                                                                              <  2:1062977/150‑1 (MQ=255)
                           aTAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaaa                                                                                                                          <  1:122766/150‑1 (MQ=255)
                            tAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAc                                                                                                                         <  1:738820/150‑1 (MQ=255)
                              aTGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCa                                                                                                                       >  1:68168/1‑150 (MQ=255)
                                gTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAAt                                                                                                                     <  2:283791/150‑1 (MQ=255)
                                gTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAAt                                                                                                                     <  2:1981326/150‑1 (MQ=255)
                                  cgcgcgcgAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTa                                                                                                                   >  1:583241/1‑150 (MQ=255)
                                     gcgcgAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAAt                                                                                                                <  1:1998012/150‑1 (MQ=255)
                                        cgAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCa                                                                                                             >  2:954925/1‑150 (MQ=255)
                                        cgAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCa                                                                                                             >  2:498101/1‑150 (MQ=255)
                                          aaGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAAt                                                                                                           <  2:1150089/150‑1 (MQ=255)
                                          aaGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAAt                                                                                                           >  1:747654/1‑150 (MQ=255)
                                             ccAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAAtt                                                                                                                    <  2:1105428/138‑1 (MQ=255)
                                             ccAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAAtt                                                                                                                    >  1:1105428/1‑138 (MQ=255)
                                                gCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTa                                                                                                     >  2:300263/1‑150 (MQ=255)
                                                   aaaTTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGTATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTaaaa                                                                                                  <  1:1225068/150‑1 (MQ=255)
                                                    aaTTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTaaaaa                                                                                                 <  1:1525503/150‑1 (MQ=255)
                                                        ctcGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAatga                                                                                             <  1:831356/150‑1 (MQ=255)
                                                         tcGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAatgat                                                                                            <  2:1627971/150‑1 (MQ=255)
                                                          cGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAatgatg                                                                                           <  2:1362035/150‑1 (MQ=255)
                                                          cGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAatgatg                                                                                           <  2:1654961/150‑1 (MQ=255)
                                                           gggTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTaaa                                                                                                                 <  1:696213/127‑1 (MQ=255)
                                                           gggTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTaaa                                                                                                                 >  2:696213/1‑127 (MQ=255)
                                                                   gATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGGGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGa                                                                                  >  2:213594/1‑150 (MQ=255)
                                                                   gATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGa                                                                                  >  1:454986/1‑150 (MQ=255)
                                                                    aTGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGAt                                                                                 <  2:768172/150‑1 (MQ=255)
                                                                     tGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATa                                                                                >  1:677589/1‑150 (MQ=255)
                                                                     tGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATa                                                                                >  1:1635605/1‑150 (MQ=255)
                                                                      gTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATAc                                                                               <  2:454986/150‑1 (MQ=255)
                                                                       tCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACa                                                                              >  2:954139/1‑150 (MQ=255)
                                                                       tCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACa                                                                              >  1:319614/1‑150 (MQ=255)
                                                                        cTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACAt                                                                             <  2:663178/150‑1 (MQ=255)
                                                                        cTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACAt                                                                             <  2:1153310/150‑1 (MQ=255)
                                                                         ttGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATa                                                                            <  2:68168/150‑1 (MQ=255)
                                                                          tGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAg                                                                           <  2:300948/150‑1 (MQ=255)
                                                                           gTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGc                                                                          <  2:1881011/150‑1 (MQ=255)
                                                                              gAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGcaca                                                                       >  1:1699126/1‑150 (MQ=255)
                                                                              gAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGcaca                                                                       >  2:883347/1‑150 (MQ=255)
                                                                              gAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGcaca                                                                       >  2:944428/1‑150 (MQ=255)
                                                                               aGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCacaa                                                                      <  2:747654/150‑1 (MQ=255)
                                                                                  aGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAg                                                                   >  2:1533787/1‑150 (MQ=255)
                                                                                   gCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAgg                                                                  >  2:227585/1‑150 (MQ=255)
                                                                                   gCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAgg                                                                  >  1:1402715/1‑150 (MQ=255)
                                                                                    cAATCGTCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGc                                                                 <  2:855147/150‑1 (MQ=255)
                                                                                    cAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGc                                                                 <  1:1278137/150‑1 (MQ=255)
                                                                                      aTCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTa                                                               <  1:300263/150‑1 (MQ=255)
                                                                                        cGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTAtt                                                             >  2:1501851/1‑150 (MQ=255)
                                                                                         ggCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTAttt                                                            >  1:1918966/1‑150 (MQ=255)
                                                                                           cTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTg                                                          <  1:1632363/150‑1 (MQ=255)
                                                                                             ttttATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCg                                                        >  1:1317856/1‑150 (MQ=255)
                                                                                              tttATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGt                                                       >  2:1713725/1‑150 (MQ=255)
                                                                                                 aTGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTa                                                    <  2:1652436/150‑1 (MQ=255)
                                                                                                   gACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGc                                                  >  2:1236341/1‑150 (MQ=255)
                                                                                                    aCGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCa                                                 >  2:1145260/1‑150 (MQ=255)
                                                                                                     cGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCAt                                                <  1:339032/150‑1 (MQ=255)
                                                                                                      gCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATc                                               <  2:310736/150‑1 (MQ=255)
                                                                                                       cTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCGTAAAAAACAAATTAAATTAAATGAATTAAAAAATAATTTGGCCGTGAAAAATGGCACAACAGCGTATTTTGGGTTAGCAttca                                             >  1:162142/1‑147 (MQ=21)
                                                                                                          gTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGtc                                           >  2:842400/1‑150 (MQ=255)
                                                                                                          gTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGta                                           >  2:1785180/1‑149 (MQ=255)
                                                                                                            ttCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGtctc                                         >  1:218007/1‑150 (MQ=255)
                                                                                                                  ttGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCg                                   <  1:227585/150‑1 (MQ=255)
                                                                                                                    ggTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGc                                 <  2:319614/150‑1 (MQ=255)
                                                                                                                     gTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCt                                >  1:578351/1‑150 (MQ=255)
                                                                                                                       gACAAAGTGCCTGTGTTTATGCCGGATGCGGGGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTcc                              <  1:213594/150‑1 (MQ=255)
                                                                                                                       gACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTcc                              >  2:349708/1‑150 (MQ=255)
                                                                                                                            aGTGCCTGTGTTTATGCCGGATGCGTCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCTGCTCCAGATa                         <  2:203656/150‑1 (MQ=255)
                                                                                                                             gTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATag                        >  1:1822907/1‑150 (MQ=255)
                                                                                                                                ccTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGc                     >  2:604806/1‑150 (MQ=255)
                                                                                                                                ccTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGc                     >  1:232402/1‑150 (MQ=255)
                                                                                                                                       ttATGCCGGATGAGTCGTTAACTCATTATCCTGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggtttt              <  1:1975229/150‑1 (MQ=255)
                                                                                                                                         aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggtttgag            <  2:583241/150‑4 (MQ=38)
                                                                                                                                         aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttgg            <  1:1533787/150‑1 (MQ=255)
                                                                                                                                          ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggt           >  1:1943023/2‑150 (MQ=255)
                                                                                                                                          ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggt           >  1:1968972/2‑150 (MQ=255)
                                                                                                                                           gCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACTTAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggtt          >  2:1790162/1‑150 (MQ=255)
                                                                                                                                            ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggttt         <  2:1292646/150‑1 (MQ=255)
                                                                                                                                            ccGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggttt         >  1:1814408/1‑150 (MQ=255)
                                                                                                                                             cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggtttga             >  2:395207/1‑143 (MQ=255)
                                                                                                                                             cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggtttga             <  1:395207/145‑3 (MQ=255)
                                                                                                                                              ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAAACAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTc       <  1:1889285/150‑1 (MQ=255)
                                                                                                                                               gATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCt      >  1:1907232/1‑150 (MQ=255)
                                                                                                                                                aTGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTg     <  2:218007/150‑1 (MQ=255)
                                                                                                                                                  gcggcgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTc   <  1:349708/150‑1 (MQ=255)
                                                                                                                                                  gcggcgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTc   <  1:954925/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGCTGCTGGAACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCT  >  NC_000913/2917749‑2918045

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: