Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,764,549 T→C D202G (GAC→GGC)  abpA ← CP4‑57 prophage; anti‑bacteriophage protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,764,5490TC98.9% 322.4 / ‑4.6 94D202G (GAC→GGC) abpACP4‑57 prophage; anti‑bacteriophage protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (46/47);  total (46/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

TGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAAT  >  NC_000913/2764417‑2764681
                                                                                                                                    |                                                                                                                                    
tGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGc                                                                                                                               <  2:589251/140‑1 (MQ=255)
 ggTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCg                                                                                                                              <  2:691110/140‑1 (MQ=255)
  gTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGa                                                                                                                             >  1:1057130/1‑140 (MQ=255)
  gTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGa                                                                                                                             >  2:113997/1‑140 (MQ=255)
  gTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGa                                                                                                                             >  1:829121/1‑140 (MQ=255)
       cTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGtt                                                                                                                        >  2:458079/1‑140 (MQ=255)
       cTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGtt                                                                                                                        >  1:16411/1‑140 (MQ=255)
       cTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGtt                                                                                                                        >  1:553222/1‑140 (MQ=255)
        ttCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGttt                                                                                                                       >  2:320233/1‑140 (MQ=255)
           gACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCgg                                                                                                                    <  2:553222/140‑1 (MQ=255)
            acacGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGc                                                                                                                               <  2:850293/128‑1 (MQ=255)
            acacGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGc                                                                                                                               >  1:850293/1‑128 (MQ=255)
             cacGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGaa                                                                                                                  <  1:159253/140‑1 (MQ=255)
                 ctgctTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTga                                                                                                              >  2:1042704/1‑140 (MQ=255)
                 ctgctTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTga                                                                                                              >  2:1076745/1‑140 (MQ=255)
                 ctgctTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTga                                                                                                              >  1:125804/1‑140 (MQ=255)
                    ctTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGGCGGTTTCGGAACTGAGaa                                                                                                           >  1:415952/1‑140 (MQ=255)
                    ctTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGaa                                                                                                           <  1:265246/140‑1 (MQ=255)
                    ctTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGaa                                                                                                           <  2:1079948/140‑1 (MQ=255)
                    ctTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGaa                                                                                                           <  2:125804/140‑1 (MQ=255)
                      tGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATg                                                                                                         <  2:829121/140‑1 (MQ=255)
                      tGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATg                                                                                                         <  2:1057130/140‑1 (MQ=255)
                      tGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATg                                                                                                         <  1:320233/140‑1 (MQ=255)
                      tGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATg                                                                                                         <  1:1076745/140‑1 (MQ=255)
                      tGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATg                                                                                                         <  2:491839/140‑1 (MQ=255)
                       gAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGa                                                                                                        >  1:516551/1‑140 (MQ=255)
                       gAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGa                                                                                                        <  1:634056/140‑1 (MQ=255)
                       gAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGa                                                                                                        <  1:458079/140‑1 (MQ=255)
                        aaCCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGAc                                                                                                       <  1:849494/140‑1 (MQ=255)
                        aaCCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGAc                                                                                                       <  1:309382/140‑1 (MQ=255)
                         aCCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGCCGGTTTCGGAACTGAGAATGAcc                                                                                                      <  2:761093/140‑1 (MQ=255)
                         aCCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGAcc                                                                                                      <  1:113997/140‑1 (MQ=255)
                         aCCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGAcc                                                                                                      <  2:516551/140‑1 (MQ=255)
                          ccACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGAcc                                                                                                      >  2:271312/1‑139 (MQ=255)
                          ccACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGAcc                                                                                                      <  1:271312/139‑1 (MQ=255)
                            aCTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGAc                                                                                                       <  1:685199/136‑1 (MQ=255)
                            aCTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGAc                                                                                                       >  2:685199/1‑136 (MQ=255)
                              tCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATg                                                                                                 >  2:32099/1‑140 (MQ=255)
                              tCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATg                                                                                                 <  2:577202/140‑1 (MQ=255)
                              tCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATg                                                                                                 >  2:293474/1‑140 (MQ=255)
                              tCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATg                                                                                                 >  1:380218/1‑140 (MQ=255)
                              tCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATg                                                                                                 >  2:363197/1‑140 (MQ=255)
                              tCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATg                                                                                                 <  1:293474/140‑1 (MQ=255)
                                     ggTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGaa                                                                                                           >  2:1008906/1‑123 (MQ=255)
                                     ggTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGaa                                                                                                           <  1:1008906/123‑1 (MQ=255)
                                          ggTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAg                                                                                              >  2:193587/1‑131 (MQ=255)
                                          ggTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAg                                                                                              <  1:193587/131‑1 (MQ=255)
                                          ggTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGa                                                                                     <  2:380218/140‑1 (MQ=255)
                                          ggTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGa                                                                                     >  1:202776/1‑140 (MQ=255)
                                           gTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAg                                                                                    >  1:18575/1‑140 (MQ=255)
                                            tCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAgg                                                                                   <  1:32099/140‑1 (MQ=255)
                                                aaGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATg                                                                                                 >  1:619342/1‑122 (MQ=255)
                                                aaGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATg                                                                                                 <  2:619342/122‑1 (MQ=255)
                                                aaGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAg                                                                               <  1:882171/140‑1 (MQ=255)
                                                    aGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGt                                                                           <  2:18575/140‑1 (MQ=255)
                                                     ggACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTa                                                                          <  2:16411/140‑1 (MQ=255)
                                                      gACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCa                                                                                >  1:345750/1‑133 (MQ=255)
                                                      gACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCa                                                                                <  2:345750/133‑1 (MQ=255)
                                                        cTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGt                                                                       >  1:690841/1‑140 (MQ=255)
                                                           gTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCt                                                                    <  1:340194/140‑1 (MQ=255)
                                                           gTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCt                                                                    <  2:690841/140‑1 (MQ=255)
                                                           gTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCt                                                                    <  2:165768/140‑1 (MQ=255)
                                                            tcatcaGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGGCGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCtt                                                                   <  2:415952/140‑1 (MQ=255)
                                                                caGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATc                                                               >  2:645562/1‑140 (MQ=255)
                                                                        tgAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAg                                                       >  2:1121373/1‑140 (MQ=255)
                                                                         gAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGt                                                      >  1:812258/1‑140 (MQ=255)
                                                                             gCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCa                                                  >  2:1141669/1‑140 (MQ=255)
                                                                             gCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCa                                                  >  2:662732/1‑140 (MQ=255)
                                                                             gCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCa                                                  >  2:131451/1‑140 (MQ=255)
                                                                               gTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGt                                                <  1:131451/140‑1 (MQ=255)
                                                                               gTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGt                                                <  1:443112/140‑1 (MQ=255)
                                                                                tCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGtt                                               >  1:341679/1‑140 (MQ=255)
                                                                                  gagaATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTcg                                             >  2:729375/1‑140 (MQ=255)
                                                                                        gCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAgg                                       <  1:363197/140‑1 (MQ=255)
                                                                                           cGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGAt                                    <  1:645562/140‑1 (MQ=255)
                                                                                                ttGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGaaaa                               <  2:812258/140‑1 (MQ=255)
                                                                                                 tGCGCGTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAac                              <  1:729375/140‑1 (MQ=255)
                                                                                                     cgTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACtt                          >  2:678740/1‑140 (MQ=255)
                                                                                                      gTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAac                              <  2:212266/135‑1 (MQ=255)
                                                                                                      gTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAac                              >  1:212266/1‑135 (MQ=255)
                                                                                                      gTCACAGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTc                         <  1:662732/140‑1 (MQ=255)
                                                                                                        cacaGTATCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTccc                       >  2:333981/1‑140 (MQ=255)
                                                                                                               tCAAGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGc                <  1:626443/140‑1 (MQ=255)
                                                                                                                  aGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCGTGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGc             >  2:777953/1‑140 (MQ=255)
                                                                                                                  aGGATGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGc             >  2:495625/1‑140 (MQ=255)
                                                                                                                     aTGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTcc                        >  1:677383/1‑126 (MQ=255)
                                                                                                                     aTGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTcc                        <  2:677383/126‑1 (MQ=255)
                                                                                                                     aTGCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATc          >  1:475264/1‑140 (MQ=255)
                                                                                                                       gCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGa        >  1:662383/1‑140 (MQ=255)
                                                                                                                       gCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGa        >  2:196951/1‑140 (MQ=255)
                                                                                                                       gCGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGa        >  2:193788/1‑140 (MQ=255)
                                                                                                                        cGGCACCAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGAt       <  1:193788/140‑1 (MQ=255)
                                                                                                                             ccAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAAt  >  1:1096410/1‑140 (MQ=255)
                                                                                                                             ccAAATGCCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAAt  >  1:1110884/1‑140 (MQ=255)
                                                                                                                                    |                                                                                                                                    
TGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAAT  >  NC_000913/2764417‑2764681

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