Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,076,155 A→C 6.1% S155A (TCC→GCC)  PP_5328 ← phosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,076,1550AC 6.1% 206.5 / 4.8 75S155A (TCC→GCC) PP_5328phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (33/37);  new base C (3/2);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

GCTGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGTCCGAAACAC  >  NC_002947/6076009‑6076301
                                                                                                                                                  |                                                                                                                                                  
gCTGACGTTGTCCAGCGCATGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCGTGTTGTCCGGGTAGGGGATCTTGTAGGTGTGCTGGAACCCCCGCGCCTGGCCATTGGACAg                                                                                                                                                 >  1:707237/1‑150 (MQ=255)
  tGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGTCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGc                                                                                                                                               <  2:538345/150‑1 (MQ=255)
  tGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGa                                                                                                                                                    >  2:1214769/1‑145 (MQ=255)
  tGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGa                                                                                                                                                    <  1:1214769/145‑1 (MQ=255)
  tGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGc                                                                                                                                               <  2:971809/150‑1 (MQ=255)
  tGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGc                                                                                                                                               <  2:947251/150‑1 (MQ=255)
      gTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGTCCAATGGCCAGGcgctg                                                                                                                                           <  1:274737/150‑5 (MQ=255)
      gTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAg                                                                                                                                           <  1:299765/150‑1 (MQ=255)
       tGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGc                                                                                                                                          <  1:709490/150‑1 (MQ=255)
           ccAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGcacca                                                                                                                                      >  2:613296/1‑150 (MQ=255)
            cAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGcaccac                                                                                                                                     <  2:691656/150‑1 (MQ=255)
            cAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGcaccac                                                                                                                                     >  2:65426/1‑150 (MQ=255)
            cAGCGCACGCCCCTGCTCGTCGAGCACCAACGCTTGATCGCCGTACGGGGAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGcaccac                                                                                                                                     <  1:408964/150‑1 (MQ=255)
              gcgcACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGCCAGCCCCAGCAccacca                                                                                                                                   >  2:418087/1‑150 (MQ=255)
                gcACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAgc                                                                                                                                 >  2:910249/1‑150 (MQ=255)
                  aCGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAgcgc                                                                                                                               >  2:1053042/1‑150 (MQ=255)
                      cccTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTg                                                                                                                           >  1:7765/1‑150 (MQ=255)
                         tGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGccc                                                                                                                        <  2:31155/150‑1 (MQ=255)
                             cgtcgAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCt                                                                                                                    <  1:967828/150‑1 (MQ=255)
                                cgAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGt                                                                                                                 >  2:597059/1‑150 (MQ=255)
                                  aGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCa                                                                                                               <  1:613296/150‑1 (MQ=255)
                                      cgaacgaTTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTAc                                                                                                           >  1:437826/1‑150 (MQ=255)
                                            aTTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGTCCAATGGCCAgg                                                                                                                                                <  2:285158/107‑1 (MQ=255)
                                            aTTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGTCCAATGGCCAgg                                                                                                                                                >  1:285158/1‑107 (MQ=255)
                                             ttGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGa                                                                                                    <  2:459946/150‑1 (MQ=255)
                                              tGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGAc                                                                                                   <  1:597059/150‑1 (MQ=255)
                                                  cGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTc                                                                                                                >  2:264445/1‑133 (MQ=255)
                                                  cGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTc                                                                                                                <  1:264445/133‑1 (MQ=255)
                                                   gCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGcacca                                                                                                                                      <  2:1166254/110‑1 (MQ=255)
                                                   gCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGcacca                                                                                                                                      >  1:1166254/1‑110 (MQ=255)
                                                   gCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGAc                                                                                              >  2:764348/1‑150 (MQ=255)
                                                   gCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGAc                                                                                              >  1:339832/1‑150 (MQ=255)
                                                            gTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTc                                                                                     >  1:366745/1‑150 (MQ=255)
                                                                  gATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCgg                                                                               >  2:123816/1‑150 (MQ=255)
                                                                  gATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCgg                                                                               <  1:65426/150‑1 (MQ=255)
                                                                       cAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGcacc                                                                                                                                       >  1:769187/1‑89 (MQ=255)
                                                                       cAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGcacc                                                                                                                                       <  2:769187/89‑1 (MQ=255)
                                                                        aGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACggg                                                                         <  1:764348/150‑1 (MQ=255)
                                                                        aGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACggg                                                                         >  1:143309/1‑150 (MQ=255)
                                                                             tGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGtt                                                                    >  2:147848/1‑150 (MQ=255)
                                                                             tGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGtt                                                                    <  1:1053042/150‑1 (MQ=255)
                                                                             cgatctttttCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGGCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGtt                                                                    <  1:296657/144‑1 (MQ=255)
                                                                              gATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTc                                                                   <  2:339832/150‑1 (MQ=255)
                                                                              gATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGACAGTACTGATGGACGTGACCTGGTTGTCTGCCGGCTCGGGGAGTTc                                                                   >  2:277950/1‑150 (MQ=255)
                                                                                tGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGa                                                                 <  2:7765/150‑1 (MQ=255)
                                                                                  gTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGaca                                                               <  2:225393/150‑1 (MQ=255)
                                                                                    tttCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGCTGTCTGCCCGCACGGGTAGTTcagatcg                                                             >  2:296657/1‑144 (MQ=255)
                                                                                       cTTGTTGTCCGGGTAGGTCATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCtt                                                          <  1:706136/150‑1 (MQ=255)
                                                                                             gtCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTg                                                                                   <  1:757623/119‑1 (MQ=255)
                                                                                             gtCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTg                                                                                   >  2:757623/1‑119 (MQ=255)
                                                                                                cGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTc                                                    >  1:1148372/1‑147 (MQ=255)
                                                                                                cGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTc                                                    <  2:1148372/147‑1 (MQ=255)
                                                                                                   gTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGgt                                                                                        <  2:802998/108‑1 (MQ=255)
                                                                                                   gTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGgt                                                                                        >  1:802998/1‑108 (MQ=255)
                                                                                                   gTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTAcc                                              <  2:707237/150‑1 (MQ=255)
                                                                                                    tAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCg                                             <  1:910249/150‑1 (MQ=255)
                                                                                                        tGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACgg                                                                          >  1:341807/1‑117 (MQ=255)
                                                                                                        tGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACgg                                                                          <  2:341807/117‑1 (MQ=255)
                                                                                                         gATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCtt                                        >  1:1280979/1‑150 (MQ=255)
                                                                                                                   gtgtGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCgaag                              >  2:1085040/1‑150 (MQ=255)
                                                                                                                          ggAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATg                       >  2:438144/1‑150 (MQ=255)
                                                                                                                            aaaaCCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGcc                     <  1:562700/150‑1 (MQ=255)
                                                                                                                             aaaCCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCt                                                                                    >  2:532995/1‑86 (MQ=255)
                                                                                                                             aaaCCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCt                                                                                    <  1:532995/86‑1 (MQ=255)
                                                                                                                             aaaCCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCt                    <  1:438144/150‑1 (MQ=255)
                                                                                                                              aaCCAGCGCCTGGCCATTGGAAAGGCCTAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGATCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTg                   <  2:223253/150‑1 (MQ=255)
                                                                                                                               tgctgggcCTGTCCAATGGCCAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTgg                  >  2:544318/7‑150 (MQ=255)
                                                                                                                               aCCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTgg                  <  2:143309/150‑1 (MQ=255)
                                                                                                                                 cAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAg                                                       <  1:98876/111‑1 (MQ=255)
                                                                                                                                 cAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAg                                                       >  2:98876/1‑111 (MQ=255)
                                                                                                                                   gcgcCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCt              >  1:1253655/1‑150 (MQ=255)
                                                                                                                                     gcCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTGCTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTg            >  1:778393/1‑150 (MQ=255)
                                                                                                                                     gcCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTg            >  1:1014495/1‑150 (MQ=255)
                                                                                                                                        tGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCATTACCGGGCTGGGTATCGAAGAACAATGCCTGGCCCTTGTcc         >  2:942433/1‑150 (MQ=255)
                                                                                                                                            cATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGTCCGaaa     >  1:283089/1‑150 (MQ=255)
                                                                                                                                            cATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGTCCGaaa     <  1:277950/150‑1 (MQ=255)
                                                                                                                                               tGGACAGGCGCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAACAACCATCGCGGGCCCTTGTCCGAAacac  <  2:437826/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCTGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGTCCGAAACAC  >  NC_002947/6076009‑6076301

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