Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,891,050 A→G 5.0% S12G (AGC→GGC)  PP_3435 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,891,0500AG 5.0% 242.5 / 2.1 79S12G (AGC→GGC) PP_3435hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (39/36);  new base G (2/2);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

ATTCGTAAATCGATTCTACATTGCGTAAAATGGTTTGTCGCGTTTTAGTTGCTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAGCCACGCCACCACCCGCCAGGCGTTGGAGCACGTCGAACACATCCTCGA  >  NC_002947/3890907‑3891198
                                                                                                                                                |                                                                                                                                                     
aTTCGTAAATCGATTCTACATTGCGTAAAATGGTTTGTCGCGTTTTAGTTGCTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTgg                                                                                                                                                  >  6:58891/1‑149 (MQ=255)
 ttCGTAAATCGATTCTACATTGCGTAAAATGGTTTGTCGCGTTTTAGTTGCTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTgg                                                                                                                                                  <  5:80605/148‑1 (MQ=255)
    gTAAATCGATTCTACATTGCGTAAAATGGTTTGTCGCGTTTTAGTTGCTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCa                                                                                                                                              >  3:329335/1‑149 (MQ=255)
       aaTCGATTCTACATTGCGTAAAATGGTTTGTCGCGTTTTAGTTGCTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGc                                                                                                                                          >  5:234807/1‑149 (MQ=255)
         tCGATACTACATTGCGTAAAATGGTTTGTCGCGTTTTAGTTGCTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGctc                                                                                                                                        <  1:359206/149‑1 (MQ=255)
             ttCTACATTGCGTAAAATGGTTTGTCGCGTTTTAGTTGCTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGc                                                                                                                                    >  3:225653/1‑149 (MQ=255)
                tACATTGCGTAAAATGGTTTGTCGCGTTTTAGTTGCTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCAAGCTCTTGCCCt                                                                                                                                 >  5:170431/1‑149 (MQ=255)
                 aCATTGCGTAAAATGGTTTGTCGCGTTTTAGTTGCTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTg                                                                                                                                >  1:172892/1‑149 (MQ=255)
                 aCATTGCGTAAAATGGTTTGTCGCGTTTTAGTTGCTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTg                                                                                                                                >  1:340595/1‑149 (MQ=255)
                  cATTGCGTAAATTGGTTTGTCGCGTTTTAGTTGCTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCt                                                                                                                                 <  3:70680/147‑1 (MQ=255)
                     tGCGTAAAATGGTTTGTCGCGTTTTAGTTGCTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGgtc                                                                                                                            >  4:124898/1‑149 (MQ=255)
                          cacgggacaggttgtcgaggatgtgttcgacgtgctccaacgcctggcgggtggtggcgtggctgctggcctgcatttcacgttcgacctgccagttcatcacggccagGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGtcgtcg                                                                                                                        <  6:123122‑M1/40‑1 (MQ=255)
                            aaTGGTTTGTCGCGTTTTAGTTGCTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCg                                                                                                                     <  8:225070/149‑1 (MQ=255)
                             aTGGTTTGTCGCGTTTTAGTTGCTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGt                                                                                                                    <  2:23420/149‑1 (MQ=255)
                                gTTTGTCGCGTTTTAGTTGCTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGtgct                                                                                                                 >  6:121436/1‑149 (MQ=255)
                                gTTTGTCGCGTTTTAGTTGCTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGtgc                                                                                                                  <  4:329335/148‑1 (MQ=255)
                                 tttGTCGCGTTTTAGTTGCTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGtgctg                                                                                                                >  4:73085/1‑149 (MQ=255)
                                     tCGCGTTTTAGTTGCTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCg                                                                                                             <  3:124898/148‑1 (MQ=255)
                                          gaggatgtgttcgacgtgctccaacgcctggcgggtggtggcgtggctgctggcctgcatttcacgttcgacctgccagttcatcacggccagGCCGCACATGACCGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCg                                                                                                             >  4:1791‑M1/94‑144 (MQ=255)
                                          gaggatgtgttcgacgtgctccaacgcctggcgggtggtggcgtggctgctggcctgcatttcacgttcgacctgccagttcatcacggccagGCCGCACATGACCGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCg                                                                                                             <  3:1791‑M1/51‑1 (MQ=255)
                                               gTTGCTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCgg                                                                                                  >  3:260492/1‑149 (MQ=255)
                                               gTTGCTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCg                                                                                                   <  6:234807/148‑1 (MQ=255)
                                                 tgttcgacgtgctccaacgcctggcgggtggtggcgtggctgctggcctgcatttcacattcgacctgccagttcatcacggccagGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGc                                                                                                 >  6:133253‑M1/87‑149 (MQ=255)
                                                   ttcgacgtgctccaacgcctggcgggtggtggcgtggctgctggcctgcatttcacgttcgacctgccagttcatcacggccagGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGcc                                                                                                >  5:123122‑M1/85‑148 (MQ=255)
                                                   cTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCGGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTg                                                                                              >  2:89838/1‑149 (MQ=255)
                                                   cTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCt                                                                                               <  2:172892/148‑1 (MQ=255)
                                                    tCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTgg                                                                                             >  4:14938/1‑149 (MQ=255)
                                                      gTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTgg                                                                                             <  5:235413/147‑1 (MQ=255)
                                                       tACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGcc                                                                                           <  1:261280/149‑1 (MQ=255)
                                                         cgcgGCGACAAACGGCTATCCTCTCGTTTTAACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCGCGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGtg                                                                                        <  7:253148/149‑1 (MQ=255)
                                                          gcggcgACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCGGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGtga                                                                                       <  1:89838/149‑1 (MQ=255)
                                                          gcggcgACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGtga                                                                                       <  5:121436/149‑1 (MQ=255)
                                                           cggcgACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGtgat                                                                                      <  8:186483/149‑1 (MQ=255)
                                                           cggcgACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGtgat                                                                                      <  4:260492/149‑1 (MQ=255)
                                                           cggcgACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGtgat                                                                                      >  1:61724/1‑149 (MQ=255)
                                                               gACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATCCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTCCACTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATCaa                                                                                   >  2:300638/1‑148 (MQ=255)
                                                                  aaaCGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTgg                                                                               >  2:83990/1‑149 (MQ=255)
                                                                              tctcGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTGCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACg                                                                    >  1:78617/1‑149 (MQ=255)
                                                                               ctcGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGt                                                                   >  5:25114/1‑148 (MQ=255)
                                                                                tcGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGtga                                                                                       >  7:300039/1‑127 (MQ=255)
                                                                                tcGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGtga                                                                                       <  8:300039/127‑1 (MQ=255)
                                                                                tcGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGt                                                                   <  8:157023/147‑1 (MQ=255)
                                                                                tcGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGt                                                                   >  7:157023/1‑147 (MQ=255)
                                                                                  gTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTg                                                                  <  5:58891/146‑1 (MQ=255)
                                                                                  gTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGaaa                                                               >  6:247071/1‑149 (MQ=255)
                                                                                  gTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGa                                                                 <  1:46893/147‑1 (MQ=255)
                                                                                      ttACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGc                                                            <  4:225653/148‑1 (MQ=255)
                                                                                        aCTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAg                                                          <  6:25114/148‑1 (MQ=255)
                                                                                          tCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAg                                                                            >  4:255830/1‑128 (MQ=255)
                                                                                            cATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGc                                                                                            >  4:234463/1‑110 (MQ=255)
                                                                                            cATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGc                                                                                            <  3:234463/110‑1 (MQ=255)
                                                                                            cATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCcag                                                     >  3:76567/1‑149 (MQ=255)
                                                                                             aTGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAAt                                                              >  5:131083/1‑139 (MQ=255)
                                                                                             aTGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAAt                                                              <  6:131083/139‑1 (MQ=255)
                                                                                               gATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGc                                                            <  1:219711/139‑1 (MQ=255)
                                                                                               gATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGc                                                            >  2:219711/1‑139 (MQ=255)
                                                                                                aTAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCcagcag                                                  <  1:212280/148‑1 (MQ=255)
                                                                                                      aGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAgccacgcc                                           >  5:280924/1‑149 (MQ=255)
                                                                                                       gATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAgccacgcca                                          <  3:73085/149‑1 (MQ=255)
                                                                                                            gCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAGCCACGCCaccacc                                     <  1:83990/149‑1 (MQ=255)
                                                                                                               tGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAGCCACGCCACCACCCGc                                  <  2:78617/149‑1 (MQ=255)
                                                                                                                 ccccTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCCGGCCAGCAGCCACGCCACCACCCGCCa                                <  6:170431/149‑1 (MQ=255)
                                                                                                                  cccTTACCGTCAAACGCCGCGCCGCACTGAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAGCCACGCCACCACCCGCCAg                               >  1:205470/1‑149 (MQ=255)
                                                                                                                  cccTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAGCCACGCCACCACCCGCCAg                               >  1:174549/1‑149 (MQ=255)
                                                                                                                  cccTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAGCCACGCCACCACCCGCCAg                               >  4:207940/1‑149 (MQ=255)
                                                                                                                         cGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAGCCACGCCACCACCCGCCAGGCGTTg                         >  3:177166/1‑148 (MQ=255)
                                                                                                                          gTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAgccac                                              >  7:140541/1‑126 (MQ=255)
                                                                                                                          gTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAgccac                                              <  8:140541/126‑1 (MQ=255)
                                                                                                                           tCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAGCCACGCCACCACCCGCCa                                >  5:88708/1‑139 (MQ=255)
                                                                                                                           tCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAGCCACGCCACCACCCGCCa                                <  6:88708/139‑1 (MQ=255)
                                                                                                                             cAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAGCCACGCCACCACCCGCCAGGCGTTGGAg                      >  1:240877/1‑148 (MQ=255)
                                                                                                                              aaaCGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAGCCACGCCACCACCCGCCAGGCGTTGGAGc                     <  4:177166/148‑1 (MQ=255)
                                                                                                                              aaaCGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAGCCACGCCACCACCCGCCAGGCGTTGGAGCa                    >  5:307376/1‑149 (MQ=255)
                                                                                                                               aaCGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAGCCACGCCACCACCCGCCAGGCGTTGGAGCAc                   >  5:232467/1‑149 (MQ=255)
                                                                                                                               aaCGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAGCCACGCCACCACCCGCCAGGCGTTGGAGCAc                   <  2:340595/149‑1 (MQ=255)
                                                                                                                                aCGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAGCCACGCCACCACCCGCCAGGCGTTGGAGCACg                  >  8:297972/1‑149 (MQ=255)
                                                                                                                                 cgccgcgccgCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAGCCACGCCACCACCCGCCAGGCGTTGGAGCACGt                 <  6:280924/149‑1 (MQ=255)
                                                                                                                                              tCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAGCCACGCCaccacc                                     >  5:200806/1‑116 (MQ=255)
                                                                                                                                              tCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAGCCACGCCaccacc                                     <  6:200806/116‑1 (MQ=255)
                                                                                                                                                aGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAGCCACGCCACCACCCGCCAGGCGTTGGAGCACGTCGAACACATCCTCg   <  2:205470/148‑1 (MQ=255)
                                                                                                                                                aGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAGCCACGCCACCACCCGCCAGGCGTTGGAGCACGTCGAACACATCCTCGa  >  7:310047/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
ATTCGTAAATCGATTCTACATTGCGTAAAATGGTTTGTCGCGTTTTAGTTGCTCGTACGCGGCGACAAACGGCTATCCTCTCGTTTTTACTCCATGATAGTCAGATCAGCATGCCCCTTACCGTCAAACGCCGCGCCGCACTCAGCTGGCGCACGCTCTTGCCCTGGGTCGTCGGCGTGCTGCCGGTCATGTGCGGCCTGGCCGTGATGAACTGGCAGGTCGAACGTGAAATGCAGGCCAGCAGCCACGCCACCACCCGCCAGGCGTTGGAGCACGTCGAACACATCCTCGA  >  NC_002947/3890907‑3891198

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: