Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,077,793 C→G 5.5% T75T (ACC→ACG cbtA → CP4‑44 prophage; cytoskeleton‑binding toxin CbtA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,077,7930CG 5.5% 351.0 / 7.1 104T75T (ACC→ACGcbtACP4‑44 prophage; cytoskeleton‑binding toxin CbtA
Reads supporting (aligned to +/- strand):  ref base C (47/51);  new base G (3/3);  total (50/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

Only 100 of 106 total aligned reads displayed.

CTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATG  >  NC_000913/2077644‑2077942
                                                                                                                                                     |                                                                                                                                                     
ctgTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTAcc                                                                                                                                                       <  2:914160/150‑1 (MQ=255)
 tgTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCg                                                                                                                                                      <  1:587694/150‑1 (MQ=255)
  gTCCCGACTGCTGGACCAGCGCTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGa                                                                                                                                                     >  2:1453841/1‑150 (MQ=255)
  gTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGa                                                                                                                                                     >  1:869335/1‑150 (MQ=255)
   tCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGAc                                                                                                                                                    <  1:415466/150‑1 (MQ=255)
    cccGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGAcc                                                                                                                                                   >  1:252545/1‑150 (MQ=255)
    cccGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGAcc                                                                                                                                                   <  2:1289420/150‑1 (MQ=255)
    cccGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGAcc                                                                                                                                                   <  1:493954/150‑1 (MQ=255)
     ccGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCa                                                                                                                                                  >  2:265060/1‑150 (MQ=255)
     ccGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCa                                                                                                                                                  >  2:259665/1‑150 (MQ=255)
     ccGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCa                                                                                                                                                  >  1:1353907/1‑150 (MQ=255)
         ctgctgGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTCGGTACGCACCAGcgc                                                                                                                                              <  1:308938/150‑3 (MQ=255)
          tgctgGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCgg                                                                                                                                             >  2:598604/1‑150 (MQ=255)
          tgctgGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCgg                                                                                                                                             >  2:311693/1‑150 (MQ=255)
          tgctgGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTCGGTACGCACCAGcgcg                                                                                                                                             <  2:360231/150‑4 (MQ=255)
           gctgGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCggg                                                                                                                                            >  1:539438/1‑150 (MQ=255)
            ctgGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGa                                                                                                                                           >  1:113792/1‑150 (MQ=255)
             tgGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGAt                                                                                                                                          <  1:371349/150‑1 (MQ=255)
               gACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTc                                                                                                                                        >  2:914354/1‑150 (MQ=255)
                       cTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTg                                                                                                                                <  1:667155/150‑1 (MQ=255)
                            gCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTAccc                                                                                                                           >  1:487971/1‑150 (MQ=255)
                             ccTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCg                                                                                                                          >  2:754312/1‑150 (MQ=255)
                               tCACACTGAATGACACAGCATTTGCCGATGAACGTGTGATTGAGCATAATATTGAGGCAGGCATTTCACTGGGTGATGCGGTGAACTTTCTCGTGGAAAATTACGCGGTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACACGct                                                                                                                        <  1:80377/150‑1 (MQ=255)
                                cacaCTGAATGACACACCTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGctct                                                                                                                      >  2:681363/1‑150 (MQ=255)
                                  cacTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCtctc                                                                                                                     <  2:1130593/150‑1 (MQ=255)
                                    cTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAg                                                                                                                   <  1:322608/150‑1 (MQ=255)
                                    cTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAg                                                                                                                   >  1:432595/1‑150 (MQ=255)
                                    cTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAg                                                                                                                   <  2:679463/150‑1 (MQ=255)
                                       aaTGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTta                                                                                                                <  2:767570/150‑1 (MQ=255)
                                          gACACACCTTTTGCCGGTGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTtaata                                                                                                             >  1:162864/1‑150 (MQ=255)
                                          gACACACCTTTTGCCGGTGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTtaata                                                                                                             >  1:157395/1‑150 (MQ=255)
                                          gACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTtaata                                                                                                             <  1:914354/150‑1 (MQ=255)
                                            cacacCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATaaa                                                                                                           >  2:932257/1‑150 (MQ=255)
                                               acCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAg                                                                                                        <  1:1163707/150‑1 (MQ=255)
                                               acCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCACGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGACTGTACCCGCTCTCAGTTAATAAACAg                                                                                                        >  1:619957/1‑150 (MQ=255)
                                                 cTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCa                                                                                                      >  2:1122966/1‑150 (MQ=255)
                                                     tGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGa                                                                                                  >  1:876379/1‑150 (MQ=255)
                                                      gCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGat                                                                                                 <  2:523876/150‑1 (MQ=255)
                                                        cGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGatat                                                                                               >  2:1294154/1‑150 (MQ=255)
                                                        cGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGatat                                                                                               <  2:1433551/150‑1 (MQ=255)
                                                           tGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATAtcct                                                                                            <  2:191691/150‑1 (MQ=255)
                                                              aCGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCg                                                                                         <  1:265060/150‑1 (MQ=255)
                                                              aCGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCg                                                                                         <  1:259665/150‑1 (MQ=255)
                                                              aCGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCg                                                                                         <  1:1310916/150‑1 (MQ=255)
                                                                gtgtGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCTGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCggg                                                                                       <  1:1088877/150‑1 (MQ=255)
                                                                gtgtGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCggg                                                                                       >  1:1259219/1‑150 (MQ=255)
                                                                gtgtGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCggg                                                                                       <  2:847510/150‑1 (MQ=255)
                                                                 tgtgATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGc                                                                                      >  2:68758/1‑150 (MQ=255)
                                                                 tgtgATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGc                                                                                      >  2:216079/1‑150 (MQ=255)
                                                                 tgtgATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGAGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGc                                                                                      >  1:980195/1‑150 (MQ=255)
                                                                    gATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGAt                                                                                                                                          >  2:1523087/1‑95 (MQ=255)
                                                                    gATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGAt                                                                                                                                          <  1:1523087/95‑1 (MQ=255)
                                                                    gATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCgggtcat                                                                                   >  2:756014/1‑146 (MQ=255)
                                                                       tGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAg                                                                                <  1:259701/150‑1 (MQ=255)
                                                                        gAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAgg                                                                               >  1:535063/1‑150 (MQ=255)
                                                                        gAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAgg                                                                               <  2:1353907/150‑1 (MQ=255)
                                                                         agcagcATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAggg                                                                              <  1:1303443/150‑1 (MQ=255)
                                                                           cagcaTATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCg                                                                            <  2:487971/150‑1 (MQ=255)
                                                                                tatTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCgg                                                                       >  2:101468/1‑150 (MQ=255)
                                                                                    gAGGCAGGCATTTCACTGTGTGATGCGGTGTACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCtg                                                                   >  2:191885/1‑150 (MQ=255)
                                                                                      ggcaggcaTTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCtgat                                                                 >  1:497142/1‑150 (MQ=255)
                                                                                        caggcaTTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCtgatga                                                               >  1:26756/1‑150 (MQ=255)
                                                                                         aggcaTTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGAc                                                              >  2:1091910/1‑150 (MQ=255)
                                                                                          ggcaTTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGAcc                                                             >  1:1133157/1‑150 (MQ=255)
                                                                                                tCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACGCGCTCTCAGTTAATAAACAGCATCGATATCCTGCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGAc                                                       <  2:919842/150‑1 (MQ=255)
                                                                                                tCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGAc                                                       <  2:1445558/150‑1 (MQ=255)
                                                                                                tCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGAc                                                       <  1:375637/150‑1 (MQ=255)
                                                                                                 cACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACa                                                      <  1:1294154/150‑1 (MQ=255)
                                                                                                      tgtgATGCGGTGTACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTAc                                                 <  1:191885/150‑1 (MQ=255)
                                                                                                           tGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAAc                                            >  1:747564/1‑150 (MQ=255)
                                                                                                             cGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACgg                                          <  2:252545/150‑1 (MQ=255)
                                                                                                             cGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACgg                                          <  1:860863/150‑1 (MQ=255)
                                                                                                              ggTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGt                                         >  2:267800/1‑150 (MQ=255)
                                                                                                                tGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTaa                                       <  1:932257/150‑1 (MQ=255)
                                                                                                                tGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTaa                                       <  1:311693/150‑1 (MQ=255)
                                                                                                                    cTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATaa                                   >  2:1355858/1‑150 (MQ=255)
                                                                                                                    cTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATaa                                   >  2:5086/1‑150 (MQ=255)
                                                                                                                         tattGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTa                              >  2:729191/4‑150 (MQ=255)
                                                                                                                              gAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTg                         <  1:1404089/150‑1 (MQ=255)
                                                                                                                               aaaaaTACGCGCTTGTTCGTACCGACGAGCCGGTATTCAGCGCCTGTACCCTCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCTGATGGCGACCTGCCTGATGAGCCGCGACAATTAGAGAACGGTAAATAACATTACCGTgg                        <  2:739176/150‑1 (MQ=255)
                                                                                                                                  aaTACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTa                     <  1:729191/150‑1 (MQ=255)
                                                                                                                                   atcccggcTGGTCGGTACGCACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGCCCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTaa                    >  1:941756/7‑150 (MQ=255)
                                                                                                                                   atcccggcTGGTCGGTACGCACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTaa                    >  2:1321954/7‑150 (MQ=255)
                                                                                                                                    tACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAg                   <  2:1349319/150‑1 (MQ=255)
                                                                                                                                    tACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAg                   <  2:113792/150‑1 (MQ=255)
                                                                                                                                     aCGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGt                  <  2:606673/150‑1 (MQ=255)
                                                                                                                                      cgcgCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTa                 <  2:1530292/150‑1 (MQ=255)
                                                                                                                                      ccggcTGGTCGGTACGCACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAggg                                                                              >  2:1143921/4‑89 (MQ=38)
                                                                                                                                      ccggcTGGTCGGTACGCACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAggg                                                                              <  1:1143921/86‑1 (MQ=38)
                                                                                                                                        cgcTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATc               >  1:1161104/1‑150 (MQ=255)
                                                                                                                                           tGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCgg            <  1:681363/150‑1 (MQ=255)
                                                                                                                                           tGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCgg            <  2:876379/150‑1 (MQ=255)
                                                                                                                                              tGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAgg         <  2:157395/150‑1 (MQ=255)
                                                                                                                                              tGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAgg         <  2:162864/150‑1 (MQ=255)
                                                                                                                                              tGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAgg         <  2:432595/150‑1 (MQ=255)
                                                                                                                                              tGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAgg         <  2:1317630/150‑1 (MQ=255)
                                                                                                                                                cGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCg       <  1:1438999/150‑1 (MQ=255)
                                                                                                                                                  tACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGaa     >  2:140983/1‑150 (MQ=255)
                                                                                                                                                   aCCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGaaa    >  1:617669/1‑150 (MQ=255)
                                                                                                                                                     cGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAatg  <  2:1161104/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATG  >  NC_000913/2077644‑2077942

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: