Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,985 C→A 18.8% intergenic (+87/‑383) yffL → / → yffM CPZ‑55 prophage; uncharacterized protein YffL/CPZ‑55 prophage; uncharacterized protein YffM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,9850CA18.8% 140.8 / 3.4 72intergenic (+87/‑383)yffL/yffMCPZ‑55 prophage; uncharacterized protein YffL/CPZ‑55 prophage; uncharacterized protein YffM
Reads supporting (aligned to +/- strand):  ref base C (25/30);  new base A (9/7);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 5.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAG  >  NC_000913/2560840‑2561132
                                                                                                                                                  |                                                                                                                                                   
cGGCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTCaaaa                                                                                                                                                 >  2:1223401/1‑150 (MQ=255)
 ggCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGTGTAaaaat                                                                                                                                                >  2:1002861/1‑149 (MQ=255)
 ggCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAaaaat                                                                                                                                                >  2:248199/1‑149 (MQ=255)
  gCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGt                                                                                                                                               <  1:1115419/150‑1 (MQ=255)
    aGAAATGGAGTGCAAAGAATACCAGTTTTTATTGAAGCCGTACAACATGATATAAAGGAATCAGGCGGGGAATTTTCCGGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGTTCAAAAGTTg                                                                                                                                             <  2:441243/150‑1 (MQ=255)
      aaaGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTggc                                                                                                                                           <  2:354265/150‑3 (MQ=255)
          ggAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTa                                                                                                                                       <  1:459646/150‑1 (MQ=255)
              tGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                                                   <  2:1049820/150‑1 (MQ=255)
              tGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCCTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAATTTGCCGTTAAGGt                                                                                                                                   >  1:1187792/1‑150 (MQ=255)
               gCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTg                                                                                                                                  <  2:93966/150‑1 (MQ=255)
               gCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGGTGGGTAAAAAGTTGACGTTAAggag                                                                                                                                  >  1:421886/1‑148 (MQ=255)
               gCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGATTATGTTGGGTCAAAATTTGAATTTAAggag                                                                                                                                  >  2:812028/1‑148 (MQ=255)
                cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAAGTTGCCGTTAaggagg                                                                                                                                 >  1:896806/1‑150 (MQ=255)
                 aaCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTggt                                                                                                                                <  1:1032324/150‑2 (MQ=255)
                 aaCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTggg                                                                                                                                <  1:427110/150‑1 (MQ=255)
                 aaCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGGCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTggg                                                                                                                                <  2:323698/150‑1 (MQ=255)
                  aCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGt                                                                                                                               <  1:1130766/150‑1 (MQ=255)
                  aCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGt                                                                                                                               >  1:663530/1‑150 (MQ=255)
                        aCCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGTGTAAAATTTTGCCTTTTAGGTGTGTGaaatt                                                                                                                         >  1:1330997/1‑148 (MQ=255)
                        aCCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTGAGTGGTGGGTAAAAAGTTGACTTAAAGGTGGGTGaatat                                                                                                                         >  1:640378/1‑147 (MQ=255)
                         ccAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGt                                                                                                                        <  1:658267/150‑1 (MQ=255)
                          cAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAATTTGCCGTTAAGGTGTGTAAAAAGtt                                                                                                                       >  1:333511/1‑150 (MQ=255)
                            ttaatttttCAAGCAGTACAACATGAAATAAAGGAATCAGGCTTGCATTTTTCCCGCCATTTCTTTTCGTTTATTTAATGAATTCAGGTGGGTAAATTGTTGCCGTTAGTTGTGGGTCAAAAGTTTCCTTTAAGGTGGGTCAAAAGTTGc                                                                                                                     <  2:441808/147‑1 (MQ=19)
                             ttattaTTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGcc                                                                                                                    <  1:28184/150‑1 (MQ=255)
                               attattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGt                                                                                                                  <  1:1397400/150‑1 (MQ=255)
                                   ttCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGGCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAg                                                                                                              <  1:744038/150‑1 (MQ=255)
                                    tCAAGCAGTAGAACATGAAATAACGGAATCAGTCTTGCAATGTTCACGCATTTTGTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGGCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAgg                                                                                                             <  2:1321369/150‑1 (MQ=37)
                                         ccGTACAACATGAAATTACGGAATCAGTATTGCAATTTTTACTCCTTTTCTTTTCCTTAATTCAATGAATTCAGGTGGGTAAAATGTTGCCGTTAGTTTTGGGTAAAAAGTTTCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGt                                                                                                        <  2:617167/150‑1 (MQ=25)
                                              caacaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGGCAAAAGTTGCCGTTAAGGGGGGTCAAAAGTTGCAGTGAAGggggcggaaaa                                                                                                   >  1:825549/1‑143 (MQ=255)
                                               aacaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGGCAAAAg                                                                                                  <  2:1641940/150‑1 (MQ=255)
                                                 caTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGtt                                                                                                >  1:360380/1‑150 (MQ=255)
                                                  aTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTg                                                                                               >  1:404026/1‑150 (MQ=255)
                                                   tGAAATAACTGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGc                                                                                              <  1:1319052/150‑1 (MQ=255)
                                                   tGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTaa                                                                                                                                      <  1:1460049/110‑1 (MQ=255)
                                                   tGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTaa                                                                                                                                      >  2:1460049/1‑110 (MQ=255)
                                                        tAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGt                                                                                         <  1:1134532/150‑1 (MQ=255)
                                                         aaCGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGGCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTa                                                                                        >  1:830356/1‑150 (MQ=255)
                                                           cGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTTAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAg                                                                                      <  1:500272/150‑1 (MQ=255)
                                                            ggAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGACAAAAGTTGCCGGTAAGGTGGGTCAACAGTTGCCGAGTAAGa                                                                                     >  1:726015/1‑150 (MQ=255)
                                                               aTCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTACTGGTGGGTCAAAAGTTGCCGTTAAGGTGAGTCAAAAGTTGCAGAAAAGGTGGGGCAAAAGTTGCCGagaaagaagt                                                                                  >  2:1572519/1‑148 (MQ=37)
                                                                  aGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTTAGGGGGGGCACATGTTGCCGAGTAAGAAGTcta                                                                               >  1:1310731/1‑150 (MQ=37)
                                                                   ggCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGTTGGGTCAAAATTTGCCGTTAAGGTGGGTCAAAATTTGAAGATAAGGTGTGTCAAAATTTGCAGAGAAAGAATTctaa                                                                              >  2:207539/1‑149 (MQ=37)
                                                                              ttCACGCATTTTCTTTCCCATTAATTAATGAATTCTGTTGTGTATAAATTAGGCTTTATTGTGGGGTAAATTCTTGCCTTTAAGTTGGGTCAAATGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAAtt                                                                   <  2:333511/150‑1 (MQ=255)
                                                                                      ttttGTTTTCCTTAATTTATTGAATTCAGTTTGGTTAAATTTTGGCTTTTTTGTTGGTTAAATTGTTGCCTTTAAGGTGGGTCAAAAGTTGCATTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGa                                                           <  2:421886/150‑1 (MQ=19)
                                                                                      ttttCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCCTTACTGGTGGGTAAAAACTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGa                                                           <  1:1457947/150‑1 (MQ=38)
                                                                                           tttGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGGTGTGTAAAAATTTGCCTTTAAGGAGGCTAAAAATTTGAAGTAAAGGTGGCGGAAAATTTGAAGAGAAACAATTATACAATTAAAAATAATACGCGACGTGa                                                      >  2:1561344/1‑150 (MQ=255)
                                                                                                cTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAACGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGACTTCTACTATTAAGAATTCTAGGCGACGGGAGtttt                                                 >  1:1228432/1‑150 (MQ=255)
                                                                                                  ttAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGAGGGTCAAAATTTGACGTTAAGGTGGCGGAAATTTTGACGAGTAACAATTCTACAATTAAAAATTCTACGCGACGAGAAttttt                                                >  2:954465/1‑150 (MQ=25)
                                                                                                    aaatCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTACTGGTGGGTCAAAACTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTa                                              <  1:419395/147‑1 (MQ=255)
                                                                                                        ccATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTTAGGTGGGGGAAATTTTGACGAGGAAGGATTCTACAATTAAGAATTCTAGGCGACGGGAGTTTTTTAGcat                                          >  1:674669/1‑148 (MQ=25)
                                                                                                           tGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGACGTTAAGGTGGGTCAAAAGTTTCCGAGTAAGAATTCTACTATTAAGAATTATAGGCGACGGGAGTTTTTTAGACtaac                                       >  2:152338/1‑147 (MQ=38)
                                                                                                           tGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGGTGGGTAAAAAGTTCCAGTTTAGCAGTCTCAAATTTTTAAGTTAACCTGTCTGAAAAGTTGACCAGTAATAAGACGACTATAAAGAAATCAAGGAGACGGGATTATTTTATCCtgat                                       >  2:829240/1‑147 (MQ=255)
                                                                                                            gAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGAGGGTCAAAATTTGCAGATTACGTGGGTCAAAATTTGACGAGTAAGCATTATAATATTAAAAATTCTAGGCGACGTGAGATTTTTAGACtaacg                                      >  2:970662/1‑146 (MQ=21)
                                                                                                             aaTTCAGGTGGGTAAAAATTAGGCCTTATTGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcgg                                     <  1:812028/150‑1 (MQ=37)
                                                                                                              aTTCAGGTGGGTCAAAAGTTGCCGTTACTGGTGGGTCAAAATTTGCCGTTAAGGAGTGTCAAAAGTTGACGTAAAGGTGGGTCAAAATTTGACGAGTAAGAATTATAAAATTAAAAATTCTAGGCCACGGGAGTTTTTTAGACTAAcggt                                    >  2:943707/1‑150 (MQ=25)
                                                                                                                 cAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTACcggtcgg                                 >  1:373857/1‑150 (MQ=255)
                                                                                                                 aaGGTGGCTCAAAAGATGTAGTTATGGGTGGTTCATTTCTTGCCTTTTATTTGGGTTAAATGTTTCATTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggtcgg                                 <  2:1169302/149‑1 (MQ=255)
                                                                                                                       tgtaaaaaTTTGCCCTTAGTGGTGGGTCAACACTTGCCGTCAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATag                           <  2:536599/146‑1 (MQ=38)
                                                                                                                       ggTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGGCGTTAAGGTGGGTCAAAATTTGCAGAGTAAGAATTCGACAATTTAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATag                           >  1:536599/1‑150 (MQ=38)
                                                                                                                        gTCAAAAGTTGCCGTTAGTGGTGGGACAAAAGTTGCCGTTAAGGTGTTTCAAAATTTGCCGTTAAGGTGGGTCAAAATTTTAAGAGTAAGAATTCTAAAATTAAGAATACAAGGCGAAGGAAGATTTTTATAATACCGCTCGGTAATaga                          >  2:495583/1‑150 (MQ=19)
                                                                                                                          aaaaaGTTGCCTTTAATGGTGGGTAATAAGTTGCCTTTAAGGTGGGTAAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATagaga                        <  2:1652434/149‑1 (MQ=38)
                                                                                                                           aattgTTGGCGTTTGTGGTGGGTCACAACCTGCCGTTATGGTGGGTCAAAAGTTTCCGTTTTGGTGGGTCAATAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTTGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATagagag                       <  2:336227/146‑1 (MQ=255)
                                                                                                                           aaaaGTAGGCGTTAGTGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATagagag                       <  1:1433419/150‑1 (MQ=255)
                                                                                                                              tctctccGTTACTGGGGGTTCAATTCTTGCCTTTAAGGTGGGTAAAAAGTTGCATTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCg                    <  2:572697/145‑1 (MQ=255)
                                                                                                                                tttccGTTAGTGGTCGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGcc                  <  1:1550707/147‑1 (MQ=255)
                                                                                                                                  ggctttATTGGTGGGTAAATAGTTGCCTTTAAGTTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTa                <  2:663530/146‑1 (MQ=38)
                                                                                                                                   ccGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTag               >  1:1040508/1‑150 (MQ=255)
                                                                                                                                          tgggggGTATAATGTTGCCTTTTAGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTaaa        <  2:825549/150‑1 (MQ=38)
                                                                                                                                           ggtggGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTaaaa       >  2:1306147/1‑149 (MQ=255)
                                                                                                                                           ggtggCTCACAACTTGCCGTTAAGGTGCGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAAc       <  1:536479/150‑1 (MQ=255)
                                                                                                                                           ggggggTCAATAGTAGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTTGGTCAAAAGGTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAAc       <  2:234327/147‑1 (MQ=39)
                                                                                                                                                gTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACagaag  >  1:173971/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGGCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAG  >  NC_000913/2560840‑2561132

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: