Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,286,793 A→G 10.2% K29E (AAG→GAG)  ychS → putative uncharacterized protein YchS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,286,7930AG10.2% 200.5 / 11.9 80K29E (AAG→GAG) ychSputative uncharacterized protein YchS
Reads supporting (aligned to +/- strand):  ref base A (21/49);  new base G (3/5);  total (26/54)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTC  >  NC_000913/1286644‑1286925
                                                                                                                                                     |                                                                                                                                    
ggTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGACGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCa                                                                                                                                      <  1:1870729/150‑1 (MQ=255)
ggTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCa                                                                                                                                      <  1:1491440/150‑1 (MQ=255)
  tcaaCGTGGGGCCATCGGGCAAAATTCGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGGGGGGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAg                                                                                                                                    <  2:1943041/148‑1 (MQ=255)
  tAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCAAg                                                                                                                                    >  1:2396029/1‑150 (MQ=255)
      cGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCCAGGCGa                                                                                                                                >  1:563017/1‑150 (MQ=255)
      cGTGGTGCCATATGCCAAAATTAGGTCGGAGTAACAGTTTTTGATTCTTTAAGTAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGa                                                                                                                                <  2:790345/150‑1 (MQ=17)
        tggtgCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAAc                                                                                                                              >  2:283633/1‑150 (MQ=255)
         ggtgCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAAcc                                                                                                                             <  2:1649928/150‑1 (MQ=255)
          gtgCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCt                                                                                                                            <  1:217646/150‑1 (MQ=255)
            gCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTg                                                                                                                          >  1:1526465/1‑150 (MQ=255)
            gCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTg                                                                                                                          >  2:2276597/1‑150 (MQ=255)
                tCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCg                                                                                                                      <  1:1039276/150‑1 (MQ=255)
                   ggCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAg                                                                                                                   <  1:2408116/150‑1 (MQ=255)
                   ggCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAg                                                                                                                   <  1:1367636/150‑1 (MQ=255)
                    gCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGc                                                                                                                  <  2:1356294/150‑1 (MQ=255)
                       aaaTTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTATTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTtc                                                                                                               <  2:1097977/150‑1 (MQ=255)
                       aaaTTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTtc                                                                                                               >  1:458686/1‑150 (MQ=255)
                         aTTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACCCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTtctc                                                                                                             <  1:1491067/150‑1 (MQ=255)
                             gCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATcc                                                                                                         <  1:676288/150‑1 (MQ=255)
                               tCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCAAGACGAACCTTAGTCGAAGCTTCTCATCCtt                                                                                                       >  2:435403/1‑150 (MQ=255)
                                 ggAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTcc                                                                                                     <  2:472484/150‑1 (MQ=255)
                                  gAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGACGAACCTTAGTCGAAGCTTCTCATCCTTccc                                                                                                    >  2:407337/1‑150 (MQ=255)
                                     tAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGc                                                                                                 <  2:1174722/150‑1 (MQ=255)
                                       aCAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCtt                                                                                               <  2:1763258/150‑1 (MQ=255)
                                          ggTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTggg                                                                                            <  2:686217/150‑1 (MQ=255)
                                           gTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCGCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGc                                                                                           >  1:1561189/1‑150 (MQ=255)
                                           gTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGc                                                                                           <  1:238374/150‑1 (MQ=255)
                                           gTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGc                                                                                           <  1:407337/150‑1 (MQ=11)
                                            ttttGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAt                                                                                                           >  1:1383533/1‑133 (MQ=255)
                                            ttttGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAt                                                                                                           <  2:1383533/133‑1 (MQ=255)
                                             tttGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGca                                                                                                >  2:468375/1‑142 (MQ=255)
                                             tttGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGca                                                                                                <  1:468375/143‑2 (MQ=255)
                                              ttGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGa                                                                                        >  1:456935/1‑150 (MQ=255)
                                                tattattTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAat                                                                                      <  2:848687/149‑1 (MQ=255)
                                                gATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAat                                                                                      <  2:1676753/150‑1 (MQ=255)
                                                 attattTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATc                                                                                                          >  1:1037116/1‑129 (MQ=255)
                                                 attattTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATc                                                                                                          <  2:1037116/129‑1 (MQ=255)
                                                    attTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATAttt                                                                                  <  2:118237/150‑1 (MQ=255)
                                                     tttAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTg                                                                                 <  1:1888677/150‑1 (MQ=255)
                                                       tAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGAt                                                                               <  1:890838/150‑1 (MQ=255)
                                                         aGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCGTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTg                                                                             <  2:458686/150‑1 (MQ=255)
                                                           gAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCAGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCg                                                                           <  2:752226/150‑1 (MQ=255)
                                                                   ggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGtt                                                                   <  1:283633/150‑1 (MQ=37)
                                                                    gtggtggtgGAATGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCTGCAACGCTCTCTCGCTGGCTCTCACGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGttt                                                                  <  2:563746/150‑1 (MQ=11)
                                                                    gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTCGATTGCGGATTCGttt                                                                  <  2:1671325/150‑1 (MQ=38)
                                                                       gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTGGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTgag                                                               <  1:82583/150‑1 (MQ=255)
                                                                              aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCg                                                        <  2:2071831/150‑1 (MQ=255)
                                                                               aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCgg                                                       <  2:1824949/150‑1 (MQ=255)
                                                                                ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCggg                                                      >  2:508924/1‑150 (MQ=255)
                                                                                ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCggg                                                      <  2:1022178/150‑1 (MQ=255)
                                                                                  aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGc                                                    <  1:2276597/150‑1 (MQ=255)
                                                                                     aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCCAGTCGAACCTTGGTCGAAGCTTCTCCTCCTTCCCCGCTTGGGCAGAAGATTTGATTGCGGATTCGTTTGAGAGTTCAGGGGCgtt                                                 >  1:2185136/1‑150 (MQ=11)
                                                                                     aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCttt                                                 <  2:737692/150‑1 (MQ=255)
                                                                                      cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGGTTGAGAGTTCCGGGGCtttt                                                >  1:1589547/1‑150 (MQ=255)
                                                                                       tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTg                                               <  2:2284165/150‑1 (MQ=255)
                                                                                       tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTg                                               <  2:1754463/150‑1 (MQ=255)
                                                                                        cAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTAGTAGAAGCTTCTCATCATTACCCGCTTGGGCAGAATATTTGATTGAGGATTCGTTTGAGAATTCCGGGGATTTTGa                                              >  2:53464/1‑150 (MQ=14)
                                                                                         aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGaa                                             <  2:1476402/150‑1 (MQ=255)
                                                                                           cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAg                                           <  2:2265029/150‑1 (MQ=255)
                                                                                           cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAg                                           <  2:1454350/150‑1 (MQ=255)
                                                                                              tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGa                                        <  1:462928/150‑1 (MQ=255)
                                                                                              tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGa                                        <  1:508924/150‑1 (MQ=255)
                                                                                                 gcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGAtgg                                     <  2:456935/150‑1 (MQ=255)
                                                                                                       cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGAtggtggtgg                               <  2:750266/150‑1 (MQ=255)
                                                                                                        gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTTATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTggg                              <  2:2239441/150‑1 (MQ=255)
                                                                                                        gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTggg                              <  2:2313283/150‑1 (MQ=255)
                                                                                                         cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCAATTCTAACTTTGGTCGAATCTTCTCATCCTTCCTCTCTTGTGCATAATATTTGATTGCGGTTTCGTTTTATAATTCGTTTTATTTTGTAAGTGATGGTGGTgggg                             <  2:1561189/150‑1 (MQ=14)
                                                                                                         cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTgggg                             >  2:1154769/1‑150 (MQ=255)
                                                                                                          ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGGTTGAGACTTCAGGGGCTTTTGAAAGTGATGGTGGTggggg                            >  1:1510448/1‑150 (MQ=255)
                                                                                                              gagagcGTTGCCGCAGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAgg                        >  1:520677/6‑150 (MQ=255)
                                                                                                               gggTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAATTCGAACCTTGGTCGAATCTTCTCATCCTTCCCCGCTTGGGCATAATCTTTCATTCCGGATGCGTTTGAGAATTCGCGGCCTTTTGAAAGTGATGGTGGTGGGGGAAGGa                       <  2:22021/150‑1 (MQ=14)
                                                                                                               gggTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGa                       <  1:17869/150‑1 (MQ=255)
                                                                                                                 gtctccgCCTCCGCCAACGCTCTCTCGCTGGCTCTCACCTCTAACCTTGTTCTAATTTTATGATCGTTCCACGGTTGGGCATATTTTTTAATTTGTTATTAGTTTGAGATTTCGTTTTTTTATGAAAGTGATGGTGGTGGGGGAAGGAtt                     <  2:338711/144‑1 (MQ=9)
                                                                                                                 gTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGAtt                     <  2:589539/150‑1 (MQ=255)
                                                                                                                 agcGTTGCCGCAGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGAtt                     >  2:463554/3‑150 (MQ=255)
                                                                                                                    gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACt                  <  1:1467321/150‑1 (MQ=255)
                                                                                                                       gCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAg               >  2:1441150/1‑150 (MQ=255)
                                                                                                                               cAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGGTTCGTTTGAGAAGTCCGGGGCTGTGGAGTGTGAGGGTGGTGGGGGAAGGATTACTCAGCGGTgaga       >  1:15832/1‑147 (MQ=17)
                                                                                                                                    ctctctCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCGGAATATTTGATTGCGGATTCGTTTGAGAGTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTc  >  1:2201821/1‑150 (MQ=255)
                                                                                                                                    ctctctCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTgggg                             <  1:1747112/123‑1 (MQ=255)
                                                                                                                                    ctctctCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTgggg                             >  2:1747112/1‑123 (MQ=255)
                                                                                                                                                  cTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAata                                                                                     >  1:1573317/1‑53 (MQ=255)
                                                                                                                                                     |                                                                                                                                    
GGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTC  >  NC_000913/1286644‑1286925

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