Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,741,564 A→T 8.7% intergenic (‑179/+49) PP_5038 ← / ← PP_5039 lipoprotein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,741,5640AT 8.7% 363.8 / 21.3 126intergenic (‑179/+49)PP_5038/PP_5039lipoprotein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (64/51);  new base T (7/4);  total (71/55)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

Only 100 of 130 total aligned reads displayed.

CTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGAAGCGGGCAA  >  NC_002947/5741423‑5741713
                                                                                                                                             |                                                                                                                                                     
cTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCg                                                                                                                                               <  2:919421/150‑1 (MQ=255)
  gCAGTGACATTCGACGTGAAAACCACGCAGAGGTGCCCGTCCCCACCCCCTATGTCCCTACACCATTGGTGGAATGGCCACTACTGTGTGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGc                                                                                                                                             <  1:1001049/150‑1 (MQ=255)
   cAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGcc                                                                                                                                            >  2:719490/1‑150 (MQ=255)
      tGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCTCCCATTGGGGGCGGCCttt                                                                                                                                         >  1:237386/1‑149 (MQ=255)
       gACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCgg                                                                                                                                              >  1:821087/1‑144 (MQ=255)
       gACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCgg                                                                                                                                              <  2:821087/144‑1 (MQ=255)
       gACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGc                                                                                                                                        <  1:195453/150‑1 (MQ=255)
        aCCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCg                                                                                                                                       >  2:1469648/1‑150 (MQ=255)
         ccTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCgg                                                                                                                                      >  2:1790481/1‑149 (MQ=255)
          cTTCGACGTGAAACTCACTCCCAGGTTCCCTGCCACACCCCCTATGTCCCTACCCCATTGGTCGTATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTc                                                                                                                                     <  1:990486/150‑1 (MQ=255)
          cTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTc                                                                                                                                     >  1:559059/1‑150 (MQ=255)
          cTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCATAACCATTTGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATTCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTc                                                                                                                                     >  2:332545/1‑150 (MQ=255)
           ttCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTcg                                                                                                                                    >  1:1141191/1‑150 (MQ=255)
           ttCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTcg                                                                                                                                    >  1:2303351/1‑150 (MQ=255)
            tCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCTCCCATTGGGGGCGGCCttttttatt                                                                                                                                   <  2:237386/150‑8 (MQ=255)
            tCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTcgc                                                                                                                                   <  2:1967659/150‑1 (MQ=255)
               aCGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGAt                                                                                                                                >  2:1869933/1‑150 (MQ=255)
               aCGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGAt                                                                                                                                >  2:1869893/1‑150 (MQ=255)
                cGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATg                                                                                                                               >  1:960070/1‑150 (MQ=255)
                 gTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGc                                                                                                                              >  1:498200/1‑150 (MQ=255)
                      aaTCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCTCCCATTggggg                                                                                                                                                 >  2:210218/1‑124 (MQ=255)
                      aaTCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCTCCCATTggggg                                                                                                                                                 <  1:210218/126‑3 (MQ=255)
                         ccctcccAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCAcag                                                                                                                      <  1:905079/150‑1 (MQ=255)
                           ctcccAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCAcagca                                                                                                                    >  2:438258/1‑150 (MQ=255)
                            tcccAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCAcagcag                                                                                                                   <  1:1869893/150‑1 (MQ=255)
                            tcccAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCAcagcag                                                                                                                   <  1:1869933/150‑1 (MQ=255)
                            tcccAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCAcagcag                                                                                                                   >  1:613825/1‑150 (MQ=255)
                            tcccAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCAcagcag                                                                                                                   >  2:562764/1‑150 (MQ=255)
                             cccAGGTTCCCTGCCCCCCCCCCTATGTCCCTTCACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAgg                                                                                                                  <  1:1613231/150‑1 (MQ=255)
                              ccAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAAt                                                                                                                                                      <  1:1807414/113‑1 (MQ=255)
                              ccAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAAt                                                                                                                                                      >  2:1807414/1‑113 (MQ=255)
                              ccAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGc                                                                                                                 <  2:519740/150‑1 (MQ=255)
                              ccAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGc                                                                                                                 <  1:719490/150‑1 (MQ=255)
                               cAGGTTCCCTGCCCCACCCCCTATGTCCGTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGcc                                                                                                                >  1:1775504/1‑150 (MQ=255)
                                  gTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGc                                                                                                                 >  2:1754066/1‑146 (MQ=255)
                                  gTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGc                                                                                                                 <  1:1754066/146‑1 (MQ=255)
                                    tCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGAGATGCCTGCACAGCAGGCCCCTTg                                                                                                           <  2:559059/150‑1 (MQ=255)
                                     cccTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCGTGCACAGCAGGCCCCTTgg                                                                                                          <  1:74114/150‑1 (MQ=255)
                                       cTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGCTGCCTGCACAGCAGGCCCCTTGGGc                                                                                                        >  1:2058019/1‑150 (MQ=255)
                                       cTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGc                                                                                                        <  2:960070/150‑1 (MQ=255)
                                        tGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGcc                                                                                                       <  1:1085692/150‑1 (MQ=255)
                                            ccACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAg                                                                                                   >  1:689140/1‑150 (MQ=255)
                                             cACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGa                                                                                                                                 >  1:1393385/1‑119 (MQ=255)
                                             cACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGa                                                                                                                                 <  2:1393385/119‑1 (MQ=255)
                                              aCCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGaa                                                                                                 >  2:1426473/1‑150 (MQ=255)
                                                  ccTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCgg                                                                                             >  1:715696/1‑150 (MQ=255)
                                                   cTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATAGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCggg                                                                                            >  2:1631102/1‑150 (MQ=255)
                                                    tATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGa                                                                                           >  2:2215020/1‑150 (MQ=255)
                                                    tATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGa                                                                                           <  2:1886738/150‑1 (MQ=255)
                                                        tCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATACCGTTTTCCCTATCTTGTGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCa                                                                                       <  2:1242961/150‑1 (MQ=255)
                                                          ccTACACCATTGGTCGGATGGCCACTACTGTGCGGTTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAgg                                                                                     >  2:271474/1‑150 (MQ=255)
                                                          ccTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCgg                                                                                             <  2:1382470/142‑1 (MQ=255)
                                                           cTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGTGACCCAGGc                                                                                    <  1:562764/150‑1 (MQ=255)
                                                           cTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGc                                                                                    >  1:1624301/1‑150 (MQ=255)
                                                            tACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCt                                                                                   >  2:399797/1‑150 (MQ=255)
                                                                   ttGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCaa                                                                            >  2:1109930/1‑150 (MQ=255)
                                                                    tGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCaaa                                                                           <  1:1426473/150‑1 (MQ=255)
                                                                         ggATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTc                                                                      <  1:1469648/150‑1 (MQ=255)
                                                                          gATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTAGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTcc                                                                     <  2:1745968/150‑1 (MQ=255)
                                                                               ccACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCt                                                                                                                                           >  2:2407709/1‑75 (MQ=255)
                                                                               ccACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCAtctc                                                                >  2:421506/1‑150 (MQ=255)
                                                                                cACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATctct                                                               >  1:330860/1‑150 (MQ=255)
                                                                                 actactGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCtt                                                              >  1:1568636/1‑150 (MQ=255)
                                                                                      tgtgCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGa                                                                                           >  1:1164074/1‑116 (MQ=255)
                                                                                      tgtgCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGa                                                                                           <  2:1164074/116‑1 (MQ=255)
                                                                                      tgtgCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTg                                                         >  1:2415340/1‑150 (MQ=255)
                                                                                           gATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTgat                                                    >  1:1900783/1‑148 (MQ=255)
                                                                                            aTTAATCACAAAGCATTAATCATGGGCCCAATAAAATTGGCCGCCTCCAATGGGAGCGGCCTTGCGTCGCACCGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTc                                                   >  1:654460/1‑150 (MQ=255)
                                                                                               aagcacaAAGCATTAATCATGCGCCGAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGTCCCTTGGGCCTTAGTACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGc                                                <  2:1003042/147‑1 (MQ=255)
                                                                                               aaTCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGc                                                <  2:613825/150‑1 (MQ=255)
                                                                                                 tCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGaa                                                                                                 >  2:918991/1‑99 (MQ=255)
                                                                                                   acaAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCt                                            >  1:886833/1‑150 (MQ=255)
                                                                                                    caAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTc                                           <  2:330860/150‑1 (MQ=255)
                                                                                                    caAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTc                                           >  2:125862/1‑150 (MQ=255)
                                                                                                     aaaGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTACGAACGTGGACCCAGGCACCAACAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTcc                                          >  1:2231678/1‑150 (MQ=255)
                                                                                                       aGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAg                                        >  2:1430774/1‑150 (MQ=255)
                                                                                                              aTCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAgg                                 >  1:1143350/1‑150 (MQ=255)
                                                                                                              aTCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAgg                                 <  2:1764789/150‑1 (MQ=255)
                                                                                                              aTCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAgg                                 >  2:1587984/1‑150 (MQ=255)
                                                                                                                   gcgcCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTtgcg                            <  1:399797/150‑1 (MQ=255)
                                                                                                                   gcgcCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTtgcg                            <  2:498200/150‑1 (MQ=255)
                                                                                                                     gcCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTg                                                                                                                           >  1:872536/1‑53 (MQ=255)
                                                                                                                     gcCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTtgcggt                          >  2:1444970/1‑150 (MQ=255)
                                                                                                                         aaTAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCATCTCCAGGTggcggtgcgg                      >  2:551951/1‑150 (MQ=255)
                                                                                                                            aaaaaaGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGaa                   >  2:185131/1‑150 (MQ=255)
                                                                                                                               aaaGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAgcg                <  1:271474/150‑1 (MQ=255)
                                                                                                                                aaGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTtgcggt                          <  1:473791/139‑1 (MQ=255)
                                                                                                                                aaGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAgcgc               <  1:125862/150‑1 (MQ=255)
                                                                                                                                  ggCCGCTCCCATTGGGGGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAGCCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGccc             >  1:1787830/1‑150 (MQ=255)
                                                                                                                                   gccgctCCCATTGGGGGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTgg                                                                                                          >  1:54407/1‑56 (MQ=25)
                                                                                                                                    ccgctCCCATTGGGGGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCAt                                                                   >  2:1161060/1‑94 (MQ=255)
                                                                                                                                     cgccCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGaa          >  2:350423/1‑150 (MQ=255)
                                                                                                                                       ctcccATTGGGGGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCATAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGAAGc        >  1:472338/3‑150 (MQ=255)
                                                                                                                                        ccccAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGAAGCg       <  2:1568636/150‑1 (MQ=255)
                                                                                                                                        ccccAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGAAGCg       <  2:2303351/150‑1 (MQ=255)
                                                                                                                                         cccAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGAAGCgg      <  1:1418085/150‑1 (MQ=255)
                                                                                                                                         cccAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGAAGCgg      <  2:708624/150‑1 (MQ=255)
                                                                                                                                           cAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGcc              >  2:51712/1‑140 (MQ=255)
                                                                                                                                           cAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGATCGGAAGATTGTTGAGCATTGTGTGGTTCTGGAAGCGCCCGAAGCGGGc    >  2:1841649/1‑150 (MQ=255)
                                                                                                                                             aTGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCAGAAGCGGGCaa  <  1:1323391/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                     
CTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGAAGCGGGCAA  >  NC_002947/5741423‑5741713

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: