Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,077,794 G→C 5.8% D76H (GAC→CAC)  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,7940GC 5.8% 357.0 / 9.3 104D76H (GAC→CAC) cbtACP4‑44 prophage; cytoskeleton‑binding toxin CbtA
Reads supporting (aligned to +/- strand):  ref base G (47/51);  new base C (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 = 5.95e-01

Only 100 of 106 total aligned reads displayed.

TGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGA  >  NC_000913/2077645‑2077943
                                                                                                                                                     |                                                                                                                                                     
tgTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCg                                                                                                                                                       <  1:587694/150‑1 (MQ=255)
 gTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGa                                                                                                                                                      >  1:869335/1‑150 (MQ=255)
  tCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGAc                                                                                                                                                     <  1:415466/150‑1 (MQ=255)
   cccGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGAcc                                                                                                                                                    <  1:493954/150‑1 (MQ=255)
   cccGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGAcc                                                                                                                                                    <  2:1289420/150‑1 (MQ=255)
   cccGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGAcc                                                                                                                                                    >  1:252545/1‑150 (MQ=255)
    ccGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCa                                                                                                                                                   >  2:265060/1‑150 (MQ=255)
    ccGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCa                                                                                                                                                   >  1:1353907/1‑150 (MQ=255)
    ccGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCa                                                                                                                                                   >  2:259665/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)
           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                                                                                                                    <  2:679463/150‑1 (MQ=255)
                                   cTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAg                                                                                                                    <  1:322608/150‑1 (MQ=255)
                                   cTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAg                                                                                                                    >  1:432595/1‑150 (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)
                                                       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                                                                                        <  2:847510/150‑1 (MQ=255)
                                                               gtgtGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCggg                                                                                        >  1:1259219/1‑150 (MQ=255)
                                                                tgtgATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGc                                                                                       >  2:68758/1‑150 (MQ=255)
                                                                tgtgATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGc                                                                                       >  2:216079/1‑150 (MQ=255)
                                                                tgtgATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGAGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGc                                                                                       >  1:980195/1‑150 (MQ=255)
                                                                   gATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGAt                                                                                                                                           <  1:1523087/95‑1 (MQ=255)
                                                                   gATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGAt                                                                                                                                           >  2:1523087/1‑95 (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)
                                                                           agcaTATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGa                                                                            >  1:235222/1‑150 (MQ=255)
                                                                            gcaTATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGAc                                                                           >  1:578182/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                                                        <  1:375637/150‑1 (MQ=255)
                                                                                               tCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGAc                                                        <  2:1445558/150‑1 (MQ=255)
                                                                                                cACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACa                                                       <  1:1294154/150‑1 (MQ=255)
                                                                                                 aCTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACaa                                                      <  2:1062615/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:311693/150‑1 (MQ=255)
                                                                                                               tGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTaa                                        <  1:932257/150‑1 (MQ=255)
                                                                                                                 aaCTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAAt                                      <  2:542769/150‑1 (MQ=255)
                                                                                                                   cTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATaa                                    >  2:5086/1‑150 (MQ=255)
                                                                                                                   cTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATaa                                    >  2:1355858/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:113792/150‑1 (MQ=255)
                                                                                                                                   tACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAg                    <  2:1349319/150‑1 (MQ=255)
                                                                                                                                    aCGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGt                   <  2:606673/150‑1 (MQ=255)
                                                                                                                                     cgcgCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTa                  <  2:1530292/150‑1 (MQ=255)
                                                                                                                                     ccggcTGGTCGGTACGCACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAggg                                                                               <  1:1143921/86‑1 (MQ=38)
                                                                                                                                     ccggcTGGTCGGTACGCACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAggg                                                                               >  2:1143921/4‑89 (MQ=38)
                                                                                                                                       cgcTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATc                >  1:1161104/1‑150 (MQ=255)
                                                                                                                                          tGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCgg             <  1:681363/150‑1 (MQ=255)
                                                                                                                                          tGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCgg             >  2:1468592/1‑150 (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:1317630/150‑1 (MQ=255)
                                                                                                                                             tGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAgg          <  2:432595/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)
                                                                                                                                                     gACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAatga  >  1:1149477/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGA  >  NC_000913/2077645‑2077943

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: