Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,076,151 T→A 5.2% N156I (AAT→ATT)  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,1510TA 5.2% 236.5 / 4.2 77N156I (AAT→ATT) PP_5328phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (35/38);  new base A (2/2);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

CGTTGACGCTGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGTCCGAAA  >  NC_002947/6076002‑6076298
                                                                                                                                                     |                                                                                                                                                   
cgttgacgCTGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCAt                                                                                                                                                     <  1:396795/150‑1 (MQ=255)
  ttgacgCTGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTg                                                                                                                                                   <  1:995415/150‑1 (MQ=255)
       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                                                                                                                                             >  1:285158/1‑107 (MQ=255)
                                                   aTTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGTCCAATGGCCAgg                                                                                                                                             <  2:285158/107‑1 (MQ=255)
                                                    ttGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGa                                                                                                 <  2:459946/150‑1 (MQ=255)
                                                     tGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGAc                                                                                                <  1:597059/150‑1 (MQ=255)
                                                         cGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTc                                                                                                             <  1:264445/133‑1 (MQ=255)
                                                         cGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTc                                                                                                             >  2:264445/1‑133 (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)
                                                                                                                                                     |                                                                                                                                                   
CGTTGACGCTGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGTCCGAAA  >  NC_002947/6076002‑6076298

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: