Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,270,641 C→A 5.1% S419S (TCC→TCA PP_4643 → xanthine/uracil permease family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,270,6410CA 5.1% 725.3 / 19.4 216S419S (TCC→TCAPP_4643xanthine/uracil permease family protein
Reads supporting (aligned to +/- strand):  ref base C (106/99);  new base A (5/6);  total (111/105)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 218 total aligned reads displayed.

AGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGG  >  NC_002947/5270492‑5270788
                                                                                                                                                     |                                                                                                                                                   
aGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTcc                                                                                                                                                     >  2:1046838/1‑150 (MQ=255)
       tCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCa                                                                                                                                              >  2:1924309/1‑150 (MQ=255)
         cGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGc                                                                                                                                            <  2:4361661/150‑1 (MQ=255)
          gCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCt                                                                                                                                           >  1:4619761/1‑150 (MQ=255)
          gCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCt                                                                                                                                           >  1:2482347/1‑150 (MQ=255)
                 gTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCggcggc                                                                                                                                    <  2:542957/150‑1 (MQ=255)
                      aGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCa                                                                                                                               >  2:4044868/1‑150 (MQ=255)
                      aGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCGCGGCGGCGATCa                                                                                                                               >  2:5392760/1‑150 (MQ=255)
                       ggTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCAt                                                                                                                              >  1:5212814/1‑150 (MQ=255)
                        gTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATg                                                                                                                             <  2:3347133/150‑1 (MQ=255)
                        gTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATg                                                                                                                             <  1:1309880/150‑1 (MQ=255)
                          cAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGc                                                                                                                           >  1:3585824/1‑150 (MQ=255)
                              tatcaGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGCTCCGCTGCATTGATGGCAATc                                                                                                                       >  1:2515138/5‑150 (MQ=255)
                              tACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATc                                                                                                                       <  1:1183456/150‑1 (MQ=255)
                               aCAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATcc                                                                                                                      >  1:596595/1‑150 (MQ=255)
                                cAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCct                                                                                                                     >  2:1459225/1‑150 (MQ=255)
                                 aaGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCctg                                                                                                                    >  1:2914486/1‑150 (MQ=255)
                                  aGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCctgc                                                                                                                   >  2:1640420/1‑150 (MQ=255)
                                   gAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCctgct                                                                                                                  >  2:2425053/1‑150 (MQ=255)
                                   gAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGTTTGAAACCATTTTTCATTCAGGCATCAGCTCCGCTGCATTGATGGCAATcgcgct                                                                                                                  >  2:286053/1‑145 (MQ=255)
                                       aacGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAAc                                                                                                              <  2:926912/150‑1 (MQ=255)
                                       aacGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAAc                                                                                                              <  1:29862/150‑1 (MQ=255)
                                        ccgccAACCTGCACATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGCCCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAAcc                                                                                                             <  2:215383/146‑1 (MQ=255)
                                         cGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCt                                                                                                            >  1:1909698/1‑150 (MQ=255)
                                          gTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTg                                                                                                           <  2:4754118/150‑1 (MQ=255)
                                             aaCCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATc                                                                                                        >  2:746435/1‑150 (MQ=255)
                                               ccTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTcgg                                                                                                                                        >  2:5083442/1‑116 (MQ=255)
                                               ccTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTcgg                                                                                                                                        <  1:5083268/116‑1 (MQ=255)
                                               ccTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCtt                                                                                                      <  1:3895145/150‑1 (MQ=255)
                                                cTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTc                                                                                                     <  2:1434578/150‑1 (MQ=255)
                                                cTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTc                                                                                                     <  1:1412893/150‑1 (MQ=255)
                                                 tGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCa                                                                                                    <  1:1589102/150‑1 (MQ=255)
                                                 tGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCa                                                                                                    <  2:3698887/150‑1 (MQ=255)
                                                  gATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCaa                                                                                                   <  2:1977743/150‑1 (MQ=255)
                                                    tcatcGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAAcc                                                                                                 >  2:172495/1‑150 (MQ=255)
                                                     catcGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCa                                                                                                >  1:1275587/1‑150 (MQ=255)
                                                     catcGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCa                                                                                                >  1:3492857/1‑150 (MQ=255)
                                                      atcGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCAt                                                                                               <  2:1698108/150‑1 (MQ=255)
                                                          tcgCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTca                                                                                           >  2:42622/1‑150 (MQ=255)
                                                                    gCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCaa                                                                                 <  2:4782939/150‑1 (MQ=255)
                                                                      tGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACt                                                                               <  1:4989341/150‑1 (MQ=255)
                                                                       gggCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTc                                                                              <  2:1497730/150‑1 (MQ=255)
                                                                          cTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGa                                                                           <  1:4453578/150‑1 (MQ=255)
                                                                           ttCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGaa                                                                          <  2:3510259/150‑1 (MQ=255)
                                                                            tCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGaaa                                                                         <  1:3452161/150‑1 (MQ=255)
                                                                            tCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGaaa                                                                         <  1:2538218/150‑1 (MQ=255)
                                                                              ggCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAAc                                                                       <  1:2028317/150‑1 (MQ=255)
                                                                              ggCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAAc                                                                       <  1:866961/150‑1 (MQ=255)
                                                                               gCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAAcc                                                                      >  2:4477191/1‑150 (MQ=255)
                                                                                cATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACca                                                                     >  1:5432447/1‑150 (MQ=255)
                                                                                     tCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACcagtcag                                                                >  2:1204390/1‑150 (MQ=255)
                                                                                     tCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACcagtcag                                                                >  1:4747879/1‑150 (MQ=255)
                                                                                         gATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGtt                                                            <  2:3059547/150‑1 (MQ=255)
                                                                                         gATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGtt                                                            <  1:4670280/150‑1 (MQ=255)
                                                                                          aTTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTc                                                           <  1:746435/150‑1 (MQ=255)
                                                                                              ccGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTg                                                       >  2:2227231/1‑150 (MQ=255)
                                                                                              ccGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTg                                                       >  2:301560/1‑150 (MQ=255)
                                                                                              ccGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTg                                                       >  1:5406238/1‑150 (MQ=255)
                                                                                              ccGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTg                                                       >  1:4119026/1‑150 (MQ=255)
                                                                                               cGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTcggcg                                                                                                                                      <  1:1070934/70‑1 (MQ=255)
                                                                                               cGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTcggcg                                                                                                                                      >  2:1070934/1‑70 (MQ=255)
                                                                                                gCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCa                                                     >  1:2112007/1‑150 (MQ=255)
                                                                                                   ccGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCACCTCGGAAAACCAGTCAGTGTTCGCTGCAGcc                                                  >  2:3388799/1‑150 (MQ=255)
                                                                                                     gACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTa                                                <  1:3076118/150‑1 (MQ=255)
                                                                                                      aCCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATTGACATCATGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCTGAAAACCAGTCAGTGTTCGCTGCAGCCTAc                                               <  2:3603723/150‑1 (MQ=255)
                                                                                                      aCCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTAc                                               >  2:4591370/1‑150 (MQ=255)
                                                                                                      aCCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTAc                                               >  2:3106002/1‑150 (MQ=255)
                                                                                                       ccTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACg                                              >  1:3189443/1‑150 (MQ=255)
                                                                                                        cTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGa                                             >  1:4235978/1‑150 (MQ=255)
                                                                                                         ttCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACAAGTCAGTGTTCGCTGCAGCCTACGAg                                            >  1:723669/1‑150 (MQ=255)
                                                                                                          tCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAgc                                           >  1:3604314/1‑150 (MQ=255)
                                                                                                            tACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAgcgc                                         >  1:4388305/1‑150 (MQ=255)
                                                                                                             accaccACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCa                                        >  2:2021909/1‑150 (MQ=255)
                                                                                                             accaccACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCa                                        <  2:1353993/150‑1 (MQ=255)
                                                                                                               caccacTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAgcg                                          <  2:1196978/146‑1 (MQ=255)
                                                                                                               caccacTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCAcc                                      >  2:3248719/1‑150 (MQ=255)
                                                                                                               caccacTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCAcc                                      >  1:229216/1‑150 (MQ=255)
                                                                                                               caccacTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCAcc                                      >  2:298445/1‑150 (MQ=255)
                                                                                                                accacTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCa                                     <  2:614806/150‑1 (MQ=255)
                                                                                                                accacTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCa                                     >  2:1645751/1‑150 (MQ=255)
                                                                                                                 ccacTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCAt                                    >  2:4109675/1‑150 (MQ=255)
                                                                                                                 ccacTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCAt                                    >  1:621545/1‑150 (MQ=255)
                                                                                                                   acTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTAc                                               >  1:3515599/1‑137 (MQ=255)
                                                                                                                   acTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATcc                                  >  2:5027221/1‑150 (MQ=255)
                                                                                                                    cTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCa                                 <  2:633028/150‑1 (MQ=255)
                                                                                                                    cTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCa                                 <  2:3585824/150‑1 (MQ=255)
                                                                                                                    cTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCa                                 <  1:1764466/150‑1 (MQ=255)
                                                                                                                     gtCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAg                                <  2:2494668/149‑1 (MQ=255)
                                                                                                                          gAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTc                           >  1:2948132/1‑150 (MQ=255)
                                                                                                                           aaCTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCg                          <  2:2814471/150‑1 (MQ=255)
                                                                                                                           aaCTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCg                          >  1:3818310/1‑150 (MQ=255)
                                                                                                                              tGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGAc                       >  1:1095334/1‑150 (MQ=255)
                                                                                                                              tGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGAc                       >  1:2551948/1‑150 (MQ=255)
                                                                                                                                    gagaCCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCGGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTcc                 >  2:5168678/1‑150 (MQ=255)
                                                                                                                                           tCTTCCACTCCGGCATCAGCTCGGCGTCGATAATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCCCCACCGGCAACTCGGAAAACCTGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGc          >  1:5350413/1‑150 (MQ=255)
                                                                                                                                             ttCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGt        >  2:3703386/1‑150 (MQ=255)
                                                                                                                                                  cTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTc                                                                                            <  1:321254/61‑1 (MQ=255)
                                                                                                                                                  cTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTc                                                                                            >  2:321254/1‑61 (MQ=255)
                                                                                                                                                  cTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATgg   <  1:2801943/150‑1 (MQ=255)
                                                                                                                                                   tCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATggg  <  2:3189443/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGG  >  NC_002947/5270492‑5270788

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