Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,455,145 C→A 6.3% A168S (GCT→TCT) ‡ yfcV ← putative fimbrial protein YfcV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,455,1450CA 6.3% 222.5 / 3.8 68A168S (GCT→TCT) ‡yfcVputative fimbrial protein YfcV
Reads supporting (aligned to +/- strand):  ref base C (36/27);  new base A (2/2);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

GGCATACACGTTTTTATCCGTTTTTATTCTTGTTCGGTCAAAATCCATTTAAATACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCAGAGAGTCACCGATCGCCAGGCTAACATTGTTAA  >  NC_000913/2454997‑2455294
                                                                                                                                                    |                                                                                                                                                     
ggCATACACGTTTTTATCCGTTTTTATTCTTGTTCGGTCAAAATCCATTTAAATACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCt                                                                                                                                                      >  1:787833/1‑150 (MQ=255)
ggCATACACGTTTTTATCCGTTTTTATTATTGTTCGGTCAAAATCCATTTAAATACCCCATAGCAGTGTCGCTATGGGGTTATGCCTAACATGTAATTAATAGGGAAGGTGGTATTTGCTTAGAAATTACCCAGATCAGGCGCATCAGca                                                                                                                                                      >  2:1194001/1‑149 (MQ=255)
    tACACGTTTTTATCCGTTTTTATTCTTGTTCGGTCAAAATCCATTTAAATACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTgcgc                                                                                                                                                  >  1:1222375/1‑150 (MQ=255)
      cacGTTTTTATCCGTTTTTATTCTTGTTCGGTCAAAATCCATTTAAATACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCg                                                                                                                                                >  2:787026/1‑150 (MQ=255)
       acGTTTTTATCCGTTTTTATTCTTGTTCGGTCAAAATCCATTTAAATACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGa                                                                                                                                               <  1:256312/150‑1 (MQ=255)
       acGTTTTTATCCGTTTTTATTCTTGTTCGGTCAAAATCCATTTAAATACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCAGCTGATGCGccag                                                                                                                                               <  2:357310/150‑3 (MQ=255)
        cGTTTTTATCCGTTTTTATTCTTGTTCGGTCAAAATCCATTTAAATACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCAGCTGTTGCGccaga                                                                                                                                              <  2:1117931/150‑4 (MQ=255)
         gTTTTTATCCGTTTTTATTCTTGTTCGGTCAAAATCCATTTAAATACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGAcc                                                                                                                                             <  2:1271230/150‑1 (MQ=255)
              tATCCGTTTTTATTCTTGTTCGGTCAAAATCCATTTAAATACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGAccagcca                                                                                                                                        <  2:908748/150‑1 (MQ=255)
                tCCGTTTTTATTCTTGTTCGGTCAAAATCCATTTAAATACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGAAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAgg                                                                                                                                      >  1:342484/1‑150 (MQ=255)
                 ccGTTTTTATTCTTGTTCGGTCAAAATCCATTTAAATACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGc                                                                                                                                     <  1:228796/150‑1 (MQ=255)
                    tttttATTCTTGTTCGGTCAAAATCCATTTAAATACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCttt                                                                                                                                  <  2:477773/150‑1 (MQ=255)
                        tATTCTTGTTCGGTCAAAATCCATTTAAATACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAg                                                                                                                              <  1:1009531/150‑1 (MQ=255)
                                     tCAAAATCCATTTAAATACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCt                                                                                                                 <  1:1342332/150‑1 (MQ=255)
                                         aaTCCATTTAAATACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGc                                                                                                             >  1:84403/1‑150 (MQ=255)
                                           tCCATTTAAATACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCtt                                                                                                           <  2:1352051/150‑1 (MQ=255)
                                              aTTTAAATACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGc                                                                                                        >  2:721171/1‑150 (MQ=255)
                                                 tAAATACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCtt                                                                                                     >  2:463349/1‑150 (MQ=255)
                                                  aaaTACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCttt                                                                                                    >  1:791795/1‑150 (MQ=255)
                                                         ccATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTg                                                                                             >  2:529217/1‑150 (MQ=255)
                                                                aGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTAt                                                                                      >  2:1161359/1‑150 (MQ=255)
                                                                  tgtCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCt                                                                                    <  2:787833/150‑1 (MQ=255)
                                                                   gtCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCtt                                                                                   >  2:803870/1‑150 (MQ=255)
                                                                     cGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGGGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTc                                                                                 <  2:904538/150‑1 (MQ=255)
                                                                      gCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCa                                                                                >  2:994017/1‑150 (MQ=255)
                                                                       cTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCAc                                                                               <  1:73023/150‑1 (MQ=255)
                                                                       cTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCAc                                                                               >  2:383835/1‑150 (MQ=255)
                                                                        tATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACg                                                                              <  1:463349/150‑1 (MQ=255)
                                                                          tGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGAt                                                                            >  1:276670/1‑150 (MQ=255)
                                                                               gTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCt                                                                       <  2:922703/150‑1 (MQ=255)
                                                                                tATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTg                                                                      >  2:1339583/1‑150 (MQ=255)
                                                                                   gCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTc                                                                   <  1:1000280/150‑1 (MQ=255)
                                                                                   gCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTc                                                                   >  1:143290/1‑150 (MQ=255)
                                                                                    ccTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCa                                                                  <  2:1222375/150‑1 (MQ=255)
                                                                                        aCAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATAc                                                              >  1:743536/1‑150 (MQ=255)
                                                                                             taagtaaTCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATg                                                         <  2:276670/150‑1 (MQ=255)
                                                                                             taagtaaTCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATg                                                         >  2:149386/1‑150 (MQ=255)
                                                                                                    tCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCttt                                                  >  1:1413466/1‑150 (MQ=255)
                                                                                                    tCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCttt                                                  <  1:1161359/150‑1 (MQ=255)
                                                                                                     cTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCtttt                                                 >  2:392155/1‑150 (MQ=255)
                                                                                                     cTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCtttt                                                 >  2:257830/1‑150 (MQ=255)
                                                                                                      tGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTg                                                >  1:170648/1‑150 (MQ=255)
                                                                                                       ggAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGt                                               >  1:140956/1‑150 (MQ=255)
                                                                                                       ggAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGt                                               >  1:146026/1‑150 (MQ=255)
                                                                                                        gAAGTTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTa                                              <  1:1194001/150‑1 (MQ=255)
                                                                                                          aGGTGGGATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCACGCTTTAAAGTTCAGAGTCTGCTTCGGTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAag                                            <  2:465869/150‑1 (MQ=255)
                                                                                                              ggtATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGt                                        <  2:1413466/150‑1 (MQ=255)
                                                                                                                tATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTAc                                      <  1:95873/150‑1 (MQ=255)
                                                                                                                 aTTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCAATGCCGCTTTTGTAAGAGGTAcc                                     >  1:674684/1‑150 (MQ=255)
                                                                                                                    tGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCag                                  >  2:806532/1‑150 (MQ=255)
                                                                                                                     gCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTATGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCaga                                 >  2:434756/1‑150 (MQ=255)
                                                                                                                      cTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCagag                                >  2:1325100/1‑150 (MQ=255)
                                                                                                                      cTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCagag                                >  1:871217/1‑150 (MQ=255)
                                                                                                                       ttCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCagaga                               >  2:579538/1‑150 (MQ=255)
                                                                                                                           aaaTTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCAGAGAGTCa                           >  2:1060161/1‑150 (MQ=255)
                                                                                                                               tACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCAGAGAGTCACCGa                       >  2:167498/1‑150 (MQ=255)
                                                                                                                                 cccAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCTCTTTTGTAAGAGGTACCCAGAGAGTCACCGATc                     >  1:777094/1‑150 (MQ=255)
                                                                                                                                 cccAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCAGAGAGTCACCGATc                     >  1:603028/1‑150 (MQ=255)
                                                                                                                                     gATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCAGAGAGTCACCGATCGCCa                 <  1:392155/150‑1 (MQ=255)
                                                                                                                                       gtcggCGCAGCTGATGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCagagag                              <  2:457585/132‑1 (MQ=255)
                                                                                                                                       gtcggCGCAGCTGATGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCagagag                              >  1:457585/4‑135 (MQ=255)
                                                                                                                                       gtcggCGCAGCTGATGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCAGAGAGTCACCGATCGCCAgg               >  1:825539/4‑150 (MQ=255)
                                                                                                                                          ggCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCAGAGAGTCACCGATCGCCAGGCTa            >  1:48336/1‑150 (MQ=255)
                                                                                                                                            cgcATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCAGAGAGTCACCGATCGCCAGGCTAAc          <  2:791795/150‑1 (MQ=255)
                                                                                                                                             gcATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCAGAGAGTCACCGATCGCCAGGCTAACa         <  1:787026/150‑1 (MQ=255)
                                                                                                                                             gcATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCAGAGAGTCACCGATCGCCAGGCTAACa         >  1:761192/1‑150 (MQ=255)
                                                                                                                                              cATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCAGAGAGTCACCGATCGCCAGGCTAACAt        <  1:1339583/150‑1 (MQ=255)
                                                                                                                                                 cAGCTGCGCCGACCAGCCATGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGATTTTGTAAGAGGTACCGAGAGAGTCACCGATCGCCAGGCTAACATTGt     <  2:1292804/150‑1 (MQ=255)
                                                                                                                                                   gCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCAGAGAGTCACCGATCGCCAGGCTAACATTGTTa   <  2:48483/150‑1 (MQ=255)
                                                                                                                                                    cTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCAGAGAGTCACCGATCGCCAGGCTAACATTGTTaa  <  1:195541/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCATACACGTTTTTATCCGTTTTTATTCTTGTTCGGTCAAAATCCATTTAAATACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCAGAGAGTCACCGATCGCCAGGCTAACATTGTTAA  >  NC_000913/2454997‑2455294

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: