Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,311,049 T→C 8.0% V267A (GTA→GCA) ‡ PP_2033 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,311,0490TC 8.0% 176.3 / 4.6 62V267A (GTA→GCA) ‡PP_2033hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/29);  new base C (2/3);  total (30/32)
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.97e-01

GGACGAGCGGCGGCATAACTTCCCGCTGGTGCGCAGTGGCCACGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTGGGCGGCGGGACGAGGCGCAAAAGCAG  >  NC_002947/2310903‑2311197
                                                                                                                                                   |                                                                                                                                                    
ggACGAGCGGCGGCATAACTTCCCGCTGGTGCGCAGTGGCCACGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAg                                                                                                                                                    <  6:204779/149‑1 (MQ=255)
    gAGCGGCGGCATAACTTCCCGCTGGTGCGCAGTGGCCACGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCAc                                                                                                                                                <  5:259529/149‑1 (MQ=255)
         gcggcATAACTTCCCGCTGGTGCGCAGTGGCCACGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTACAGGCAcgcact                                                                                                                                           <  7:97517/149‑4 (MQ=255)
                   ttCCCGCTGGTGCGCAGTGGCCACGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCAcgc                                                                                                                                              >  5:63608/1‑136 (MQ=255)
                   ttCCCGCTGGTGCGCAGTGGCCACGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCAcgc                                                                                                                                              <  6:63608/136‑1 (MQ=255)
                   ttCCCGCTGGTGCGCAGTGGCCACGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTACAGGCAcgcac                                                                                                                                            <  1:138497/138‑3 (MQ=255)
                   ttCCCGCTGGTGCGCAGTGGCCACGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTACAGGCAcgcac                                                                                                                                            >  2:138497/1‑136 (MQ=255)
                   ttCCCGCTGGTGCGCAGTGGCCACGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGCGCGTGCCTACAGGCAcgca                                                                                                                                             <  4:105679/137‑2 (MQ=255)
                   ttCCCGCTGGTGCGCAGTGGCCACGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGCGCGTGCCTACAGGCAcgca                                                                                                                                             >  3:105679/1‑136 (MQ=255)
                    tCCCGCTGGTGCGCAGTGGCCACGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGc                                                                                                                                >  3:149860/1‑149 (MQ=255)
                                      gCCACGACATCGTGCTGCCGCGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGt                                                                                                               >  7:221571/1‑148 (MQ=255)
                                        cACGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAg                                                                                                            >  4:282965/1‑149 (MQ=255)
                                        cACGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAg                                                                                                            >  5:326822/1‑149 (MQ=255)
                                         aCGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGt                                                                                                           <  1:129562/149‑1 (MQ=255)
                                           gACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTTGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGc                                                                                                         <  4:57379/149‑1 (MQ=255)
                                             cATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCt                                                                                                       <  7:42582/149‑1 (MQ=255)
                                              aTCGTGCTGCCGGGCGTTCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGTGCAGGTGTTTCTGTGCAAGCCGCAGTCCCAGTGCGTTCCTGTAGGAACGCGTGTCGATCAGTGCCGCGTACATGCCAAGGTGAGTGCCtt                                                                                                      <  5:282787/149‑1 (MQ=255)
                                                  tgctgcCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTg                                                                                                  <  5:41840/149‑1 (MQ=255)
                                                    gtgcCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGc                                                                                                <  1:271320/148‑1 (MQ=255)
                                                              tGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTTCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGttt                                                                                      <  4:149860/149‑1 (MQ=255)
                                                               gCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTg                                                                                     >  5:114732/1‑149 (MQ=255)
                                                                 aCTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGcc                                                                                   >  3:145854/1‑149 (MQ=255)
                                                                   tCGAATATTGGGGGTGGTTACCCGGAGACCCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGcc                                                                                   >  1:193381/1‑148 (MQ=255)
                                                                   tCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGttt                                                                                      >  4:93621/1‑144 (MQ=255)
                                                                   tCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGttt                                                                                      <  3:93621/144‑1 (MQ=255)
                                                                   tCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGa                                                                                 >  4:91163/1‑149 (MQ=255)
                                                                      aaTATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGatat                                                                              >  8:17918/1‑149 (MQ=255)
                                                                      aaTATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTGTGCCGatat                                                                              >  8:13308/1‑149 (MQ=255)
                                                                                   gtTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAgg                                                                 <  1:262550/149‑1 (MQ=255)
                                                                                    ttACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGc                                                                >  3:102401/1‑149 (MQ=255)
                                                                                       cccGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCcg                                                              <  2:193381/148‑1 (MQ=255)
                                                                                        ccGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCcgcg                                                            >  6:279412/1‑149 (MQ=255)
                                                                                            agaCCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGtg                                                        <  6:114732/149‑1 (MQ=255)
                                                                                             gaCCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGtga                                                       >  2:204610/1‑149 (MQ=255)
                                                                                                 ccTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAg                                                                  >  1:221833/1‑135 (MQ=255)
                                                                                                 ccTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAg                                                                  <  2:221833/135‑1 (MQ=255)
                                                                                                 ccTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGtgatga                                                    >  7:159851/1‑149 (MQ=255)
                                                                                                   ttCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAg                                                                  <  7:238008/133‑1 (MQ=255)
                                                                                                   ttCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAg                                                                  >  8:238008/1‑133 (MQ=255)
                                                                                                    tCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCt                                                 >  4:297313/1‑149 (MQ=255)
                                                                                                       ggAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTgg                                               >  3:11008/1‑148 (MQ=255)
                                                                                                             ggTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGc                                        >  2:170276/1‑149 (MQ=255)
                                                                                                                    cTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGAtt                                   >  7:19499/1‑147 (MQ=255)
                                                                                                                    cTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGAtt                                   <  8:19499/147‑1 (MQ=255)
                                                                                                                       tgCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGcc                                <  4:145854/147‑1 (MQ=255)
                                                                                                                       tgCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCg                               <  6:326822/148‑1 (MQ=255)
                                                                                                                          aaGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTgg                           >  2:300703/1‑149 (MQ=255)
                                                                                                                          aaGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTgg                           <  1:300703/149‑1 (MQ=255)
                                                                                                                            gCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTGggc                         <  8:159851/149‑1 (MQ=255)
                                                                                                                              cGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTGggcgg                       >  8:88538/1‑149 (MQ=255)
                                                                                                                               gCAGTCCCAGTGCATGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGc                                 >  7:238654/1‑138 (MQ=255)
                                                                                                                               gCAGTCCCAGTGCATGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGc                                 <  8:238654/138‑1 (MQ=255)
                                                                                                                                cAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCg                               >  6:82522/1‑139 (MQ=255)
                                                                                                                                cAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCg                               <  5:82522/139‑1 (MQ=255)
                                                                                                                                cAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTGGgcggcg                     <  3:282965/149‑1 (MQ=255)
                                                                                                                                   tCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTGGGCGGCGGGa                  <  7:13308/149‑1 (MQ=255)
                                                                                                                                    cccAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTGGGCGGCGGGa                  <  7:17918/148‑1 (MQ=255)
                                                                                                                                         tgcgtgcCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTGGGCGGCGGGACGAGgc            <  4:102401/149‑1 (MQ=255)
                                                                                                                                              gcCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGGCCTGGGAGGAGCCGAtt                                   >  3:44329/1‑121 (MQ=255)
                                                                                                                                              gcCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGGCCTGGGAGGAGCCGAtt                                   >  3:48428/1‑121 (MQ=255)
                                                                                                                                              gcCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGGCCTGGGAGGAGCCGAtt                                   <  4:48428/121‑1 (MQ=255)
                                                                                                                                              gcCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGGCCTGGGAGGAGCCGAtt                                   <  4:44329/121‑1 (MQ=255)
                                                                                                                                                   tAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGTCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTGGGCGGCGTGACGAGGCGCAAAAGCAg  <  8:221571/149‑1 (MQ=255)
                                                                                                                                                   tAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTGGGCGGCGGGACGAGGCGCAAAAGc    <  3:91163/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGACGAGCGGCGGCATAACTTCCCGCTGGTGCGCAGTGGCCACGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTGGGCGGCGGGACGAGGCGCAAAAGCAG  >  NC_002947/2310903‑2311197

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: