Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,770,310 T→A 7.2% V203D (GTT→GAT) ‡ 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,3100TA 7.2% 284.1 / 10.4 98V203D (GTT→GAT) ‡yfjRCP4‑57 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (39/51);  new base A (3/4);  total (43/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.14e-01

TCCCTGCCGGTTACTCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTG  >  NC_000913/2770164‑2770457
                                                                                                                                                  |                                                                                                                                                   
tCCCTGCCGGTTACTCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGa                                                                                                                                                  <  2:180939/150‑1 (MQ=255)
        ggTTACTCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCgg                                                                                                                                          <  1:769919/150‑1 (MQ=255)
          ttACTCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTc                                                                                                                                        <  2:2340/150‑1 (MQ=255)
          ttACTCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTc                                                                                                                                        <  1:1087810/150‑1 (MQ=255)
              tCATCCACCAGCAGACCTTGTATCTGGTGGCTGAACACGAAGGGAATATCGCCGTATTCACACTTTATGAAATCAATCTGATTCAGCCTTTGCAGGAGACTTTTAGCCTCAATGACATTCTGTGTCGCCTGGTTGAAGACCCGGTCttaa                                                                                                                                    <  2:266993/150‑3 (MQ=255)
              tCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttca                                                                                                                                    <  1:962314/150‑1 (MQ=255)
              tCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttca                                                                                                                                    <  1:1235094/150‑1 (MQ=255)
               cATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttcat                                                                                                                                   >  2:441748/1‑150 (MQ=255)
               cATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttcat                                                                                                                                   <  1:981280/150‑1 (MQ=255)
               cATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttcat                                                                                                                                   >  2:267566/1‑150 (MQ=255)
               cATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttcat                                                                                                                                   >  1:1250021/1‑150 (MQ=255)
               cATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCGGGGTGGAGGCCggggggtgcc                                                                                                                                   >  1:641029/1‑143 (MQ=255)
                 tCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCGGGTCttcaacc                                                                                                                                 <  1:811175/150‑4 (MQ=255)
                 tCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttcattc                                                                                                                                 >  1:1386587/1‑150 (MQ=255)
                 tCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttcattc                                                                                                                                 >  2:1373646/1‑150 (MQ=255)
                  ccaccaGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttcattca                                                                                                                                <  2:371041/150‑1 (MQ=255)
                    accaGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCGGGTCTTCAACCAgg                                                                                                                              <  2:165363/150‑1 (MQ=255)
                     ccaGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGc                                                                                                                             >  1:1010762/1‑150 (MQ=255)
                     ccaGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGc                                                                                                                             <  2:1194646/150‑1 (MQ=255)
                     ccaGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGAATGAAGACCGGGTCTTCAACCAGGc                                                                                                                             <  1:182926/150‑1 (MQ=255)
                     ccaGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGAATGAAGACCGGGTCTTCAACCAGGc                                                                                                                             <  1:19940/150‑1 (MQ=255)
                      cagcagACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGcc                                                                                                                            <  2:418005/150‑1 (MQ=255)
                      cagcagACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGcc                                                                                                                            <  2:412776/150‑1 (MQ=255)
                      cagcagACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGAATGAAGACCGGGTCTTCAACCAGGcg                                                                                                                            <  1:715000/150‑2 (MQ=255)
                        gcagACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGAATGAAGACCGGGTCTTCAACCAGGcgac                                                                                                                          <  1:1298342/150‑4 (MQ=255)
                           gACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTAcc                                                                                                                       <  2:221835/150‑1 (MQ=255)
                           gACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTAcc                                                                                                                       >  2:1275740/1‑150 (MQ=255)
                              cTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCa                                                                                                                    >  2:1250420/1‑150 (MQ=255)
                                ggTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCAtt                                                                                                                  <  1:1373646/150‑1 (MQ=255)
                                    tCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCg                                                                                                              <  2:1251465/150‑1 (MQ=255)
                                     cTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGc                                                                                                             <  2:35655/150‑1 (MQ=255)
                                      tggtggCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCt                                                                                                            <  2:1121894/150‑1 (MQ=255)
                                      aggtggCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCt                                                                                                            >  2:721809/2‑150 (MQ=255)
                                       ggtggCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCtt                                                                                                           <  2:321508/150‑1 (MQ=255)
                                       ggtggCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCtt                                                                                                           >  1:456964/1‑150 (MQ=255)
                                       ggtggCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCtt                                                                                                           >  1:243587/1‑150 (MQ=255)
                                       ggtggCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCtt                                                                                                           <  1:573854/150‑1 (MQ=255)
                                         tggCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTa                                                                                                         <  1:267566/150‑1 (MQ=255)
                                         tggCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTa                                                                                                         <  1:441748/150‑1 (MQ=255)
                                         tggCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTa                                                                                                         <  2:591138/150‑1 (MQ=255)
                                          ggCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTAt                                                                                                        >  2:927107/1‑150 (MQ=255)
                                          ggCTGAACACGAAGGGCATATCGCCGTATGCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTAt                                                                                                        >  1:393988/1‑150 (MQ=255)
                                            cTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATcc                                                                                                      <  1:1007143/150‑1 (MQ=255)
                                             tGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAgg                                                                                                                              <  2:1456699/125‑1 (MQ=255)
                                             tGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAgg                                                                                                                              >  1:1456699/1‑125 (MQ=255)
                                             tGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCtt                                                                                                           >  2:949759/1‑144 (MQ=255)
                                             tGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCtt                                                                                                           <  1:949759/144‑1 (MQ=255)
                                             tGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCa                                                                                                     >  1:1296222/1‑150 (MQ=255)
                                              gAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAg                                                                                                    >  2:1074918/1‑150 (MQ=255)
                                                acacGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCa                                                                                                  >  1:767121/1‑150 (MQ=255)
                                                acacGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCa                                                                                                  >  1:371947/1‑150 (MQ=255)
                                                    gAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCa                                                                                              >  2:708834/1‑150 (MQ=255)
                                                    gAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCa                                                                                              >  1:403166/1‑150 (MQ=255)
                                                          cATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCtt                                                                                        <  2:1010762/150‑1 (MQ=255)
                                                          cATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCtt                                                                                        <  2:125926/150‑1 (MQ=255)
                                                                 ccGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTg                                                                                 <  2:343205/150‑1 (MQ=255)
                                                                 ccGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTg                                                                                 <  1:441037/150‑1 (MQ=255)
                                                                      ttCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCg                                                                            <  2:1295088/150‑1 (MQ=255)
                                                                      ttCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCg                                                                            >  2:1179720/1‑150 (MQ=255)
                                                                        cacaCTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGGCGCGTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGc                                                                          <  1:114760/150‑1 (MQ=255)
                                                                          cacTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCg                                                                        >  1:48560/1‑150 (MQ=255)
                                                                               gatgaAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCt                                                                                         <  1:1102351/128‑1 (MQ=255)
                                                                               gatgaAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCt                                                                                         >  2:1102351/1‑128 (MQ=255)
                                                                               gatgaAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCATGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCAATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGc                                                                   >  1:1009016/1‑150 (MQ=255)
                                                                                    aaTCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCAcc                                                              >  1:1182021/1‑150 (MQ=255)
                                                                                    aaTCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATCACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCAcc                                                              >  1:1157149/1‑150 (MQ=255)
                                                                                         aTCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCGGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAAt                                                         >  2:1359370/1‑150 (MQ=255)
                                                                                          tCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATa                                                        >  1:291109/1‑150 (MQ=255)
                                                                                            tGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCAGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAgc                                                      <  2:835890/150‑1 (MQ=255)
                                                                                                     cTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCgtt                                             >  1:750270/1‑150 (MQ=255)
                                                                                                     cTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCgtt                                             <  2:1182021/150‑1 (MQ=255)
                                                                                                     cTTTGCAGGAGACTTTTCGCCGCAATCACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCgtt                                             <  2:1157149/150‑1 (MQ=255)
                                                                                                       ttGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCgttgt                                           <  2:882945/150‑1 (MQ=255)
                                                                                                           aGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATc                                       <  2:1296222/150‑1 (MQ=255)
                                                                                                            ggAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCa                                      >  1:589957/1‑150 (MQ=255)
                                                                                                               gaCTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAAcc                                   >  2:110310/1‑150 (MQ=255)
                                                                                                                 cTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCg                                 >  1:233028/1‑150 (MQ=255)
                                                                                                                 cTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCg                                 >  1:1029245/1‑150 (MQ=255)
                                                                                                                   tttCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGc                               <  1:614729/150‑1 (MQ=255)
                                                                                                                    ttCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCa                              <  1:338164/150‑1 (MQ=255)
                                                                                                                         cgcAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATg                         <  1:721809/150‑1 (MQ=255)
                                                                                                                         cgcAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATg                         >  1:112658/1‑150 (MQ=255)
                                                                                                                          gcAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATgg                        <  1:499165/150‑1 (MQ=255)
                                                                                                                           cAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATgga                       <  1:1074918/150‑1 (MQ=255)
                                                                                                                                 cAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGcc                 <  2:1029245/150‑1 (MQ=255)
                                                                                                                                  aGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGccc                <  2:1486084/150‑1 (MQ=255)
                                                                                                                                    tCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAg              <  1:397102/150‑1 (MQ=255)
                                                                                                                                     gggtaaggCCTGAATGAAGACCGGGTCTTCAACCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGt             >  2:184217/8‑150 (MQ=38)
                                                                                                                                      ggtaaggCCTGAATGAAGACCGGGTCTTCAACCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTg            >  1:557646/7‑150 (MQ=39)
                                                                                                                                       gtgtCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGc           >  1:143295/1‑150 (MQ=255)
                                                                                                                                         gtCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGccc         <  1:1140321/150‑1 (MQ=255)
                                                                                                                                         gtCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGccc         >  1:892535/1‑150 (MQ=255)
                                                                                                                                         gtCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGccc         >  2:1408348/1‑150 (MQ=255)
                                                                                                                                          tCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCg        <  2:747072/150‑1 (MQ=255)
                                                                                                                                          aggCCTGAATGAAGACCGGGTCTTCAACCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCg        >  1:1126606/3‑150 (MQ=255)
                                                                                                                                            gCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGtt      <  2:371947/150‑1 (MQ=255)
                                                                                                                                               tGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTAtt   >  1:259507/1‑150 (MQ=255)
                                                                                                                                                ggTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTg  <  1:110310/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCCCTGCCGGTTACTCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTG  >  NC_000913/2770164‑2770457

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