Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,563,261 T→G 6.3% intergenic (+144/‑331) yffP → / → yffQ CPZ‑55 prophage; uncharacterized protein YffP/CPZ‑55 prophage; uncharacterized protein YffQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,563,2610TG 6.3% 213.9 / 2.0 71intergenic (+144/‑331)yffP/yffQCPZ‑55 prophage; uncharacterized protein YffP/CPZ‑55 prophage; uncharacterized protein YffQ
Reads supporting (aligned to +/- strand):  ref base T (34/32);  new base G (2/3);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACA  >  NC_000913/2563114‑2563406
                                                                                                                                                   |                                                                                                                                                 
cTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTa                                                                                                                                                 >  1:834622/1‑150 (MQ=255)
 tGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAg                                                                                                                                                <  1:1087604/150‑1 (MQ=255)
 tGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAg                                                                                                                                                <  1:1178261/150‑1 (MQ=255)
 tGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAg                                                                                                                                                <  1:921011/150‑1 (MQ=255)
       ttCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAAcca                                                                                                                                          >  2:1502373/1‑149 (MQ=255)
         cGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTAc                                                                                                                                        >  1:1471027/1‑150 (MQ=255)
          ggACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTAcc                                                                                                                                       <  1:918891/150‑1 (MQ=255)
                   cgcgcgTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAACCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACCTACACTCGCCAAAAATAAACGCCCGGGTAGGTTGCTAACTGCCTGGGCCttt                                                                                                                              <  2:563627/150‑1 (MQ=255)
                     cgcgTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTtct                                                                                                                            >  2:1411380/1‑150 (MQ=255)
                       cgTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTca                                                                                                                          <  2:1196905/150‑1 (MQ=255)
                        gTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTcac                                                                                                                         <  2:97639/150‑1 (MQ=255)
                         tAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTcaca                                                                                                                        >  1:556230/1‑150 (MQ=255)
                          aGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACaa                                                                                                                       <  2:467176/150‑1 (MQ=255)
                                cGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCg                                                                                                                 >  1:355428/1‑150 (MQ=255)
                                     gCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGAc                                                                                                            <  2:195234/150‑1 (MQ=255)
                                       ggCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTg                                                                                                              <  1:1274496/146‑1 (MQ=255)
                                       ggCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTg                                                                                                              >  2:1274496/1‑146 (MQ=255)
                                       ggCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACAt                                                                                                          <  1:204255/150‑1 (MQ=255)
                                        gCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACAtt                                                                                                         <  2:765429/150‑1 (MQ=255)
                                             aCTGGCGGGCAATTCCCCCCTCCAACCCCAAAAATCCCACACAACAACCCTTCAAATTTTCCACATTTCAACTACACGCCCCAAAAAAAAACGCCCACGCACGTACCAACCCACCCCGGCCTTTCTCACAATCCCCGCTGACATTATATc                                                                                                    <  2:441344/150‑1 (MQ=255)
                                                ggCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCAcc                                                                                                 >  2:363360/1‑150 (MQ=255)
                                                ggCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCAcc                                                                                                 >  2:851534/1‑150 (MQ=255)
                                                ggCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCAcc                                                                                                 <  2:556230/150‑1 (MQ=255)
                                                 gCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCAccc                                                                                                <  1:1526791/150‑1 (MQ=255)
                                                      cATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCtt                                                                                           >  1:1479817/1‑150 (MQ=255)
                                                        tttCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTc                                                                                         <  1:683393/150‑1 (MQ=255)
                                                           ccccGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGt                                                                                      >  2:887839/1‑150 (MQ=255)
                                                           ccccGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGt                                                                                      >  1:1159575/1‑150 (MQ=255)
                                                               gCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTc                                                                                  <  2:85606/150‑1 (MQ=255)
                                                                cTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCa                                                                                 <  2:616644/150‑1 (MQ=255)
                                                                  ccTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAgg                                                                               <  1:1502373/150‑1 (MQ=255)
                                                                   cTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTGAGGt                                                                              <  2:355428/150‑1 (MQ=255)
                                                                    tACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTg                                                                             <  2:834622/150‑1 (MQ=255)
                                                                      cgcgAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTggg                                                                           <  2:1479817/150‑1 (MQ=255)
                                                                      cgcgAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTggg                                                                           <  1:327100/150‑1 (MQ=255)
                                                                      cgcgAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTggg                                                                           >  2:507041/1‑150 (MQ=255)
                                                                         gAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGgaga                                                                        >  1:676916/1‑150 (MQ=255)
                                                                         gAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGgaga                                                                        >  2:1226105/1‑150 (MQ=255)
                                                                          aaaaaTCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGAc                                                                       <  1:375399/150‑1 (MQ=255)
                                                                             aaTCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGa                                                                    <  2:905794/150‑1 (MQ=255)
                                                                                  cTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGgcc                                                               >  2:1248039/1‑150 (MQ=255)
                                                                                    cacaACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGgccgc                                                             >  2:1608375/1‑150 (MQ=255)
                                                                                      caACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGccgccg                                                           >  2:701602/1‑150 (MQ=255)
                                                                                      caACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGccgccg                                                           >  2:156038/1‑150 (MQ=255)
                                                                                            cccTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTg                                                     >  1:121100/1‑150 (MQ=255)
                                                                                               tGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGcac                                                  >  1:847363/1‑150 (MQ=255)
                                                                                                   aaTTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAAtt                                              <  1:887839/150‑1 (MQ=255)
                                                                                                   aaTTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAAtt                                              <  2:1162369/150‑1 (MQ=255)
                                                                                                      tttCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCCCAATCCCCGCTGACATTATATAACCCATCTTTGCGTTTTCAGGTGGGAGACTGATGGCCGCCGCCGTTGCACAAATGAt                                           >  2:950984/1‑150 (MQ=255)
                                                                                                        tCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTa                                         >  1:547599/1‑150 (MQ=255)
                                                                                                         ccAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTaa                                        >  1:1108329/1‑150 (MQ=255)
                                                                                                           aGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAAtt                                      <  1:1226105/150‑1 (MQ=255)
                                                                                                             aTTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTc                                    <  1:424938/150‑1 (MQ=255)
                                                                                                               ttCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGCCTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTg                                  >  1:1064547/1‑150 (MQ=255)
                                                                                                                     tACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACAtt                            >  1:961225/1‑150 (MQ=255)
                                                                                                                      acacGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTc                           >  1:452946/1‑150 (MQ=255)
                                                                                                                       ctcggggcAAACATAAAGGCCCATGCAGTTAGCAACCTACCCGGGCGTTTCTCACAATCCCCGCTGAGTTTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGGCGAGGTTGCACAATTGATTATTTTCTGAACATTCa                          <  2:1200060/143‑1 (MQ=255)
                                                                                                                       cacGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCa                          >  2:138682/1‑150 (MQ=255)
                                                                                                                        acGCGCCAAAAATAAAGGCACAGGCAGTTAGCAAACTACCCGGGCCTTTCTCACAATCACCGCTGACATTATATCACCCATCTTTCGGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCaa                         >  2:769370/1‑150 (MQ=255)
                                                                                                                            gcCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCa                     >  2:38505/1‑150 (MQ=255)
                                                                                                                                 tgagaaaGGCCCGGGTAGGTTGCTAACTGCCTGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTcc                >  1:142616/5‑150 (MQ=38)
                                                                                                                                  gagaaaGGCCCGGGTAGGTTGCTAACTGCCTGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTccc               >  1:1591060/4‑150 (MQ=38)
                                                                                                                                    tAAAGGCCCAGGCTGTTTGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACTTTATTTCACCGATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAATTCAACTCccac             <  2:1471027/150‑1 (MQ=255)
                                                                                                                                     aaaGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCccacc            <  1:507041/150‑1 (MQ=255)
                                                                                                                                      aaGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCccacca           >  1:950682/1‑150 (MQ=255)
                                                                                                                                         gCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAg        <  1:156038/150‑1 (MQ=255)
                                                                                                                                          cccAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGc       >  2:1339319/1‑150 (MQ=255)
                                                                                                                                            cAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCca     >  2:1459990/1‑150 (MQ=255)
                                                                                                                                              ggCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCcaac   <  2:1134122/150‑1 (MQ=255)
                                                                                                                                               gCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCcaaca  >  1:1491278/1‑150 (MQ=255)
                                                                                                                                               gCAGGTCGCCACCTACCCCGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCCCCGACGTTGCACAATTGATTAATTTCTGCCCATTCAAAGCACCTCCCCCCAAAGCcaacc  <  2:1457074/150‑2 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACA  >  NC_000913/2563114‑2563406

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