Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,449,296 G→C 8.4% R165T (AGG→ACG) ‡ PP_3067 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,449,2960GC 8.4% 215.0 / 9.2 72R165T (AGG→ACG) ‡PP_3067hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/34);  new base C (4/2);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

TTTGGGGGAATTTGAATCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTGAAGATGCAAACAAACACAAAGGTTCG  >  NC_002947/3449150‑3449441
                                                                                                                                                  |                                                                                                                                                  
tttGGGGGAATTTGAATCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAgg                                                                                                                                                   <  3:358917/148‑1 (MQ=255)
   gggggAATTTGAATCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGAc                                                                                                                                                <  3:162979/148‑1 (MQ=255)
   gggggAATTTGAATCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGAc                                                                                                                                                >  7:53305/1‑148 (MQ=255)
   gggggAATTTGAATCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCAATCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGAc                                                                                                                                                <  3:331847/148‑1 (MQ=255)
               aTCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGAGAGCGTCCCTGACGCTGTCGCtcacc                                                                                                                                   <  5:112415/149‑5 (MQ=255)
                    ttACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTc                                                                                                                              <  8:53305/149‑1 (MQ=255)
                        tAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAgg                                                                                                                          >  1:94696/1‑149 (MQ=255)
                          aaCATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGa                                                                                                                            >  1:131807/1‑145 (MQ=255)
                          aaCATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGa                                                                                                                            <  2:131807/145‑1 (MQ=255)
                             aTCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAg                                                                                                                     <  2:94696/149‑1 (MQ=255)
                               cccTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGATCGAGGGGAAggt                                                                                                                   >  2:192311/1‑149 (MQ=255)
                                    ccTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAggtggg                                                                                                                <  5:29658/147‑1 (MQ=255)
                                          tAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCTCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGCTGGGTGGGGaa                                                                                                         >  1:49787/1‑148 (MQ=255)
                                            aGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAatat                                                                                                      <  2:232030/149‑1 (MQ=255)
                                                 ttGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTTTGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAggtgggtg                                                                                                              >  1:360771/1‑136 (MQ=255)
                                                 ttGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAggtgggtg                                                                                                              <  2:360771/136‑1 (MQ=255)
                                                      ttCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAaga                                                                                            <  5:34361/149‑1 (MQ=255)
                                                      ttCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAaga                                                                                            <  7:20706/149‑1 (MQ=255)
                                                       tCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGgtgggtgg                                                                                                             <  4:356587/131‑1 (MQ=255)
                                                         aCGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGc                                                                                         <  5:26608/149‑1 (MQ=255)
                                                          cGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCCCTGACGCTGTCGCtc                                                                                                                                      >  3:287808/1‑102 (MQ=255)
                                                          cGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCCCTGACGCTGTCGCtc                                                                                                                                      <  4:287808/103‑2 (MQ=255)
                                                            ttttCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTtg                                                                                      >  4:192199/1‑149 (MQ=255)
                                                               tcttcAGGGGCCCAACAGGGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTtgat                                                                                    <  4:57945/148‑1 (MQ=255)
                                                                      gtgtCCAACAGTGTCATTGATTAACTTCAAGCGAAAGGTGGCGCTTATGATGAAATGGTGGTGAGCGAGAGCGTCAGGGACGCTCTCGCTGAATTTCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATGAt                                                                            <  4:103049/149‑1 (MQ=255)
                                                                         tcgaaCAGTGTCATTGAGGAACTTCAAGCGAAAGGTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTaa                                                                         <  2:49787/146‑1 (MQ=255)
                                                                         tCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGAGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTa                                                                          <  4:79998/148‑1 (MQ=255)
                                                                            aaCAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTtgat                                                                                    <  5:244168/135‑1 (MQ=255)
                                                                            aaCAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTtgat                                                                                    >  6:244168/1‑135 (MQ=255)
                                                                             actgtgtCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATTATGAAATGGTTGTTAGCGACAGCGTCAGGGACGCTCTCGCTGAATTTCTCGAGGTGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTaaggaa                                                                     <  5:348893/146‑1 (MQ=255)
                                                                               aGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCAtt                                                                           <  8:202460/141‑1 (MQ=255)
                                                                               aGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCAtt                                                                           >  7:202460/1‑141 (MQ=255)
                                                                                gtgtCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAAAGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGCGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAggg                                                                  >  7:186759/1‑149 (MQ=255)
                                                                                        gAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTc                                                             >  5:207319/1‑146 (MQ=255)
                                                                                         aGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTc                                                             >  5:77762/1‑145 (MQ=255)
                                                                                         aGGAACTTCAAGCGAAAGCTGCCGATTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTc                                                             <  6:77762/145‑1 (MQ=255)
                                                                                          ggAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGAt                                                        <  3:192199/149‑1 (MQ=255)
                                                                                             aCTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACg                                                     >  6:243108/1‑149 (MQ=255)
                                                                                              cTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACg                                                     >  6:325830/1‑148 (MQ=255)
                                                                                               ttCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGAt                                                        <  2:216016/144‑1 (MQ=255)
                                                                                               ttCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGAt                                                        >  1:216016/1‑144 (MQ=255)
                                                                                               ttCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGc                                                    >  5:168201/1‑148 (MQ=255)
                                                                                                    gCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCa                                              >  1:166136/1‑149 (MQ=255)
                                                                                                       aaaGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCatat                                           >  6:322357/1‑149 (MQ=255)
                                                                                                        aaGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGa                                                                >  7:279584/1‑127 (MQ=255)
                                                                                                        aaGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGa                                                                <  8:279584/127‑1 (MQ=255)
                                                                                                         aGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATtc                                         <  2:166136/149‑1 (MQ=255)
                                                                                                            tGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATtcctc                                     >  5:93275/1‑149 (MQ=255)
                                                                                                             gCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCa                                                                             >  3:368111/1‑109 (MQ=255)
                                                                                                             gCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCa                                                                             <  4:368111/109‑1 (MQ=255)
                                                                                                              ccGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATtcctc                                     <  5:322357/147‑1 (MQ=255)
                                                                                                                  ttATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGAt                               >  2:236244/1‑149 (MQ=255)
                                                                                                                      gatgaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGTTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGc                                                    >  3:118389/1‑125 (MQ=255)
                                                                                                                      gatgaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGTTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGc                                                    <  4:118389/125‑1 (MQ=255)
                                                                                                                      gatgaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATAc                                                      >  6:101315/1‑123 (MQ=255)
                                                                                                                      gatgaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATAc                                                      <  5:101315/123‑1 (MQ=255)
                                                                                                                       atgaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTGaa                          <  1:192311/149‑1 (MQ=255)
                                                                                                                        tgaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTg                            >  3:239327/1‑146 (MQ=255)
                                                                                                                         gaAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTg                            >  2:164645/1‑145 (MQ=255)
                                                                                                                             tggtggTGAGCGACATCTTCAGGGACGCTCTCGCTTAATTGCTCGAGGGGAAGGTGGGTGGGGAATATTATAAACAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTGAAGATGca                    <  4:239327/149‑1 (MQ=255)
                                                                                                                                tggtgAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAATTCATATTCCTCATGGATGTGAAGAGGaaaa                  >  4:248954/1‑148 (MQ=255)
                                                                                                                                   tgAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTGAAGATGcaaacaaa              <  1:164645/149‑1 (MQ=255)
                                                                                                                                      gCGAGAGCGTCCCTGACGCTGTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCATCCTGGATTTGAAGATGCAaacaaacaa           >  7:329372/1‑148 (MQ=255)
                                                                                                                                       cGAGAGCGTCCCTGACGCTGTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTGAAGATGCAAACAAacac           >  1:10216/1‑148 (MQ=255)
                                                                                                                                       cGAGAGCGTCCCTGACGCTGTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTGAAGATGCAAACAAacac           >  8:5524/1‑148 (MQ=255)
                                                                                                                                       cGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTGAAGATGCAAACAAacac           >  5:60692/1‑148 (MQ=255)
                                                                                                                                       cGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTGAAGATGCAAACAAacac           >  4:48555/1‑148 (MQ=255)
                                                                                                                                        gACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTGAAGATGCAAACAAACACaa         <  6:93275/149‑1 (MQ=255)
                                                                                                                                              gTCAGGGACGCTCTCGCTGAATTTCTCGAGGTGAAGGTGGGTTTTTAATATGATAAAAAGAGGCTTTATGAAATCATTAAGGAAGGTGAGTATAGATACGCAAATCATATTCCTCCTGGATTTGAAGATTCAAACAATCACAAAGTTTc   <  6:168201/149‑1 (MQ=255)
                                                                                                                                              gTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTGAAGATGcaaac                 >  4:203495/1‑135 (MQ=255)
                                                                                                                                              gTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTGAAGATGcaaac                 <  3:203495/135‑1 (MQ=255)
                                                                                                                                               tCATGGAAGCTCTTGCTTTATTTCTAGATGGGAATGTGGGTGGGGAATATGATAAAAATAGGCTTTATTAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCCTGGATTTGAAGATGCAAACAAACACAAAGGTTc   <  1:236244/149‑1 (MQ=255)
                                                                                                                                               tCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGGTTTGAAGATGCAAACAAACACAAAGGTTCg  >  8:325355/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TTTGGGGGAATTTGAATCAATTACTAAACATCCCTTCCTCAGTAAGGACTTGCATTCACGTTTTCTTCAGGTGTCCAACAGTGTCATTGAGGAACTTCAAGCGAAAGCTGCCGCTTATGATGAAATGGTGGTGAGCGACAGCGTCAGGGACGCTCTCGCTGAATTGCTCGAGGGGAAGGTGGGTGGGGAATATGATAAAAAGAGGCTTGATGAAATCATTAAGGAAGGGGAGTCTAGATACGCAAATCATATTCCTCCTGGATTTGAAGATGCAAACAAACACAAAGGTTCG  >  NC_002947/3449150‑3449441

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