Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,952,947 T→G 6.9% Q32P (CAG→CCG)  PP_3483 ← type II secretion system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,952,9470TG 6.9% 653.7 / 3.4 239Q32P (CAG→CCG) PP_3483type II secretion system protein
Reads supporting (aligned to +/- strand):  ref base T (117/97);  new base G (13/12);  total (130/109)
Fisher's exact test for biased strand distribution p-value = 8.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 248 total aligned reads displayed.

GTGTCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGG  >  NC_002947/3952800‑3953096
                                                                                                                                                   |                                                                                                                                                     
    cGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTGTGTGGGc                                                                                                                                                 >  1:4029531/1‑150 (MQ=255)
    cGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGc                                                                                                                                                 <  2:332966/150‑1 (MQ=255)
     gCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACCGGTACATCACAGGTCTGTGGGCt                                                                                                                                                >  1:3275842/1‑150 (MQ=255)
     gCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCt                                                                                                                                                <  1:945827/150‑1 (MQ=255)
     gCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGGACATCACAGGTCGGTGGGCt                                                                                                                                                >  1:2109823/1‑150 (MQ=255)
         gccaggcGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCCTCACAGGTCTGTGGGCTGCAt                                                                                                                                            >  2:4026325/1‑150 (MQ=255)
         gccaggcGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCt                                                                                                                                                       <  2:5777810/139‑1 (MQ=255)
           caggcGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACGGGTACATCACAGGTCTGTGGGCTGCATAc                                                                                                                                          >  2:969074/1‑150 (MQ=255)
           caggcGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACCGGTACAGCACAGGTCTGTGGGCTGCATAc                                                                                                                                          >  2:424676/1‑150 (MQ=255)
            gccgagccagccggcgacgcgcccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACa                                                                                                                                         >  1:5851861‑M1/95‑150 (MQ=255)
            aggcGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACCGGTACCTCACAGGTCTGGGGGCTGCATACa                                                                                                                                         >  2:2347726/1‑150 (MQ=255)
              gCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGGAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATc                                                                                                                                       >  1:5087159‑M2/1‑149 (MQ=255)
                gTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcct                                                                                                                                     >  1:2606163‑M2/1‑147 (MQ=255)
                 ttCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctg                                                                                                                                    >  1:5199442‑M2/1‑146 (MQ=255)
                    cAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgg                                                                                                                                 >  2:6404708‑M2/1‑143 (MQ=255)
                    cAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgg                                                                                                                                 >  1:3511136‑M2/1‑143 (MQ=255)
                    cAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgg                                                                                                                                 >  2:606430‑M2/1‑143 (MQ=255)
                    agcctgcgacgcgcccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTgg                                                                                                                                 >  2:4193606‑M1/87‑150 (MQ=255)
                     aGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCATCACAGGTCTGTGGGCTGCATACATcctgtggg                                                                                                                                <  2:3928668‑M2/150‑9 (MQ=255)
                     aGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtggg                                                                                                                                <  2:2275564‑M2/150‑9 (MQ=255)
                      ggCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtggga                                                                                                                               >  1:2510165‑M2/1‑141 (MQ=255)
                      cctgcgacgcgcccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGa                                                                                                                               >  2:1705072‑M1/85‑150 (MQ=255)
                       gCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagc                                                                                                                             >  1:691552‑M2/1‑139 (MQ=255)
                         aGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcg                                                                                                                            <  2:732369‑M2/150‑13 (MQ=255)
                         aGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCATCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcg                                                                                                                            <  1:5009334‑M2/150‑13 (MQ=255)
                          gTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgg                                                                                                                           >  1:134605‑M2/1‑137 (MQ=255)
                          cgacgcgcccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCgg                                                                                                                           >  1:4339007‑M1/81‑150 (MQ=255)
                          cgacgcgcccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCgg                                                                                                                           >  1:6752954‑M1/81‑150 (MQ=255)
                          cgacgcgcccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCgg                                                                                                                           >  1:666279‑M1/81‑150 (MQ=255)
                           tCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACCTCACAGGTCTGTGGGCTGCATACATcctgggggagcggg                                                                                                                          <  2:3371611‑M2/150‑15 (MQ=255)
                           tCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcggg                                                                                                                          >  2:1976718‑M2/1‑136 (MQ=255)
                           tCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcggg                                                                                                                          >  1:2027436‑M2/1‑136 (MQ=255)
                             gAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACCGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgggca                                                                                                                        <  2:4736418‑M2/150‑17 (MQ=255)
                             gAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACCTCACAGGTCGGGGGGCTGCATACATcgggtgggagcgggcg                                                                                                                        <  1:5525255‑M2/150‑17 (MQ=255)
                              aGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgggcat                                                                                                                       <  2:3608727‑M2/150‑18 (MQ=255)
                              aGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgggcat                                                                                                                       >  1:5737424‑M2/1‑133 (MQ=255)
                               gCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgggcatg                                                                                                                      >  1:1166128‑M2/1‑132 (MQ=255)
                               gCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACAAcctgtgggagcggggctg                                                                                                                      <  1:6508891‑M2/150‑19 (MQ=255)
                                gcccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCAGGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGc                                                                                                                     >  2:512597‑M1/75‑150 (MQ=255)
                                cAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgggcatgc                                                                                                                     >  2:5203591‑M2/1‑131 (MQ=255)
                                cAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcggg                                                                                                                          <  2:190170‑M2/145‑15 (MQ=255)
                                  ggCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgggcatgccc                                                                                                                   >  2:5870505‑M2/1‑129 (MQ=255)
                                   gcgcTCGTCAGTGGCTTGCAGACGGGCCTGGGCCCGGAGGTCGCGGGGGTAGGCGGAAATTGTTGCTGGGTTCTTCGCGGGCATGCCCGCGCCCACAGGGAGCTCACGGGTCTGTGGGCTGCCTACAGccggggggagcgggccgggccg                                                                                                                  <  1:784020‑M2/150‑23 (MQ=255)
                                     gcTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCGGTGGGCTGCATACATcctgtgggagcgggcctgcccccg                                                                                                                <  1:6435694‑M2/150‑25 (MQ=255)
                                      cTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATGTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCGGTGGGCTGCATACATccggtgggagcgggcatgcccgcga                                                                                                               <  2:6507181‑M2/150‑26 (MQ=255)
                                      atcGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTACCACAGGTACATCACAGGTCTGTGGGCTGCATACATccggtgggagcgggcatgcccccga                                                                                                               <  2:5628315‑M2/149‑26 (MQ=255)
                                       tCGTCTGTGGCTTGCAGTGGGGCCTGGGCCTGGAGGTCGGGGGGTTTGGCGTAAATTTTTGCTGGTTTGTTCGCGGCCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgggcatgcccgcgaa                                                                                                              >  2:4342866‑M2/1‑124 (MQ=255)
                                       tCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACCTcctgtgggagcgggcatgcccgcgac                                                                                                              <  2:1179620‑M2/150‑27 (MQ=255)
                                       gggtcTGTGGGTGGGAGACGGGCGTGGGCCCGGCGGTCGCGGGGGTCGGCGGAAATTTTGGCTGGGTTCTTCGCGGGCTTGCCCGCTCCCACTGTTACATCACTGGTCTGTGGGCTGCATACATcctgtgggagcgggcttgcccgcgaa                                                                                                              >  2:3008607‑M2/3‑124 (MQ=255)
                                           cAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCGGTGGGCGGCATAAATcctgtgggagcgggcctgcccgcgccgaag                                                                                                          <  1:101051‑M2/150‑31 (MQ=255)
                                                  ttGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGACGACGTCGGTGc                                                                                                   >  2:5763407/1‑150 (MQ=255)
                                                  ttGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCGCCCACCGGTCCCTCCCAGGGCGGGGGGCGGCCTACATCCGGTGGGTGCGGGGCTTCCCGCGCGGACGTCGGGGc                                                                                                   >  1:3068126/1‑150 (MQ=255)
                                                   tGCAGCCGGGCCTGGGCCTGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTGCGCGGGCATGCCCGCTCCCACAGGTACATCACCGGTCTGGGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCa                                                                                                  <  2:1491480/150‑1 (MQ=39)
                                                   tGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCGGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCa                                                                                                  >  2:6785449/1‑150 (MQ=255)
                                                    gCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCCCAGGGACATCACCGGTCTGGGGGCGGGTTTCATCCGGTGGGTGCGGGGGTGGGCGCGAAGACGGCGGTGCAg                                                                                                 >  1:5326320/1‑150 (MQ=255)
                                                    gCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGGTGCATAGATCCGGTGGGAGCGGGGATGACCGCGCCGACGTCGGTGCAg                                                                                                 >  1:4293360/1‑150 (MQ=39)
                                                        aCGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGGGCCGGcgt                                                                                             >  2:4292298/1‑148 (MQ=255)
                                                                    cAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtg                                                                                                                                  >  2:4316101‑M2/1‑95 (MQ=255)
                                                                    cAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCgg                                                                                 >  2:5566982/1‑150 (MQ=255)
                                                                        tggTCGCGGGGGTTGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCATCACCGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGtt                                                                             <  2:1057387/149‑1 (MQ=255)
                                                                         ggTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTg                                                                            >  2:1250029/1‑150 (MQ=255)
                                                                          gTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGGGGGAGCGGGCATGCCCGCGACGACGTCGGTGCCGGCTGTGCAGATTCCGGGGTTGc                                                                           >  2:5308918/1‑150 (MQ=255)
                                                                            cgggggggTCGGCTTAATTTTTGGCGGGGTTCGTCGCGGGCCGGCCCGCCCCCCCAGGTACATCACCGTTCGGTGGGCTGAATCCATCCTGTGGCCGCGGTCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCAt                                                                         <  2:955535/147‑1 (MQ=255)
                                                                                 gggTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCTGAGGTCTGTGGGCTTCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAAt                                                                    <  2:1553499/150‑1 (MQ=255)
                                                                                    tAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATAACAGGTCTGTGGGCTGCATACATCCTGAGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGAGTCCGGGGTTGCATCGAATGCt                                                                 >  2:2750516/1‑150 (MQ=255)
                                                                                                tttGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGGGCATCACCGGGCTGTGGGCTGCATACATCCTGTGGGAGCGGGCAGGCCCCCGACGACGGCGGTGCAGGCTTTGCAGATTCCGGGGGTGCAACGAAGGCTTCTGGCAGGGcc                                                     >  1:2120138/1‑150 (MQ=255)
                                                                                                   gcggggTTCTTCGGGGGCTGGCCCGCTCCCACCGGTACCTCCCCGGTCTGTGGGCGGCATCCTGCGGGTGGGCGCGGGCTTGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGc                                                  <  2:1778258/147‑1 (MQ=255)
                                                                                                   gCTGGGTTCTTCGCGGGCCTGCCCCCTCCCACAGGTACACCCCAGGTCTGTGGGCTGCATACATCCGGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGc                                                  <  2:6326197/150‑1 (MQ=255)
                                                                                                   gCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCGGTGGGCTGCATACATCCGGGGGGAGCGGGCAAGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGc                                                  >  2:2171072/1‑150 (MQ=255)
                                                                                                     tGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCCTCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCCTCGAATGCTTCTGGCCTGGCCTGCGa                                                >  2:4697930/1‑150 (MQ=255)
                                                                                                            TTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGCTTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCt                                         <  2:4938034‑M1/150‑1 (MQ=255)
                                                                                                                        cgCGCTCCCACAGGTACATCCCTGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTc                             >  2:6246088‑M1/3‑150 (MQ=255)
                                                                                                                        CCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTc                             <  1:2896601‑M1/150‑1 (MQ=255)
                                                                                                                         CCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCg                            <  1:1774247‑M1/150‑1 (MQ=255)
                                                                                                                         CCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCg                            <  1:3877436‑M1/150‑1 (MQ=255)
                                                                                                                         CCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCg                            <  1:1901459‑M1/150‑1 (MQ=255)
                                                                                                                         CCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCg                            <  2:4447408‑M1/150‑1 (MQ=255)
                                                                                                                         CCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCg                            <  1:6561098‑M1/150‑1 (MQ=255)
                                                                                                                          cGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagccactgacgagcgcctgctcgact                           <  2:146777‑M2/150‑110 (MQ=255)
                                                                                                                            CTCCCACCGGTACCTCACCGGTCGGTGGGCTGCATACATCGGGGGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAg                         >  2:2457871‑M1/1‑150 (MQ=255)
                                                                                                                               ccACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgggcatgcccgcgaagaatccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagccactgacgagcgcctgctcgactgcctg                      <  1:5048621‑M2/150‑115 (MQ=255)
                                                                                                                               ccACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagccactgacgagcgcctgctcgactgcctg                      <  1:1667195‑M2/150‑115 (MQ=255)
                                                                                                                               ccACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctgg                                                                     >  1:468643‑M2/1‑36 (MQ=255)
                                                                                                                               ccACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctgg                                                                     <  2:3888470‑M2/103‑68 (MQ=255)
                                                                                                                               ccACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctgg                                                                     >  1:3888470‑M2/1‑36 (MQ=255)
                                                                                                                               CCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAggct                      <  2:4572389‑M1/150‑1 (MQ=255)
                                                                                                                                         aTCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCa            >  2:3671776/1‑150 (MQ=255)
                                                                                                                                          tCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAg           >  1:220285/1‑150 (MQ=255)
                                                                                                                                          tCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAg           >  1:4186547/1‑150 (MQ=255)
                                                                                                                                          tCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAg           >  1:5630043/1‑150 (MQ=255)
                                                                                                                                           cacaGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAgg          >  1:4069169/1‑150 (MQ=255)
                                                                                                                                            acaGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAggg         >  1:3205982/1‑150 (MQ=255)
                                                                                                                                             caGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAgggg        >  1:5993529/1‑150 (MQ=255)
                                                                                                                                              aGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGt       >  1:5012442/1‑150 (MQ=255)
                                                                                                                                              aGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGt       >  1:6290861/1‑150 (MQ=255)
                                                                                                                                                  cTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAggg         <  1:6216738/144‑1 (MQ=255)
                                                                                                                                                  cTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCgg   >  1:5222983/1‑150 (MQ=255)
                                                                                                                                                   tgtgGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCggg  >  2:4616468/1‑150 (MQ=255)
                                                                                                                                                   tgtgGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCggg  >  1:493604/1‑150 (MQ=255)
                                                                                                                                                   tgtgGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCggg  >  1:6344018/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTGTCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGG  >  NC_002947/3952800‑3953096

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