Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,770,319 C→G 6.6% P206R (CCG→CGG)  yfjR → CP4‑57 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,770,3190CG 6.6% 422.0 / 16.6 133P206R (CCG→CGG) yfjRCP4‑57 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (79/45);  new base G (4/5);  total (83/50)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 139 total aligned reads displayed.

GGTTACTCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCC  >  NC_000913/2770172‑2770468
                                                                                                                                                   |                                                                                                                                                     
ggTTACTCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCgg                                                                                                                                                     >  1:1221547/1‑150 (MQ=255)
   tACTCATCCTCCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCt                                                                                                                                                  <  1:80219/150‑1 (MQ=255)
   tACTCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGCCTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCt                                                                                                                                                  >  1:1566515/1‑150 (MQ=255)
    aCTCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCtt                                                                                                                                                 <  1:908761/150‑1 (MQ=255)
     cTCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttc                                                                                                                                                >  1:1765505/1‑150 (MQ=255)
       cATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttcat                                                                                                                                              <  2:1712662/150‑1 (MQ=255)
         tCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttcattc                                                                                                                                            <  1:971514/150‑1 (MQ=255)
          ccaccaGCAGCCCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGCCCCGGTCttcattca                                                                                                                                           >  2:1000051/1‑150 (MQ=255)
           caccaGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAg                                                                                                                                          >  2:707137/1‑150 (MQ=255)
            accaGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCGGGTCTTCAACCAgg                                                                                                                                         <  1:1699216/150‑1 (MQ=255)
             ccaGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACCGTCTGTGTCGCCTGGTTGAAGACCGGGTCTTCAACCAGGc                                                                                                                                        <  1:817501/150‑1 (MQ=255)
              cagcagACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGcc                                                                                                                                       <  2:160762/150‑1 (MQ=255)
              cagcagACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACATGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGcc                                                                                                                                       <  2:1899395/150‑1 (MQ=255)
                   gACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTAcc                                                                                                                                  <  1:950296/150‑1 (MQ=255)
                    aCCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTAccc                                                                                                                                 <  1:59327/150‑1 (MQ=255)
                    aCCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGAATGAAGACCGGGTCTTCAACCAGGcgacacag                                                                                                                                 <  1:646421/150‑8 (MQ=38)
                       tGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCAt                                                                                                                              >  1:456294/1‑150 (MQ=255)
                         gTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCAttt                                                                                                                            <  2:839088/150‑1 (MQ=255)
                           aTCTGGTGGCTGAACACGAAGGGCCTATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGTCCCGGTCTTCATTCAGGCCTTACCCCATTTTc                                                                                                                          <  1:34327/150‑1 (MQ=255)
                           aTCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTc                                                                                                                          >  2:633326/1‑150 (MQ=255)
                             cTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGc                                                                                                                        >  1:523942/1‑150 (MQ=255)
                               gggggCTGAACACGCAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGGCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCtt                                                                                                                      <  2:1350331/147‑1 (MQ=255)
                                  ggCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTAt                                                                                                                   >  1:1546555/1‑150 (MQ=255)
                                  ggCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTAt                                                                                                                   <  1:422666/150‑1 (MQ=255)
                                  ggCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTAt                                                                                                                   >  2:1039829/1‑150 (MQ=255)
                                    cTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGATTTATcc                                                                                                                 <  2:1380884/150‑1 (MQ=255)
                                     tGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCTTTTTCGCTTTATCCa                                                                                                                >  1:992346/1‑150 (MQ=255)
                                        acacGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCCTTCAGGCCTTACCCCATTTTCGCTTTATCCAGCa                                                                                                             >  2:1066720/1‑150 (MQ=255)
                                           cGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTc                                                                                                          >  2:1319861/1‑150 (MQ=255)
                                                ggCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGc                                                                                                     <  1:1404678/150‑1 (MQ=255)
                                                   atatCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTa                                                                                                  <  2:651388/150‑1 (MQ=255)
                                                    tatCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTAc                                                                                                 >  1:167472/1‑150 (MQ=255)
                                                    tatCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTAc                                                                                                 >  1:145072/1‑150 (MQ=255)
                                                     atCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACg                                                                                                >  1:492297/1‑150 (MQ=255)
                                                     atCGCCGTATTCACACTCGATGAAATCCATCTGATACAGCCTTTGCAGGAGACTTTTCGCCCCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACACCATTTTCGATATATCCAGCATTCACTGCTTAct                                                                                                >  2:1595923/1‑149 (MQ=255)
                                                       cGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGtt                                                                                              >  2:961077/1‑150 (MQ=255)
                                                       cGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGtt                                                                                              <  1:1491100/150‑1 (MQ=255)
                                                        gCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGttt                                                                                             >  1:1796927/1‑150 (MQ=255)
                                                         ccGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTg                                                                                            <  2:273360/150‑1 (MQ=255)
                                                           gTATTCACACTTTATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGtt                                                                                          <  1:1381522/150‑1 (MQ=255)
                                                            tATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGGCTTCATTCAGGCCGTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTc                                                                                         >  2:485259/1‑150 (MQ=255)
                                                             aTTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGCCCCGGGCTTCCTTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTcc                                                                                        >  2:167364/1‑150 (MQ=255)
                                                             aTTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCGGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTcc                                                                                        <  2:1564870/150‑1 (MQ=255)
                                                                cacaCTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGc                                                                                     >  2:1571693/1‑150 (MQ=255)
                                                                cacaCTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGc                                                                                     <  1:1273368/150‑1 (MQ=255)
                                                                 acacTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGcc                                                                                    <  1:73083/150‑1 (MQ=255)
                                                                 acacTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGcc                                                                                    <  1:2027819/150‑1 (MQ=255)
                                                                  cacTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCGTTATCCAGCATTCACTGCTTACGGTTGTTCCGGCCg                                                                                   >  2:1110012/1‑150 (MQ=255)
                                                                     ttGATGAAATCCATCTGATTAAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACa                                                                                <  2:1039002/150‑1 (MQ=255)
                                                                       gatgaAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGTCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGc                                                                              >  2:1527820/1‑150 (MQ=255)
                                                                        atgaAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGcc                                                                             >  2:1612492/1‑150 (MQ=255)
                                                                        atgaAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGcc                                                                             <  2:1159848/150‑1 (MQ=255)
                                                                           aaaTCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGACTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTGAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCAc                                                                          <  1:807530/150‑1 (MQ=255)
                                                                              tCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTGGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCgg                                                                       <  1:1466496/150‑1 (MQ=255)
                                                                               ccATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGa                                                                      >  1:1122504/1‑150 (MQ=255)
                                                                                 aTCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAAt                                                                    >  2:1960427/1‑150 (MQ=255)
                                                                                    tGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAgc                                                                 >  2:1902234/1‑150 (MQ=255)
                                                                                     gATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAgcg                                                                >  1:620254/1‑150 (MQ=255)
                                                                                       ttCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACGCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGACCACCGGAATAGCGCa                                                              >  1:821843/1‑150 (MQ=255)
                                                                                       ttCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCa                                                              >  1:912733/1‑150 (MQ=255)
                                                                                          aGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGc                                                           >  2:228863/1‑150 (MQ=255)
                                                                                          aGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGc                                                           >  2:1396424/1‑150 (MQ=255)
                                                                                          aGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGc                                                           >  1:818084/1‑150 (MQ=255)
                                                                                          aGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACACGCTCTTCATTCACGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCAACCGGAATAGCGCAGGc                                                           >  2:1414130/1‑150 (MQ=255)
                                                                                           gCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCg                                                          <  1:1710604/150‑1 (MQ=255)
                                                                                           gCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCg                                                          >  1:461614/1‑150 (MQ=255)
                                                                                             cTTTGCAGGAGACTTTTCGCCGCAATGCCAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGGTTACGATTGTTCCGGCCGACAGCCCACCGCAATTCCGCAGGCgtt                                                        >  2:148050/1‑150 (MQ=255)
                                                                                             cTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTAAGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCgtt                                                        >  2:866775/1‑150 (MQ=255)
                                                                                                tGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCgttgtt                                                     >  1:1423183/1‑150 (MQ=255)
                                                                                                   aGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATc                                                  <  2:1163391/150‑1 (MQ=255)
                                                                                                    ggAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCCGCATTCACTGGTTCCGTGTGTTCCGGCCGACAGCCCACAGGAATAGAGAAGGCGTTGTTATCa                                                 >  2:1246396/1‑150 (MQ=255)
                                                                                                    ggAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTAGCTTTATCCAGCATTCACTGCTTACGTATGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCgttgtaagca                                                 >  1:210816/1‑145 (MQ=255)
                                                                                                         cTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTAGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGTCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCg                                            >  2:1097698/1‑150 (MQ=255)
                                                                                                           tttCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGc                                          <  1:1883303/150‑1 (MQ=255)
                                                                                                            ttCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCa                                         <  1:1865483/150‑1 (MQ=255)
                                                                                                               gccgcAATGACAGTCTGTGTCGCCTGGTTGAAGACACGGTCTTCATTCAGGCCTTACCCCATTTTCGCTGTATCCAGCAGTCACTGCTTACGTTTGTTCCGGGCGTCAGCCCACCGGATTAGCGCAGGCGGTGATATCAACCCGGAAACa                                      >  2:360241/1‑150 (MQ=255)
                                                                                                                  gcAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATgg                                   <  1:780858/150‑1 (MQ=255)
                                                                                                                         cAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTAATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGcc                            <  2:785720/150‑1 (MQ=255)
                                                                                                                           gTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCCGGCGTTGTTATCAACCCGGCAACATGGAGGCGCCCa                          >  1:1008152/1‑150 (MQ=255)
                                                                                                                           gTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCa                          <  1:1607136/150‑1 (MQ=255)
                                                                                                                                 gtCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGccc                    >  2:565150/1‑150 (MQ=255)
                                                                                                                                    gCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCCCTTTATCCAGCAGTCACTGCTTACGTTTGTTCCGGCCGCCAGCCCACCGGAATAGAGCAGGCGTTGTTATCAACCCGGCAACATGTAGGAGCCCAGTGCCAGtt                 >  1:900900/1‑150 (MQ=255)
                                                                                                                                     ccTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTa                >  2:495517/1‑150 (MQ=255)
                                                                                                                                     ccTGGTTGAAGACCCGGTCTTCATTCAGGCATTACCCCAATTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGATATCAACCCCGCAACATGGAGGAGCCCAGTGCCCGTTa                >  1:1514039/1‑150 (MQ=255)
                                                                                                                                      cTGGTTGAACACCCGGTCTTCATTCAGGCCTTACCACATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTGCCGGCCGACAGCCCACCGGAATAGCGCAGGAGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTAt               <  1:344936/150‑1 (MQ=255)
                                                                                                                                        ggTTGACGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTg             <  1:93104/150‑1 (MQ=255)
                                                                                                                                        ggTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACCGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTg             >  1:1773391/1‑150 (MQ=255)
                                                                                                                                        ggTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTg             >  2:549659/1‑150 (MQ=255)
                                                                                                                                        ggTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTATCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCCGTGCCCGTTATTg             >  1:1318671/1‑150 (MQ=255)
                                                                                                                                        gaatGAAGCCCGGGTCTTCATTCAGGCATTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCATCATGGAGGAGCCCAGTGCCCGTTATTg             >  2:876190/4‑150 (MQ=255)
                                                                                                                                         gTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACAGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGc            >  2:624620/1‑150 (MQ=255)
                                                                                                                                          ttGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTAATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGcc           <  2:93881/150‑1 (MQ=255)
                                                                                                                                           tGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCa          >  1:271215/1‑150 (MQ=255)
                                                                                                                                           tGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCACAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGACCGTTATTGCCa          >  1:885671/1‑150 (MQ=255)
                                                                                                                                             aaGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCGACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCAtt        >  1:1294461/1‑150 (MQ=255)
                                                                                                                                             aaGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCACCATGGAGGCGCCCAGTGCCCGTTATTGCCAtt        >  2:1796472/1‑150 (MQ=255)
                                                                                                                                               gACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTAt      >  1:700780/1‑150 (MQ=255)
                                                                                                                                                 cccGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCg    >  1:1094477/1‑150 (MQ=255)
                                                                                                                                                 cccGGTCTTCATTCAGGACTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCg    >  1:1149842/1‑150 (MQ=255)
                                                                                                                                                   cGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGcc  <  1:961078/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGTTACTCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCC  >  NC_000913/2770172‑2770468

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: