Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 177,396 C→T 6.0% V77M (GTG→ATG)  yadS ← conserved inner membrane protein YadS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913177,3960CT 6.0% 254.0 / 2.9 81V77M (GTG→ATG) yadSconserved inner membrane protein YadS
Reads supporting (aligned to +/- strand):  ref base C (44/32);  new base T (3/2);  total (47/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GACGCCCATACAAACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTCC  >  NC_000913/177253‑177544
                                                                                                                                               |                                                                                                                                                    
gACGCCCATACAAACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGt                                                                                                                                                >  2:6877/1‑150 (MQ=255)
       aTACAAACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGCCCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATg                                                                                                                                         >  2:653462/1‑150 (MQ=255)
        tACAAACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGc                                                                                                                                        >  2:114768/1‑150 (MQ=255)
           aaaCCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTgg                                                                                                                                     <  1:280305/150‑1 (MQ=255)
           aaaCCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTgg                                                                                                                                     <  1:1665354/150‑1 (MQ=255)
            aaCCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGt                                                                                                                                    <  1:785312/150‑1 (MQ=255)
              ccGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGa                                                                                                                                  <  2:793144/150‑1 (MQ=255)
               cgcgATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCAGCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGAc                                                                                                                                 <  2:982948/150‑1 (MQ=255)
               cgcgATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCATGCTGACCATCGTGCTGGTgcg                                                                                                                                 <  1:416694/150‑3 (MQ=255)
                  gATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCAt                                                                                                                              >  1:1165911/1‑150 (MQ=255)
                   aTTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCAtt                                                                                                                             >  1:820986/1‑150 (MQ=255)
                   aTTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCACCAGCATGCTGACCAtc                                                                                                                             <  1:454256/150‑2 (MQ=255)
                     tAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGc                                                                                                                           <  1:1509844/150‑1 (MQ=255)
                     tAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGc                                                                                                                           >  1:1630085/1‑150 (MQ=255)
                      aaCGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCa                                                                                                                          >  2:1407434/1‑150 (MQ=255)
                       aCGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCaa                                                                                                                         >  2:924034/1‑150 (MQ=255)
                        cGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAAc                                                                                                                        >  2:838333/1‑150 (MQ=255)
                        cGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTAGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATAGCAAc                                                                                                                        >  2:920403/1‑150 (MQ=255)
                          gACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCGTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGa                                                                                                                      <  1:147531/150‑1 (MQ=255)
                            ccGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGAcc                                                                                                                    >  1:1637310/1‑150 (MQ=255)
                             cGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCa                                                                                                                   <  2:826651/150‑1 (MQ=255)
                              ggCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAg                                                                                                                  <  2:736221/150‑1 (MQ=255)
                              ggCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAg                                                                                                                  <  1:114768/150‑1 (MQ=255)
                                   ccGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCCGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCgg                                                                                                             >  2:1652018/1‑150 (MQ=255)
                                    cGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGt                                                                                                            >  2:841169/1‑150 (MQ=255)
                                     gCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTg                                                                                                           >  1:931428/1‑150 (MQ=255)
                                     gCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTg                                                                                                           >  1:185673/1‑150 (MQ=255)
                                          aaaGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATc                                                                                                      >  1:1021722/1‑150 (MQ=255)
                                           aaGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCt                                                                                                     <  2:1409034/150‑1 (MQ=255)
                                               cTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCa                                                                                                 >  2:398731/1‑150 (MQ=255)
                                               cTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCa                                                                                                 >  1:379464/1‑150 (MQ=255)
                                                  tATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCAccc                                                                                              >  1:966234/1‑150 (MQ=255)
                                                   aTTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCa                                                                                             <  1:1355606/150‑1 (MQ=255)
                                                      cACGCCAATGCCGACAAACACCGCCAGCCCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCCCGCTGGTCAGCATGCTGGTGCCCATTGCAACGCCCAGCTCGGTGGGATCTTTCACTCaaaa                                                                                          >  2:927094/1‑150 (MQ=255)
                                                      cACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTg                                                                                                                                   >  2:1456673/1‑109 (MQ=255)
                                                      cACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTg                                                                                                                                   <  1:1456673/109‑1 (MQ=255)
                                                      cACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCaaaa                                                                                          >  1:503213/1‑150 (MQ=255)
                                                      cACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCaaaa                                                                                          >  2:1139103/1‑150 (MQ=255)
                                                            aaTGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCgg                                                                                    <  1:119990/150‑1 (MQ=255)
                                                                 cGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGt                                                                               >  1:1613394/1‑150 (MQ=255)
                                                                   aCAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGa                                                                             <  2:1054577/150‑1 (MQ=255)
                                                                    cAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGAt                                                                            <  1:1407434/150‑1 (MQ=255)
                                                                      aaCACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATcc                                                                          >  2:1594725/1‑150 (MQ=255)
                                                                      aaCACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATcc                                                                          <  1:838333/150‑1 (MQ=255)
                                                                      aaCACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATcc                                                                          >  2:506517/1‑150 (MQ=255)
                                                                       acacCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCa                                                                         <  2:820986/150‑1 (MQ=255)
                                                                        cacCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAg                                                                        >  1:120087/1‑150 (MQ=255)
                                                                            gCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCCCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAGATACCGGGCCGTGATCCAGCGcc                                                                    >  1:55460/1‑150 (MQ=255)
                                                                                  ccAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGa                                                                                                        >  2:1061400/1‑108 (MQ=255)
                                                                                  ccAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGa                                                                                                        <  1:1061400/108‑1 (MQ=255)
                                                                                  ccAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTcg                                                              >  1:926488/1‑150 (MQ=255)
                                                                                    aaCGGCGTCCAGCACCGGCACCATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGGGGGATCATTCACCCAAAATACCGGGCCGTGATCCCGCGCCATGTcgcg                                                            >  2:1620558/1‑150 (MQ=255)
                                                                                      cGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACTACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGaa                                                          <  2:452801/150‑1 (MQ=255)
                                                                                         cGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTg                                                       <  2:1637310/150‑1 (MQ=255)
                                                                                           tCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTc                                                     >  2:1666214/1‑150 (MQ=255)
                                                                                           tCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTc                                                     <  1:42763/150‑1 (MQ=255)
                                                                                            ccAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTcc                                                    <  2:953919/150‑1 (MQ=255)
                                                                                             cAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTccc                                                   >  2:40849/1‑150 (MQ=255)
                                                                                              aGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCccg                                                  <  2:926488/150‑1 (MQ=255)
                                                                                                 gccGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCccgccg                                               >  1:434258/2‑150 (MQ=255)
                                                                                                 aCCGGCAAGATCCTTTTTGGTAAGCGTCTTGGCTGGCGCTCGAGCATGATGGACTGCATGCTGGTGACGATTGCAACGACCAGATCGGTGGGATCTTTGACGAAAAATTCCGGGCCGTGTTCCATCGCCATTTCGCGAATTGTCccgccg                                               <  2:228488/150‑1 (MQ=255)
                                                                                                    ggCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCt                                            >  1:1038592/1‑150 (MQ=255)
                                                                                                         cATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGc                                       >  1:512947/1‑150 (MQ=255)
                                                                                                           tCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCgg                                     >  2:679990/1‑150 (MQ=255)
                                                                                                            ccATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGt                                    >  2:489960/1‑150 (MQ=255)
                                                                                                            ccATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGt                                    >  2:544360/1‑150 (MQ=255)
                                                                                                                   ggTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCaag                             >  2:1084409/1‑148 (MQ=255)
                                                                                                                      aaGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGccc                          <  1:544360/150‑1 (MQ=255)
                                                                                                                         cGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGt                       >  1:1004062/1‑150 (MQ=255)
                                                                                                                            cTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTAcc                    <  2:379464/150‑1 (MQ=255)
                                                                                                                              tggctggcGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAg                  >  1:33008/1‑150 (MQ=255)
                                                                                                                              tggctggcGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAg                  <  2:1630085/150‑1 (MQ=255)
                                                                                                                               ggctggcGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGa                 <  1:924034/150‑1 (MQ=255)
                                                                                                                                    ggtcACCAGCATGCTGACCATCGTGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACAcc            >  2:445253/4‑150 (MQ=255)
                                                                                                                                     gtcACCAGCATGCTGACCATCGTGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCa           >  1:1331705/3‑150 (MQ=255)
                                                                                                                                     gtcACCAGCATGCTGACCATCGTGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCa           >  1:167618/3‑150 (MQ=255)
                                                                                                                                           aGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAggg     >  1:553190/1‑150 (MQ=255)
                                                                                                                                           aGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAggg     <  1:489960/150‑1 (MQ=255)
                                                                                                                                            gCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGt    >  1:229496/1‑150 (MQ=255)
                                                                                                                                             cACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTc   <  1:653462/150‑1 (MQ=255)
                                                                                                                                              aCGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTcc  <  2:973838/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GACGCCCATACAAACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTCC  >  NC_000913/177253‑177544

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: