Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,313,582 A→T 5.2% K78M (AAG→ATG) ‡ rpiB → allose‑6‑phosphate isomerase/ribose‑5‑phosphate isomerase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,313,5820AT 5.2% 439.1 / 9.6 134K78M (AAG→ATG) ‡rpiBallose‑6‑phosphate isomerase/ribose‑5‑phosphate isomerase B
Reads supporting (aligned to +/- strand):  ref base A (63/64);  new base T (3/4);  total (66/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01

Only 100 of 137 total aligned reads displayed.

TGAAGTGATTGATAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCAC  >  NC_000913/4313433‑4313731
                                                                                                                                                     |                                                                                                                                                     
tGAAGTGATTGATAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACaa                                                                                                                                                       >  2:323359/1‑150 (MQ=255)
     tgatttatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGTGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTg                                                                                                                                                  <  2:1250014/150‑1 (MQ=255)
     tgattgatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTg                                                                                                                                                  <  1:930959/150‑1 (MQ=255)
     tgattgatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTg                                                                                                                                                  <  1:935686/150‑1 (MQ=255)
      gattgatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGc                                                                                                                                                 >  1:581829/1‑150 (MQ=255)
      gattgatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGc                                                                                                                                                 >  1:151359/1‑150 (MQ=255)
        ttgatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCg                                                                                                                                               >  1:1330829/1‑150 (MQ=255)
        ttgatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCAAACTTGTTCGCCg                                                                                                                                               >  2:886781/1‑150 (MQ=255)
         tgatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCgg                                                                                                                                              <  1:1087974/150‑1 (MQ=255)
          gatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGa                                                                                                                                             >  1:780578/1‑150 (MQ=255)
          gatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGa                                                                                                                                             >  1:1141640/1‑150 (MQ=255)
          gatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGa                                                                                                                                             >  2:593594/1‑150 (MQ=255)
          gatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGa                                                                                                                                             >  2:101561/1‑150 (MQ=255)
          gatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCAAACTTGTTCGCCgct                                                                                                                                             <  1:303283/150‑3 (MQ=255)
            tAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAAt                                                                                                                                           >  1:605410/1‑150 (MQ=255)
             aaaGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCAAACTTGTTCGCCgctatc                                                                                                                                          <  1:246366/150‑6 (MQ=255)
             aaaGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCAAACTTGTTCGCCgctatc                                                                                                                                          <  1:886781/150‑6 (MQ=255)
                 gAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTcgcg                                                                                                                                      <  1:508778/150‑1 (MQ=255)
                  aaCCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTcgcgc                                                                                                                                     <  1:665659/150‑1 (MQ=255)
                    ccTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCgg                                                                                                                                   >  2:1165996/1‑150 (MQ=255)
                    ccTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCgg                                                                                                                                   <  2:247395/150‑1 (MQ=255)
                       ggTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGgtcg                                                                                                                                >  2:250499/1‑150 (MQ=255)
                       ggTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGAGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGgtcg                                                                                                                                >  1:564629/1‑150 (MQ=255)
                        gtcgtcAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGgtcgt                                                                                                                               >  1:582569/1‑150 (MQ=255)
                         tcgtcAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTATCGCGgtcgtc                                                                                                                              <  2:1343183/150‑1 (MQ=255)
                             cAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTa                                                                                                                          >  1:699043/1‑150 (MQ=255)
                                     aCTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCt                                                                                                                  <  2:176225/150‑1 (MQ=255)
                                      cTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCtt                                                                                                                 >  1:1062/1‑150 (MQ=255)
                                            aTCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCcg                                                                                                           >  1:1177551/1‑150 (MQ=255)
                                             tCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCcgc                                                                                                          >  1:1118368/1‑150 (MQ=255)
                                               cacaTTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCgcgc                                                                                                        >  2:722744/1‑150 (MQ=255)
                                               cacaTTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCgcgc                                                                                                        >  2:781712/1‑150 (MQ=255)
                                                acaTTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCa                                                                                                       <  1:154531/150‑1 (MQ=255)
                                                 caTTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCaa                                                                                                      >  1:1192637/1‑150 (MQ=255)
                                                          aGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCgg                                                                                             >  1:390653/1‑150 (MQ=255)
                                                             cAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGgcag                                                                                          >  2:1127996/1‑150 (MQ=255)
                                                             cAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGgcag                                                                                          >  2:1428620/1‑150 (MQ=255)
                                                              aaGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGgcagc                                                                                         <  1:593594/150‑1 (MQ=255)
                                                              aaGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGgcagc                                                                                         <  2:877374/150‑1 (MQ=255)
                                                               aGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGgcagca                                                                                        >  2:346505/1‑150 (MQ=255)
                                                               aGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGgcagca                                                                                        >  2:674162/1‑150 (MQ=255)
                                                               aGTCGCACTGGCTGTTGCTGGCGGAGAGGGTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGGCTGTAGCGAACCTTATTCCGCGAAAGTTGCGCggaataa                                                                                        <  1:1243111/150‑6 (MQ=255)
                                                                gTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCAt                                                                                       >  2:521347/1‑150 (MQ=255)
                                                                        ggCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGacac                                                                               <  2:1184592/150‑1 (MQ=255)
                                                                         gCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACAcc                                                                              <  1:564920/150‑1 (MQ=255)
                                                                         gCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACAcc                                                                              <  2:926488/150‑1 (MQ=255)
                                                                         gCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACAcc                                                                              <  2:660754/150‑1 (MQ=255)
                                                                           tGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCaa                                                                            >  1:689272/1‑150 (MQ=255)
                                                                            gTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAAc                                                                           <  1:323359/150‑1 (MQ=255)
                                                                              tGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGt                                                                         >  1:1408159/1‑150 (MQ=255)
                                                                               gCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTg                                                                        <  2:1071653/150‑1 (MQ=255)
                                                                                 tGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCt                                                                      <  2:780578/150‑1 (MQ=255)
                                                                                  ggCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTg                                                                     <  1:521347/150‑1 (MQ=255)
                                                                                   gCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTTCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTgg                                                                    >  1:1387010/1‑150 (MQ=255)
                                                                                    cGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGc                                                                   <  2:1192637/150‑1 (MQ=255)
                                                                                       agagGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCttt                                                                >  2:235709/1‑150 (MQ=255)
                                                                                        gagGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCtttt                                                               <  2:1142162/150‑1 (MQ=255)
                                                                                          ggTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTgg                                                             >  1:190619/1‑150 (MQ=255)
                                                                                          ggTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTgg                                                             <  1:771129/150‑1 (MQ=255)
                                                                                          ggTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTgg                                                             >  2:963469/1‑150 (MQ=255)
                                                                                          ggTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTgg                                                             <  2:1330829/150‑1 (MQ=255)
                                                                                           gTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGt                                                            >  1:649858/1‑150 (MQ=255)
                                                                                              gATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGATCCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTTGttga                                                         <  1:250499/150‑3 (MQ=255)
                                                                                              gATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCa                                                         <  1:722744/150‑1 (MQ=255)
                                                                                              gATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCa                                                         <  2:601646/150‑1 (MQ=255)
                                                                                                tGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTTGCTTTTGGTTCACg                                                       <  2:1177551/150‑1 (MQ=255)
                                                                                                  gCGTGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAg                                                     <  1:858708/150‑1 (MQ=255)
                                                                                                  gCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGGCGTCTGTAGCGAACCTTAATCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAg                                                     <  2:1187566/150‑1 (MQ=255)
                                                                                                        ttttGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTg                                               <  2:6284/150‑1 (MQ=255)
                                                                                                            gATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCt                                           >  1:85363/1‑150 (MQ=255)
                                                                                                             aTTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTc                                          <  1:674162/150‑1 (MQ=255)
                                                                                                              tttGTGGTACTGGCGTCGGTCTTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCg                                         <  2:4899/150‑1 (MQ=255)
                                                                                                               ttGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGGCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGa                                        <  1:885560/150‑1 (MQ=255)
                                                                                                                 gtgGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTTCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAAc                                      <  2:1387010/150‑1 (MQ=255)
                                                                                                                   ggTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTg                                    >  2:528203/1‑150 (MQ=255)
                                                                                                                    gTACTGGCGTCGGTATTTCGCTAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTgg                                   <  2:149398/150‑1 (MQ=255)
                                                                                                                     tACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGc                                  <  1:937103/150‑1 (MQ=255)
                                                                                                                     tACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGc                                  <  2:231293/150‑1 (MQ=255)
                                                                                                                       cTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCaa                                >  1:424245/1‑150 (MQ=255)
                                                                                                                         ggCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCaaaa                              >  2:713411/1‑150 (MQ=255)
                                                                                                                         ggCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAAACTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCaaaa                              >  2:285304/1‑150 (MQ=255)
                                                                                                                              cGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGAt                         <  1:983115/150‑1 (MQ=255)
                                                                                                                               ggTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGAtt                        >  1:465114/1‑150 (MQ=255)
                                                                                                                               ggTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGAGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGAtt                        <  1:101561/150‑1 (MQ=255)
                                                                                                                                  aTTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATtgtg                     <  1:781712/150‑1 (MQ=255)
                                                                                                                                  aTTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATtgtg                     <  1:711023/150‑1 (MQ=255)
                                                                                                                                   tttCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTgg                    <  2:831896/150‑1 (MQ=255)
                                                                                                                                   tttCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTgg                    <  2:969351/150‑1 (MQ=255)
                                                                                                                                    ttCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGa                   <  2:859642/150‑1 (MQ=255)
                                                                                                                                     tCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGAt                  >  2:102204/1‑150 (MQ=255)
                                                                                                                                       aattccggcaAACTTGTTCGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGc                >  1:1301122/6‑150 (MQ=255)
                                                                                                                                        aTAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCg               >  1:1254867/1‑150 (MQ=255)
                                                                                                                                             ggcgAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCt          >  2:349400/1‑150 (MQ=255)
                                                                                                                                             ggcgAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCt          >  1:1032670/1‑150 (MQ=255)
                                                                                                                                               cgAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTgg        <  2:1254867/150‑1 (MQ=255)
                                                                                                                                               cgAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTgg        >  1:398618/1‑150 (MQ=255)
                                                                                                                                                 aaCAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGgc      >  2:716550/1‑150 (MQ=255)
                                                                                                                                                  aCAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGgcg     <  1:346505/150‑1 (MQ=255)
                                                                                                                                                    aaGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGGTCACGAGTGGGTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGca   <  2:897059/150‑1 (MQ=255)
                                                                                                                                                     aGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGcac  <  1:877628/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGAAGTGATTGATAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCAC  >  NC_000913/4313433‑4313731

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: