Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,321,516 C→A 5.4% K407N (AAG→AAT) ‡ PP_1154 ← sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,321,5160CA 5.4% 435.2 / 8.7 130K407N (AAG→AAT) ‡PP_1154sensory box protein
Reads supporting (aligned to +/- strand):  ref base C (62/61);  new base A (2/5);  total (64/66)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 132 total aligned reads displayed.

CCGCGCAACAAGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTGGTATCGCA  >  NC_002947/1321367‑1321665
                                                                                                                                                     |                                                                                                                                                     
  gcgcAACAAGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCtt                                                                                                                                                     <  2:127016/150‑1 (MQ=255)
  gcgcAACAAGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCtt                                                                                                                                                     <  1:3935141/150‑1 (MQ=255)
  gcgcAACAAGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCtt                                                                                                                                                     >  2:1656727/1‑150 (MQ=255)
  gcgcAACAAGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCtt                                                                                                                                                     <  2:2418317/150‑1 (MQ=255)
  gcgcAACAAGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCAGGAAACCGATGGCCtt                                                                                                                                                     >  1:2426597/1‑150 (MQ=255)
  gcgcAACAAGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTGGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCtt                                                                                                                                                     <  2:4535383/150‑1 (MQ=255)
    gcAACAAGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGc                                                                                                                                                   <  2:3794332/150‑1 (MQ=255)
      aacaaGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCa                                                                                                                                                 >  2:4816341/1‑150 (MQ=255)
      aacaaGGCCTGCCCCAGGTGCCGGGCGAGCTGGTTGCGGTTGGCCAGTCCGGTCAGCGGGTCGTCCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGGGATCTCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCa                                                                                                                                                 <  2:969549/150‑1 (MQ=255)
       acaaGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGTAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCAt                                                                                                                                                >  1:809066/1‑150 (MQ=255)
         aaGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATcg                                                                                                                                              >  2:3078929/1‑150 (MQ=255)
         aaGGACTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGACGTGCAGCGAAAGGTGCTGAATGCGCGAGAGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATcg                                                                                                                                              >  2:4456154/1‑150 (MQ=255)
             ccTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGcc                                                                                                                                          >  2:1897673/1‑150 (MQ=255)
              cTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCa                                                                                                                                         <  2:1965858/150‑1 (MQ=255)
              cTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCa                                                                                                                                         >  2:1943652/1‑150 (MQ=255)
              cTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCAAGGCCATcggtttc                                                                                                                                         <  1:2752459/150‑6 (MQ=255)
              cTGCTCCAGGTGCCGGGCGAGCTTGGTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCa                                                                                                                                         <  2:3937625/150‑1 (MQ=255)
              cTGCTCCAGGGGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGACAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCAAGGCCATcggtttc                                                                                                                                         <  1:194282/150‑6 (MQ=255)
                  tCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATa                                                                                                                                     >  1:1845016/1‑150 (MQ=255)
                     aGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAgcg                                                                                                                                  >  1:3777477/1‑150 (MQ=255)
                           cGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCgcgggcg                                                                                                                            <  2:2137488/150‑1 (MQ=255)
                            gggCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCgcgggcgg                                                                                                                           >  1:2313987/1‑150 (MQ=255)
                            gggCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCgcgggcgg                                                                                                                           >  2:3890699/1‑150 (MQ=255)
                            gggCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCgcgggcgg                                                                                                                           >  1:2968576/1‑150 (MQ=255)
                             ggCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTTCCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGa                                                                                                                          <  1:4383042/150‑1 (MQ=255)
                             ggCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGa                                                                                                                          <  1:4383047/150‑1 (MQ=255)
                                   cTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCg                                                                                                                    >  1:1183054/1‑150 (MQ=255)
                                    ttgttgCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAATGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCgg                                                                                                                   >  2:4254566/1‑150 (MQ=255)
                                    ttgttgCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGTCTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCgg                                                                                                                   >  2:4254318/1‑150 (MQ=255)
                                    ttgttgCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCgg                                                                                                                   <  2:1702112/150‑1 (MQ=255)
                                    ttgttgCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCgg                                                                                                                   <  2:1702035/150‑1 (MQ=255)
                                    ttgttgCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCgg                                                                                                                   >  2:4252063/1‑150 (MQ=255)
                                        tgCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGt                                                                                                               >  1:3850068/1‑150 (MQ=255)
                                        tgCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGt                                                                                                               >  1:624437/1‑150 (MQ=255)
                                         gCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTa                                                                                                              <  1:1663908/150‑1 (MQ=255)
                                         gCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTa                                                                                                              <  2:1537883/150‑1 (MQ=255)
                                          cGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAg                                                                                                             <  1:1115865/150‑1 (MQ=255)
                                          cGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAg                                                                                                             <  1:1115784/150‑1 (MQ=255)
                                                  cAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCTCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCTGTACACCGGAAAACGATTGCCTTGCCATCGCGCCAGATAGCTCGGGCGGATATCCGGCAGTAGCGGTTCTg                                                                                                     <  1:135368/150‑1 (MQ=255)
                                                  cAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTg                                                                                                     <  1:4966583/150‑1 (MQ=255)
                                                          tCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCt                                                                                             <  2:4771478/150‑1 (MQ=255)
                                                           cAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTg                                                                                            <  2:2637085/150‑1 (MQ=255)
                                                           cAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTg                                                                                            <  2:2636823/150‑1 (MQ=255)
                                                           cAGTGGGTCGTACAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTg                                                                                            <  1:4326850/150‑1 (MQ=255)
                                                             gTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGt                                                                                          <  2:2482591/150‑1 (MQ=255)
                                                                ggTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTcgc                                                                                       >  2:2101294/1‑150 (MQ=255)
                                                                ggTCGTGCAGCGAAAGGCGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTcgc                                                                                       >  2:3923707/1‑150 (MQ=255)
                                                                ggTCGTGCAGCGAAAGGCGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTcgc                                                                                       >  2:3924105/1‑150 (MQ=255)
                                                                     tGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTaa                                                                                  <  1:4708935/150‑1 (MQ=255)
                                                                       cAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTGCTGGCCATGCAGGTCGCGGTAAGt                                                                                >  2:3391910/1‑150 (MQ=255)
                                                                        aGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTa                                                                               <  2:4955891/150‑1 (MQ=255)
                                                                          cGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACa                                                                             >  2:5047285/1‑150 (MQ=255)
                                                                          cGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACa                                                                             <  2:3248375/150‑1 (MQ=255)
                                                                          cGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACa                                                                             >  2:5048447/1‑150 (MQ=255)
                                                                           gAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAg                                                                            >  1:165902/1‑150 (MQ=255)
                                                                            aaaGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAgc                                                                           >  1:532211/1‑150 (MQ=255)
                                                                                    tGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTg                                                                   >  2:5111401/1‑150 (MQ=255)
                                                                                      aaTGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGc                                                                 <  2:2679059/150‑1 (MQ=255)
                                                                                      aaTGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGc                                                                 <  1:1897673/150‑1 (MQ=255)
                                                                                        tGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCa                                                               >  1:3349954/1‑150 (MQ=255)
                                                                                                     gTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGcggctc                                                  >  2:5044387/1‑150 (MQ=255)
                                                                                                     gTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGcggctc                                                  <  1:788964/150‑1 (MQ=255)
                                                                                                       cGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGcggctcgg                                                <  2:4337871/150‑1 (MQ=255)
                                                                                                       cGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGcggctcgg                                                <  1:3421385/150‑1 (MQ=255)
                                                                                                         aCTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGcggctcggct                                              <  1:2702507/150‑1 (MQ=255)
                                                                                                          cTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTg                                             <  2:2872049/150‑1 (MQ=255)
                                                                                                            tcgtcgGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCa                                           <  2:233824/150‑1 (MQ=255)
                                                                                                              gtcgGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGt                                         <  1:9680/150‑1 (MQ=255)
                                                                                                               tcgGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTc                                        >  2:2827411/1‑150 (MQ=255)
                                                                                                               tcgGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTc                                        >  2:2827189/1‑150 (MQ=255)
                                                                                                                    gATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCCTTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGc                                   >  2:3131036/1‑150 (MQ=255)
                                                                                                                    gATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGc                                   >  2:1574141/1‑150 (MQ=255)
                                                                                                                    gATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGc                                   >  2:4675163/1‑150 (MQ=255)
                                                                                                                     atatCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCg                                  >  1:2065550/1‑150 (MQ=255)
                                                                                                                     atatCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCg                                  >  1:2065552/1‑150 (MQ=255)
                                                                                                                      tatCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGt                                 >  1:3182593/1‑150 (MQ=255)
                                                                                                                      tatCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGt                                 >  1:3182513/1‑150 (MQ=255)
                                                                                                                       atCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTc                                >  2:5120758/1‑150 (MQ=255)
                                                                                                                         cGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCa                              <  2:652367/150‑1 (MQ=255)
                                                                                                                          gCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAg                             <  1:1656727/150‑1 (MQ=255)
                                                                                                                              gCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCaa                         <  1:4257018/150‑1 (MQ=255)
                                                                                                                               cGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAg                        <  1:2147752/150‑1 (MQ=255)
                                                                                                                                ggTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAgg                       >  2:1324060/1‑150 (MQ=255)
                                                                                                                                ggTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAgg                       >  1:2416585/1‑150 (MQ=255)
                                                                                                                                     cccGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCTGTCAGGGCGTCCAGCCAAGGTGACa                  >  1:2452760/1‑150 (MQ=255)
                                                                                                                                     cccGGAAACCGATGGCCGTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACa                  <  2:679267/150‑1 (MQ=255)
                                                                                                                                      ccGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTCAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAg                 >  2:4991051/1‑150 (MQ=255)
                                                                                                                                      ccGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAg                 >  1:2463931/1‑150 (MQ=255)
                                                                                                                                       cGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAgc                >  2:3392128/1‑150 (MQ=255)
                                                                                                                                       cGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGAAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAgc                >  1:2456916/1‑150 (MQ=255)
                                                                                                                                          aaaCCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAgcggc             <  1:4254318/150‑1 (MQ=255)
                                                                                                                                          aaaCCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAgcggc             <  1:4254566/150‑1 (MQ=255)
                                                                                                                                            aCCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGt           >  1:3157782/1‑150 (MQ=255)
                                                                                                                                            aCCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGt           >  1:601518/1‑150 (MQ=255)
                                                                                                                                             gcGATGGCAAGGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGGTGCGGCTCGGCTGCAGTCAGGGCGTCGAGCCAAGGTGACAGCGGCGTg          >  2:856688/2‑150 (MQ=255)
                                                                                                                                             gcGATGGCAAGGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTg          >  1:314820/2‑150 (MQ=255)
                                                                                                                                             ccGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTg          <  1:4456154/150‑1 (MQ=255)
                                                                                                                                              cGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTgg         >  2:4850559/1‑150 (MQ=255)
                                                                                                                                                   gCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTGGTATCg    >  2:2097387/1‑150 (MQ=255)
                                                                                                                                                     cTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTGGTATCGCa  <  1:3078929/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCGCAACAAGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTGGTATCGCA  >  NC_002947/1321367‑1321665

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