Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 160,519 C→A 5.4% A63S (GCA→TCA) ‡ PP_0150 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947160,5190CA 5.4% 306.6 / 5.0 93A63S (GCA→TCA) ‡PP_0150hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (42/46);  new base A (3/2);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

CACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTCGGCCAG  >  NC_002947/160373‑160663
                                                                                                                                                  |                                                                                                                                                 
cacgcacgCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAg                                                                                                                                                 <  5:241337/149‑1 (MQ=255)
 acgcacgcGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGc                                                                                                                                                >  5:162011/1‑149 (MQ=255)
 acgcacgcGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGc                                                                                                                                                >  3:358120/1‑149 (MQ=255)
   gcacgcGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTAGTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGcgt                                                                                                                                              >  4:234970/1‑149 (MQ=255)
         gcTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGACGCTGCAAGcttcgtcgc                                                                                                                                        <  3:342/149‑2 (MQ=255)
         gcTTCACAGTGCGAGGATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGcgtcgtcgt                                                                                                                                        <  7:58490/149‑1 (MQ=255)
          cTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTg                                                                                                                                       <  1:186051/149‑1 (MQ=255)
           ttCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTgg                                                                                                                                      >  6:237234/1‑149 (MQ=255)
           ttCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTgg                                                                                                                                      >  4:84047/1‑149 (MQ=255)
           ttCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTgg                                                                                                                                      >  4:116552/1‑149 (MQ=255)
              acaGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGACGCTGCAAGcttcgtcgc                                                                                                                                        <  7:268111/144‑2 (MQ=255)
              acaGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGACGCTGCAAGcttcgtcgc                                                                                                                                        >  8:268111/1‑143 (MQ=255)
              acaGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGAAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGATGGCGACGAAGCTTGCAGCGTCGTCGTGGGCt                                                                                                                                   >  3:248018/1‑149 (MQ=255)
                  tGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTcgc                                                                                                                               >  2:104796/1‑149 (MQ=255)
                     gAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGca                                                                                                                            <  3:116552/149‑1 (MQ=255)
                     gAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGca                                                                                                                            <  6:322565/149‑1 (MQ=255)
                           ggCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTAGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTAGCAGCGTCGTCGTGGGCTTCGCGCACCACtt                                                                                                                      >  5:197208/1‑149 (MQ=255)
                            gCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACtt                                                                                                                      <  4:319038/148‑1 (MQ=255)
                             cTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTg                                                                                                                    <  6:52450/149‑1 (MQ=255)
                             cTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTg                                                                                                                    <  3:234970/149‑1 (MQ=255)
                             cTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTg                                                                                                                    <  3:169354/149‑1 (MQ=255)
                             cTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTg                                                                                                                    >  7:39573/1‑149 (MQ=255)
                             cTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTg                                                                                                                    <  1:104796/149‑1 (MQ=255)
                             cTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTg                                                                                                                    >  8:221580/1‑149 (MQ=255)
                              ttCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTg                                                                                                                    >  8:127482/1‑148 (MQ=255)
                                cAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGa                                                                                                                  <  1:279629/148‑1 (MQ=255)
                                   ccAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcg                                                                                                              >  8:241965/1‑149 (MQ=255)
                                    cAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGcgtcgtcgt                                                                                                                                        <  2:305615/122‑1 (MQ=255)
                                    cAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGcgtcgtcgt                                                                                                                                        >  1:305615/1‑122 (MQ=255)
                                    cAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcgc                                                                                                             <  6:162011/149‑1 (MQ=255)
                                          cTGGTCGCTGGCGCCGCGCGCTTCCGGTTCTCGCTCGCGCAGGGTATTGAAGGACGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATc                                                                                                        <  7:241965/148‑1 (MQ=255)
                                              tCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggt                                                                                                   >  6:40087/1‑149 (MQ=255)
                                                 cTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggtggt                                                                                                <  3:84047/149‑1 (MQ=255)
                                                  tGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggtggtg                                                                                               >  1:181976/1‑149 (MQ=255)
                                                  tGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggtggtg                                                                                               <  3:318585/149‑1 (MQ=255)
                                                          cgcgcTTCCGGTACTCGCTGGCGCAGGGTTTTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGTGGAGGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTTGTG‑GTGTCTGa                                                                                      <  6:197208/149‑1 (MQ=255)
                                                           gcgcTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAg                                                                                     <  5:237234/149‑1 (MQ=255)
                                                            cgcTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggtggtg                                                                                               >  3:246674/1‑139 (MQ=255)
                                                            cgcTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggtggtg                                                                                               <  4:246674/139‑1 (MQ=255)
                                                            cgcTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAgg                                                                                    >  3:368091/1‑149 (MQ=255)
                                                               ttCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTg                                                                                  >  4:145924/1‑148 (MQ=255)
                                                                 ccGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGa                                                                                 >  5:244798/1‑147 (MQ=255)
                                                                 ccGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGa                                                                                 <  6:244798/147‑1 (MQ=255)
                                                                 caggTACTCGCTCGCGCTGGGTAGTGAAGGCCGCCTCGATTTGTACGCCGCGGATGTCACAGTTGCTGGAGACGAATCTTGCATCTTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCTCGGATCGAGGTTGTG‑GTGTCTGAgtttaa                                                                                <  4:248018/146‑6 (MQ=255)
                                                                  cGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCAGCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGt                                                                              <  4:358120/149‑1 (MQ=255)
                                                                  cGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGt                                                                              >  2:194122/1‑149 (MQ=255)
                                                                  cGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGt                                                                              >  7:44/1‑149 (MQ=255)
                                                                    gTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGtc                                                                             >  7:242863/1‑148 (MQ=255)
                                                                       ctcgctcgCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGtcgat                                                                          <  1:194122/148‑1 (MQ=255)
                                                                          gctcgcGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcgc                                                                                                             >  2:150701/1‑111 (MQ=255)
                                                                          gctcgcGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcgc                                                                                                             <  1:150701/111‑1 (MQ=255)
                                                                          gctcgcGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGtcgatcg                                                                        <  4:368091/147‑1 (MQ=255)
                                                                              gcgcAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGc                                                                  >  2:281246/1‑149 (MQ=255)
                                                                                gcAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑gtgt                                                                                          >  7:296477/1‑123 (MQ=255)
                                                                                gcAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑gtgt                                                                                          <  8:296477/123‑1 (MQ=255)
                                                                                   gggTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAg                                                             <  5:40087/149‑1 (MQ=255)
                                                                                   gggTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAg                                                             <  8:114615/149‑1 (MQ=255)
                                                                                         tGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGgcacgc                                                       >  8:117283/1‑149 (MQ=255)
                                                                                           aaGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGgcacgcac                                                     >  2:87473/1‑149 (MQ=255)
                                                                                           aaGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GGGTCTGAGGTGAAGTCGATCGAACGGCCAAAGgcacgcac                                                     >  1:124797/1‑149 (MQ=255)
                                                                                           aaGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTC‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGgcacgcac                                                     >  5:317731/1‑149 (MQ=255)
                                                                                                cgccTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACgatga                                                <  5:142747/149‑1 (MQ=255)
                                                                                                 gccTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACgatgat                                               <  1:87473/149‑1 (MQ=255)
                                                                                                  ccTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGatgatg                                              <  6:317731/149‑1 (MQ=255)
                                                                                                   cTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGt                                             >  1:44612/1‑149 (MQ=255)
                                                                                                     cGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGt                                                                              <  8:205942/114‑1 (MQ=255)
                                                                                                      gAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGaa                                                                                <  4:218416/111‑1 (MQ=255)
                                                                                                      gAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGaa                                                                                >  3:218416/1‑111 (MQ=255)
                                                                                                      gAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGtt                                            <  3:152506/147‑1 (MQ=255)
                                                                                                      gAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGc                                          <  8:39573/149‑1 (MQ=255)
                                                                                                        gTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCtg                                        <  7:221580/149‑1 (MQ=255)
                                                                                                         ttGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGCGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCtg                                        <  1:281246/149‑1 (MQ=255)
                                                                                                         ttGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCtg                                        <  2:181976/148‑1 (MQ=255)
                                                                                                           gTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCtggt                                      <  7:127482/148‑1 (MQ=255)
                                                                                                            tACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCtggtg                                     <  5:228457/148‑1 (MQ=255)
                                                                                                              cgccgcGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGc                                   >  7:259589/1‑148 (MQ=255)
                                                                                                                ccgcGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGccac                                >  3:32078/1‑149 (MQ=255)
                                                                                                                ccgcGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGccac                                >  1:254082/1‑149 (MQ=255)
                                                                                                                           aCCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGTGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCtgtg                     >  1:202521/1‑149 (MQ=255)
                                                                                                                            ccGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGc                    <  4:32078/149‑1 (MQ=255)
                                                                                                                             cGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggtggtg‑gt                                                                                            <  7:58951/76‑1 (MQ=255)
                                                                                                                             cGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggtggtg‑gt                                                                                            >  8:58951/1‑76 (MQ=255)
                                                                                                                                 gCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTCGCTGGTGGCCACCAGGGTCTGTGCCTgg                <  3:145924/148‑1 (MQ=255)
                                                                                                                                   ccacgacgaCGCTGCAAGCTTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGc              >  4:121276/4‑148 (MQ=255)
                                                                                                                                     gCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAg           <  8:242863/149‑1 (MQ=255)
                                                                                                                                     acgacgaCGCTGCAAGCTTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAg           >  8:289444/2‑149 (MQ=255)
                                                                                                                                       gacgaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCtggtg                                     >  3:188777/1‑121 (MQ=255)
                                                                                                                                       gacgaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCtggtg                                     <  4:188777/121‑1 (MQ=255)
                                                                                                                                       gacgaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGt         >  2:195807/1‑149 (MQ=255)
                                                                                                                                       gacgaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGGCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGt         >  2:159414/1‑149 (MQ=255)
                                                                                                                                        acgaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTc        <  8:44/149‑1 (MQ=255)
                                                                                                                                            aGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTCGGc     >  7:263459/1‑148 (MQ=255)
                                                                                                                                               ttGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTCGGCcag  >  4:132173/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTCGGCCAG  >  NC_002947/160373‑160663

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