Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,903,887 G→A 6.0% intergenic (+87/+20) PP_5176 → / ← spuH hypothetical protein/spermidine/putrescine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,903,8870GA 6.0% 337.2 / 7.1 117intergenic (+87/+20)PP_5176/spuHhypothetical protein/spermidine/putrescine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (52/56);  new base A (4/3);  total (58/59)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 117 total aligned reads displayed.

TGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTAC  >  NC_002947/5903746‑5904033
                                                                                                                                             |                                                                                                                                                  
   aTTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                                                                                         <  2:2401512/150‑1 (MQ=255)
    ttGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAg                                                                                                                                        >  1:1878550/1‑150 (MQ=255)
     tGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAgg                                                                                                                                       <  1:2170220/150‑1 (MQ=255)
       aacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccccttcgcgggcatGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAgg                                                                                                                                       <  1:422797‑M1/74‑1 (MQ=255)
        tCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCAGTGCGGTGCCTGCGGGATCTCGCGTGCCCGCTAAAGGGcc                                                                                                                                    >  2:2322403/1‑150 (MQ=255)
          aGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAg                                                                                                                                  <  1:2389261/150‑1 (MQ=255)
              gTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCa                                                                                                                              <  1:2262650/150‑1 (MQ=255)
               tATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAg                                                                                                                             >  2:1180687/1‑150 (MQ=255)
                 tCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGc                                                                                                                           >  2:1572777/1‑150 (MQ=255)
                   aGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCg                                                                                                                         >  1:1370469/1‑150 (MQ=255)
                           tCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGcc                                                                                                                 >  1:1784102/1‑150 (MQ=255)
                           tCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGCCCAGGGCAGGCCGACGCGGcc                                                                                                                 >  1:1793206/1‑150 (MQ=255)
                                            tATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGACAGGGTTTATCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCtt                                                                                                >  2:1872984/1‑150 (MQ=255)
                                              tGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGACTCACCACAACCCTTAAGGCTGTGAGGTCCCTGTAGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGCGCGCTTGa                                                                                              >  2:432136/1‑150 (MQ=255)
                                                 cTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTgcgc                                                                                                         <  2:857986/136‑1 (MQ=255)
                                                 cTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTgcgc                                                                                                         >  1:857986/1‑136 (MQ=255)
                                                 cTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCt                                                                                           <  1:1572777/150‑1 (MQ=255)
                                                  tttGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAATGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTg                                                                                          >  1:1098905/1‑150 (MQ=255)
                                                  tttGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGt                                                                                                                                >  1:1890001/1‑112 (MQ=255)
                                                  tttGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGt                                                                                                                                <  2:1890001/112‑1 (MQ=255)
                                                         cGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCACGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGc                                                                                   <  1:556909/150‑1 (MQ=255)
                                                            tGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGaa                                                                                <  2:233680/150‑1 (MQ=255)
                                                            tGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGaa                                                                                <  1:2322403/150‑1 (MQ=255)
                                                             gCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAg                                                                               >  2:698778/1‑150 (MQ=255)
                                                              cTGGAGCGTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGc                                                                              <  1:672220/150‑1 (MQ=255)
                                                              cTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGc                                                                              <  2:1966736/150‑1 (MQ=255)
                                                              cTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGc                                                                              <  1:1185090/150‑1 (MQ=255)
                                                                ggAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGt                                                                            >  1:508617/1‑150 (MQ=255)
                                                                  aGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTg                                                                          <  1:1180687/150‑1 (MQ=255)
                                                                   gCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTgg                                                                         <  2:653530/150‑1 (MQ=255)
                                                                   cctttcgcgggcacgcccgctcccacagggacctcacagccttcagggttgtggtgagccctgtgggagcgggcATGCCCGCGAAGGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGtt                                                                           >  2:422797‑M2/75‑148 (MQ=255)
                                                                   cctttcgcgggcacgcccgctcccacagggacctcacagccttcagggttgtggtgagccctgtgggagcgggcATGCCCGCGAAGGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGtt                                                                           <  1:422797‑M2/74‑1 (MQ=255)
                                                                    ccTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGc                                                                        <  2:2090984/150‑1 (MQ=255)
                                                                      ttCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTg                                                                      >  2:1438425/1‑150 (MQ=255)
                                                                        cgctGGCATGCCCGCTCCCACAGGGCTCACCGCAACCGTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTgca                                                                    <  2:508617/150‑1 (MQ=255)
                                                                             gCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCAcc                                                               <  2:347297/150‑1 (MQ=255)
                                                                               aTGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGACCCGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCct                                                             >  2:1540178/1‑150 (MQ=255)
                                                                                   CCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTc                                                         <  2:1340098‑M1/150‑1 (MQ=255)
                                                                                   CCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTc                                                         <  1:764544‑M1/150‑1 (MQ=255)
                                                                                   CCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTc                                                         <  1:52722‑M1/150‑1 (MQ=255)
                                                                                   CCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTc                                                         <  1:1441947‑M1/150‑1 (MQ=255)
                                                                                   CCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTc                                                         <  2:637708‑M1/150‑1 (MQ=255)
                                                                                   CCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTc                                                         <  1:1333263‑M1/150‑1 (MQ=255)
                                                                                     GCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGc                                                                                                                  <  1:606779‑M1/91‑1 (MQ=255)
                                                                                     GCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGc                                                                                                                  >  2:606779‑M1/1‑91 (MQ=255)
                                                                                     GCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGa                                                       <  2:1909454‑M1/150‑1 (MQ=255)
                                                                                       tcccacagggacctcacagccttcagggttgtggtgagccctgtgggagcgggcGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCtt                                                                                                <  1:67857‑M2/53‑1 (MQ=255)
                                                                                       tcccacagggacctcacagccttcagggttgtggtgagccctgtgggagcgggcGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCtt                                                                                                >  2:67857‑M2/55‑107 (MQ=255)
                                                                                       tcccacagggacctcacagccttcagggttgtggtgagccctgtgggagcgggcATGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATc                                                     >  2:190423‑M2/55‑150 (MQ=255)
                                                                                         ccacagggacctcacagccttcagggttgtggtgagccctgtgggagcgggcTGAACCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGc                                                   >  1:1978224‑M2/53‑150 (MQ=255)
                                                                                         ccacagggacctcacagccttcagggttgtggtgagccctgtgggagcgggcATGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCtt                                                                                                >  2:1962401‑M2/53‑105 (MQ=255)
                                                                                         ccacagggacctcacagccttcagggttgtggtgagccctgtgggagcgggcATGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCtt                                                                                                <  1:1962401‑M2/53‑1 (MQ=255)
                                                                                         CCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGc                                                   <  1:945749‑M1/150‑1 (MQ=255)
                                                                                          cacagggacctcacagccttcagggttgtggtgagccctgtgggagcgggcGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCgg                                                                                                      <  2:884385‑M2/47‑1 (MQ=255)
                                                                                          cacagggacctcacagccttcagggttgtggtgagccctgtgggagcgggcGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGcc                                                  >  1:668892‑M2/52‑150 (MQ=255)
                                                                                          CACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGc                                                                                                                           <  2:2020390‑M1/77‑1 (MQ=255)
                                                                                          CACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGc                                                                                                                           >  1:2020390‑M1/1‑77 (MQ=255)
                                                                                          CACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGc                                                                                         <  2:1279345‑M1/111‑1 (MQ=255)
                                                                                          CACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCAACCTGCGCCGGGCGCTTGACCTGc                                                                                         >  1:1279345‑M1/1‑111 (MQ=255)
                                                                                           acagggacctcacagccttcagggttgtggtgagccctgtgggagcgggcATGCCCGCGAAAGGGCCAGGTCAGGCCGAAGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCct                                                 >  1:1764258‑M2/51‑150 (MQ=255)
                                                                                           ACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGc                                                                              >  1:1361759‑M1/1‑121 (MQ=255)
                                                                                           ACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGc                                                                              <  2:1361759‑M1/121‑1 (MQ=255)
                                                                                             aGGGCTCACCACAACGCTGAAGGTTGTGAGGTGCCTGTGGGGGAGGGCGTGACAGCGAGAGGGCCAGGCCAGGCCGATGCAGCCATCCTCGCCGAGCGCTTGTCCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTGTTCGATCGCctgc                                               <  2:2079291/150‑1 (MQ=255)
                                                                                                gacctcacagccttcagggttgtggtgagccctgtgggagcgggcGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTg                                                                          >  1:2089068‑M2/46‑120 (MQ=255)
                                                                                                gacctcacagccttcagggttgtggtgagccctgtgggagcgggcGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTg                                                                          <  2:2089068‑M2/75‑1 (MQ=255)
                                                                                                  ggacctcacagccttcagggttgtggtccctgtgggagcgggcGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAAt                                          >  1:591560‑M2/44‑150 (MQ=255)
                                                                                                  acctcacagccttcagggttgtggggtccctgtgggagcgggcGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTCGGc                                                                        <  1:2194467‑M2/77‑1 (MQ=255)
                                                                                                  acctcacagccttcagggttgtggggtccctgtgggagcgggcGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTCGGc                                                                        >  2:2194467‑M2/44‑120 (MQ=255)
                                                                                                     ccacAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGc                                       <  1:21179/150‑1 (MQ=255)
                                                                                                      cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGcc                                      >  2:2348044/1‑150 (MQ=255)
                                                                                                       acaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCt                                     >  1:916218/1‑150 (MQ=255)
                                                                                                         aaCCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTg                                   >  1:197239/1‑150 (MQ=255)
                                                                                                          ccccTGAAGGCTGTAAGGGCCGTGTGGGAGGGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGc                                  <  1:378541/149‑1 (MQ=255)
                                                                                                          cagccttcagggttgtggtccctgtgggagcgggcGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGc                                  <  2:591560‑M2/115‑1 (MQ=255)
                                                                                                           cccTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcg                                 >  2:353858/1‑150 (MQ=255)
                                                                                                            ccTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcga                                <  2:495288/150‑1 (MQ=255)
                                                                                                             cTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGc                                       <  2:960188/142‑1 (MQ=255)
                                                                                                             cTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGc                                       >  1:960188/1‑142 (MQ=255)
                                                                                                              tGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGc                                                                                                                  >  2:1750452/1‑66 (MQ=255)
                                                                                                              tGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGtt                                                                           <  2:1934035/105‑1 (MQ=255)
                                                                                                              tGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGtt                                                                           >  1:1934035/1‑105 (MQ=255)
                                                                                                              tGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcgacg                              >  2:921648/1‑150 (MQ=255)
                                                                                                              tGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcgacg                              >  2:1211971/1‑150 (MQ=255)
                                                                                                              tGAAGGCTGTGACGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGc                                                                                                                  <  1:1750452/66‑1 (MQ=255)
                                                                                                               taaGGCTGTGAGGTCCCTGTGGGAGCGGCCGTGCACGAGAAAGGGCCAGGTCCGCCCGATGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcgacga                             <  1:398214/149‑1 (MQ=255)
                                                                                                                  ggCTGTGAGGTCCCTGTGGGAACGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGc                          >  2:498885/1‑150 (MQ=255)
                                                                                                                     tgtgAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGGCGctct                       >  1:1307700/1‑150 (MQ=255)
                                                                                                                           gTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGcc                 <  1:2255779/150‑1 (MQ=255)
                                                                                                                            tCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGc                                  <  2:512064/132‑1 (MQ=255)
                                                                                                                            tCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGc                                  >  1:512064/1‑132 (MQ=255)
                                                                                                                              cctgtgggagcgggcGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTg                                                                                          <  1:1260066‑M2/59‑1 (MQ=255)
                                                                                                                              cctgtgggagcgggcGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTg                                                                                          >  2:1260066‑M2/16‑74 (MQ=255)
                                                                                                                                  tgggagcgggcGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTggcggc          >  1:1528967‑M2/12‑150 (MQ=255)
                                                                                                                                    ggagcgggcGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAATCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTg        >  1:2098475‑M2/10‑150 (MQ=255)
                                                                                                                                      agcgggcGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGaa      <  2:916218‑M2/143‑1 (MQ=255)
                                                                                                                                      agcgggcATGCCCGCGAAAGGGCCAGGTCAGGCCGAAGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGaa      <  2:1764258‑M2/143‑1 (MQ=255)
                                                                                                                                       gcgggcGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAg     <  2:1333263‑M2/144‑1 (MQ=255)
                                                                                                                                          ggcGTGCCCGCGAAATGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTAc  <  2:1098905‑M2/147‑1 (MQ=255)
                                                                                                                                          ggcGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTAc  >  1:747936‑M2/4‑150 (MQ=255)
                                                                                                                                          ggcGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTAc  <  2:1607641‑M2/147‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
TGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTAC  >  NC_002947/5903746‑5904033

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