Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,776,132 A→C 6.7% S72R (AGC→CGC)  yfjX → CP4‑57 prophage; putative antirestriction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,776,1320AC 6.7% 271.1 / 10.4 90S72R (AGC→CGC) yfjXCP4‑57 prophage; putative antirestriction protein
Reads supporting (aligned to +/- strand):  ref base A (56/28);  new base C (3/3);  total (59/31)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

CGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCC  >  NC_000913/2775983‑2776279
                                                                                                                                                     |                                                                                                                                                   
cGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCa                                                                                                                                                     >  2:518643/1‑150 (MQ=255)
   aGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCa                                                                                                                                                  <  2:944543/150‑1 (MQ=255)
       aCCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATgg                                                                                                                                              >  2:22164/1‑150 (MQ=255)
          cGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCgg                                                                                                                                           >  1:210650/1‑150 (MQ=255)
                aCAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCCGCCATTGCTGAGCGtg                                                                                                                                     <  2:233144/150‑2 (MQ=255)
                   gCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTAt                                                                                                                                  >  1:606972/1‑150 (MQ=255)
                     atatCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTTATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCCGCCATTGCTGAGCGTGTaaaa                                                                                                                                <  2:1458573/150‑4 (MQ=37)
                      tatCCGCTTCTGTCCGCAGCACTTTGGCACCATCCCGCAGTGTATAACGCTGTAAACGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTTGTCCTTTTACACGCTCCGCCATTGCTGAGCGTGTaaaag                                                                                                                               <  2:1001048/150‑5 (MQ=37)
                      tatCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGc                                                                                                                               >  2:766187/1‑150 (MQ=255)
                       atCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGCGCGTTTATGGcg                                                                                                                              >  2:1428115/1‑149 (MQ=255)
                       atCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCt                                                                                                                              >  1:713584/1‑150 (MQ=255)
                       atCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCt                                                                                                                              >  2:760076/1‑150 (MQ=255)
                         ccGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTcc                                                                                                                            <  2:568060/150‑1 (MQ=255)
                          cGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCt                                                                                                                           <  2:1467398/150‑1 (MQ=255)
                              tCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATg                                                                                                                       >  2:1331657/1‑150 (MQ=255)
                               cTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGc                                                                                                                      >  2:1111244/1‑150 (MQ=255)
                               cTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGc                                                                                                                      <  1:1002521/150‑1 (MQ=255)
                                      cagcaCTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGt                                                                                                               <  2:919464/150‑1 (MQ=255)
                                      cagcaCTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGt                                                                                                               <  2:213589/150‑1 (MQ=255)
                                         caCTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGAc                                                                                                            <  1:1329444/150‑1 (MQ=255)
                                         caCTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGAc                                                                                                            >  1:39657/1‑150 (MQ=255)
                                                gCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAAt                                                                                                     >  1:1431381/1‑150 (MQ=255)
                                                 cACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCg                                                                                                                                       >  2:559417/1‑115 (MQ=255)
                                                 cACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCg                                                                                                                                       <  1:559417/115‑1 (MQ=255)
                                                 cACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATg                                                                                                    <  2:683174/150‑1 (MQ=255)
                                                 cACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATg                                                                                                    >  1:475449/1‑150 (MQ=255)
                                                   ccATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGc                                                                                                  <  2:713584/150‑1 (MQ=255)
                                                   ccATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGc                                                                                                  >  1:1156472/1‑150 (MQ=255)
                                                      tCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCAtt                                                                                               >  1:788834/1‑150 (MQ=255)
                                                         cGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGc                                                                                            <  1:1111244/150‑1 (MQ=255)
                                                          gCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCt                                                                                           >  2:981504/1‑150 (MQ=255)
                                                           cAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTc                                                                                          <  1:283633/150‑1 (MQ=255)
                                                            aGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCa                                                                                         >  1:508457/1‑150 (MQ=255)
                                                              tGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGGGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAAc                                                                                       >  1:1492971/1‑150 (MQ=255)
                                                               ggATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACg                                                                                      >  2:16045/1‑150 (MQ=255)
                                                               ggATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCGCAAcc                                                                                      >  2:683021/1‑149 (MQ=255)
                                                                 aTAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGc                                                                                    >  1:481189/1‑150 (MQ=255)
                                                                  tAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGAc                                                                                                                  >  2:1082825/1‑119 (MQ=255)
                                                                  tAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGAc                                                                                                                  <  1:1082825/119‑1 (MQ=255)
                                                                  tAACGCTCGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGCTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCa                                                                                   >  2:43818/1‑150 (MQ=255)
                                                                    aCGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATg                                                                                 >  2:246951/1‑150 (MQ=255)
                                                                       cTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAAt                                                                              >  1:1453546/1‑150 (MQ=255)
                                                                            aCCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCaa                                                                         >  1:911967/1‑150 (MQ=255)
                                                                            aCCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCaa                                                                         >  1:566782/1‑150 (MQ=255)
                                                                              cgcgCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATg                                                                       >  2:1294687/1‑150 (MQ=255)
                                                                              cgcgCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATg                                                                       >  1:900685/1‑150 (MQ=255)
                                                                                    ttttCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCgg                                                                 >  2:582462/1‑150 (MQ=255)
                                                                                    ttttCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCgg                                                                 >  2:1374133/1‑150 (MQ=255)
                                                                                            tggatggaCCGCTTCTGTGATGAGTACTTCGGGGGTATCTGGTCCTTTTACACCCTCAGCAATGGCGGAGCGTTTATCGCTCCTGATGCTGACGGTGACGATAAATGGCATTTCCTCAACGGCATGAATGTCAATCGTGACGAAATAAgc                                                         >  2:1156381/1‑150 (MQ=37)
                                                                                             ggatggaCCGCTTCTGGGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAgcg                                                        <  1:747292/150‑1 (MQ=255)
                                                                                              gatggaCCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAgcgc                                                       >  1:203495/1‑150 (MQ=255)
                                                                                                tggaCCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCgg                                                     >  1:754289/1‑150 (MQ=255)
                                                                                                   aCCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAg                                                  >  2:902711/1‑150 (MQ=255)
                                                                                                      gCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCg                                               >  2:571756/1‑150 (MQ=255)
                                                                                                       cTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGc                                              <  2:754289/150‑1 (MQ=255)
                                                                                                        ttCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCa                                             <  2:1302466/150‑1 (MQ=255)
                                                                                                         tCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAg                                            >  1:724232/1‑150 (MQ=255)
                                                                                                            gtgATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTa                                         >  1:1114685/1‑150 (MQ=255)
                                                                                                            gtgATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTa                                         <  2:1431381/150‑1 (MQ=255)
                                                                                                             tgatgaGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTAt                                        <  1:766187/150‑1 (MQ=255)
                                                                                                                   gTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGt                                  >  2:742543/1‑150 (MQ=255)
                                                                                                                      cTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTg                               >  2:56166/1‑150 (MQ=255)
                                                                                                                       tGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGc                              >  1:712627/1‑150 (MQ=255)
                                                                                                                        gCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGcc                             <  2:1453546/150‑1 (MQ=255)
                                                                                                                             ggTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGAtt                        <  2:210650/150‑1 (MQ=255)
                                                                                                                             ggTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGAtt                        <  1:16045/150‑1 (MQ=255)
                                                                                                                               tATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGa                      >  1:508930/1‑150 (MQ=255)
                                                                                                                                 tCTGGTCCTTTTACACGCGCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAat                    <  2:1219877/150‑1 (MQ=255)
                                                                                                                                  cTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAata                   >  2:1408473/1‑150 (MQ=255)
                                                                                                                                  cTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAata                   <  2:788834/150‑1 (MQ=255)
                                                                                                                                   tGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAatat                  >  2:714439/1‑150 (MQ=255)
                                                                                                                                    ggTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAAtata                 >  2:551799/1‑150 (MQ=255)
                                                                                                                                    agTCCTTTTAAAAGCTCAACAATTGCGGATCTTTTATGGATACTGATGCTGACGTTTAATATAAATAGAATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAAtata                 >  2:735328/2‑150 (MQ=255)
                                                                                                                                        gccataaacGCTCCGCCATTGCTGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCa             >  2:896955/8‑150 (MQ=255)
                                                                                                                                        cTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCa             >  2:301211/1‑150 (MQ=255)
                                                                                                                                         ttttACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCAt            <  1:474756/150‑1 (MQ=255)
                                                                                                                                          tttACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATc           >  2:1032625/1‑150 (MQ=255)
                                                                                                                                          tttACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATc           >  1:507035/1‑150 (MQ=255)
                                                                                                                                          cataaacGCTCCGCCATTGCTGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATc           >  1:1193294/6‑150 (MQ=255)
                                                                                                                                          cataaacGCTCCGCCATTGCTGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATc           >  1:1531510/6‑150 (MQ=255)
                                                                                                                                             acacGCTCAGCAATTGCTGAGCGTTTATATCTCCTGATTCTCACTATCAAGAAAAAAAAAATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCAGCCTGATTGAATATAGACATCAcc        >  2:1179007/1‑149 (MQ=255)
                                                                                                                                             acacGCTCAGCAATGGCGGAGCGTTTATGTCTCCTGATTCTGACGGTCAAGATAAAAAGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACg        >  2:421698/1‑150 (MQ=255)
                                                                                                                                             acacGCTCAGCAATGGCGGAGCGTTTATGTCTCCTGATGCTGACGTTGACGATAAATGTCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACg        >  2:1050678/1‑150 (MQ=255)
                                                                                                                                             acacGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACg        <  1:325927/150‑1 (MQ=255)
                                                                                                                                             acacGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACg        <  2:39657/150‑1 (MQ=255)
                                                                                                                                             acacGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACCGTGACGATAAATGCCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACg        >  2:385136/1‑150 (MQ=255)
                                                                                                                                                cGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcct     >  1:1509148/1‑150 (MQ=255)
                                                                                                                                                cGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcct     >  1:1353432/1‑150 (MQ=255)
                                                                                                                                                cGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcct     >  1:304849/1‑150 (MQ=255)
                                                                                                                                                  cTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcctgc   <  1:990875/150‑1 (MQ=255)
                                                                                                                                                   tCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCATCGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcctgcc  <  1:518643/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCC  >  NC_000913/2775983‑2776279

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