Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,455,143 A→T 7.4% A168A (GCT→GCA) ‡ yfcV ← putative fimbrial protein YfcV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,455,1430AT 7.4% 196.8 / 4.4 68A168A (GCT→GCA) ‡yfcVputative fimbrial protein YfcV
Reads supporting (aligned to +/- strand):  ref base A (37/26);  new base T (2/3);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

TCAGGCATACACGTTTTTATCCGTTTTTATTCTTGTTCGGTCAAAATCCATTTAAATACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCAGAGAGTCACCGATCGCCAGGCTAACATTGT  >  NC_000913/2454994‑2455291
                                                                                                                                                     |                                                                                                                                                    
tCAGGCATACACGTTTTTATCCGTTTTTATTCTTGTTCGGTCAAAATCCATTTAAATACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCa                                                                                                                                                      <  1:297854/150‑1 (MQ=255)
   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:1161359/150‑1 (MQ=255)
                                                                                                       tCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCttt                                               >  1:1413466/1‑150 (MQ=255)
                                                                                                        cTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCtttt                                              >  2:257830/1‑150 (MQ=255)
                                                                                                        cTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCtttt                                              >  2:392155/1‑150 (MQ=255)
                                                                                                         tGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTg                                             >  1:170648/1‑150 (MQ=255)
                                                                                                          ggAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGt                                            >  1:146026/1‑150 (MQ=255)
                                                                                                          ggAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGt                                            >  1:140956/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                             >  1:871217/1‑150 (MQ=255)
                                                                                                                         cTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCagag                             >  2:1325100/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                           >  1:457585/4‑135 (MQ=255)
                                                                                                                                          gtcggCGCAGCTGATGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCagagag                           <  2:457585/132‑1 (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:761192/1‑150 (MQ=255)
                                                                                                                                                gcATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCAGAGAGTCACCGATCGCCAGGCTAACa      <  1:787026/150‑1 (MQ=255)
                                                                                                                                                 cATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCAGAGAGTCACCGATCGCCAGGCTAACAt     <  1:1339583/150‑1 (MQ=255)
                                                                                                                                                    cAGCTGCGCCGACCAGCCATGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGATTTTGTAAGAGGTACCGAGAGAGTCACCGATCGCCAGGCTAACATTGt  <  2:1292804/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCAGGCATACACGTTTTTATCCGTTTTTATTCTTGTTCGGTCAAAATCCATTTAAATACCCCATAGCAGTGTCGCTATGGGGGTATGCCTTACAGGTAAGTAATCTGGAAGGTGGTATTTGCTTCGAAATTACCCAGATCTGGCGCATCAGCTGCGCCGACCAGCCAGGCTTTAAAGTTCAGAGTCTGCTTCGCTTTGCCTTTATTAGTGGCAGTATCTTTCACGATTTTCTGGTCAATACCCATGCCGCTTTTGTAAGAGGTACCCAGAGAGTCACCGATCGCCAGGCTAACATTGT  >  NC_000913/2454994‑2455291

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: