Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,042,321 C→A 7.4% L137I (CTC→ATC) ‡ PP_3563 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,042,3210CA 7.4% 286.0 / 8.7 89L137I (CTC→ATC) ‡PP_3563hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (39/43);  new base A (4/3);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

AGCCTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCGGG  >  NC_002947/4042177‑4042464
                                                                                                                                                 |                                                                                                                                               
aGCCTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCa                                                                                                                                              <  2:12143/148‑1 (MQ=255)
  ccTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCg                                                                                                                                           >  2:329459/1‑149 (MQ=255)
   cTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCgg                                                                                                                                          >  7:279030/1‑149 (MQ=255)
      tCGGGGCCACGCGACGGACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGTTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGgcc                                                                                                                                        <  8:262542/148‑1 (MQ=255)
      tCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGgccg                                                                                                                                       >  2:58324/1‑149 (MQ=255)
      tCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGgc                                                                                                                                         <  6:31888/147‑1 (MQ=255)
       cGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGgccgc                                                                                                                                      >  2:274026/1‑149 (MQ=255)
        ggggCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGAGCAGCGgccgt                                                                                                                                      <  5:401648/148‑2 (MQ=255)
            ccACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTggg                                                                                                                                 <  1:29089/149‑1 (MQ=255)
            ccACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTggg                                                                                                                                 <  1:317236/149‑1 (MQ=255)
                 cgACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGTATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACt                                                                                                                            <  1:78064/149‑1 (MQ=255)
                     gcacgcGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCa                                                                                                                        <  7:331838/149‑1 (MQ=255)
                        cgcgCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAgg                                                                                                                      >  3:137595/1‑148 (MQ=255)
                        cgcgCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAgg                                                                                                                      <  4:137595/148‑1 (MQ=255)
                         gcgcCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGcc                                                                                                                    >  1:288639/1‑149 (MQ=255)
                           gccgccATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGa                                                                                                                  <  2:288639/149‑1 (MQ=255)
                            ccgccATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAg                                                                                                                 <  7:290646/149‑1 (MQ=255)
                            ccgccATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAg                                                                                                                 >  1:178214/1‑149 (MQ=255)
                               ccATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCGGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGa                                                                                                              >  7:106341/1‑149 (MQ=255)
                                 atgcatgcCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACa                                                                                                            <  8:378559/149‑1 (MQ=255)
                                   gcatgcCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAg                                                                                                           >  7:24004/1‑148 (MQ=255)
                                     atgcCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGGGGCCGCCTGGGCGAGTGGCAGGCCGAGGGACAGAAg                                                                                                        <  5:91867/149‑1 (MQ=255)
                                     atgcCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAg                                                                                                        <  2:178214/149‑1 (MQ=255)
                                         cATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGtgc                                                                                                    >  3:53196/1‑149 (MQ=255)
                                          aTCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGtgct                                                                                                   >  6:56141/1‑149 (MQ=255)
                                           tCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGtgctg                                                                                                  >  1:16118/1‑149 (MQ=255)
                                             gAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCt                                                                                                                                    <  3:198746/113‑1 (MQ=255)
                                             gAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCt                                                                                                                                    >  4:198746/1‑113 (MQ=255)
                                              aGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAgg                                                                                                                <  5:136224/132‑1 (MQ=255)
                                              aGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAgg                                                                                                                >  6:136224/1‑132 (MQ=255)
                                                 aaaGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCaggag                                                                                            >  5:54931/1‑149 (MQ=255)
                                                  aaGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCaggagg                                                                                           <  2:346258/149‑1 (MQ=255)
                                                   aGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGt                                                                                          >  3:26582/1‑149 (MQ=255)
                                                    gCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGt                                                                                          >  4:304640/1‑148 (MQ=255)
                                                    gCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTc                                                                                         <  1:58324/149‑1 (MQ=255)
                                                         tGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAAcgc                                                                                    >  4:61513/1‑149 (MQ=255)
                                                         tGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAAcgc                                                                                    >  4:303234/1‑149 (MQ=255)
                                                                  gtgCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAg                                                                                                                                             <  6:259679/83‑1 (MQ=255)
                                                                  gtgCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAg                                                                                                                                             >  5:259679/1‑83 (MQ=255)
                                                                    gCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAg                                                                          <  1:329459/149‑1 (MQ=255)
                                                                    gCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAgg                                                                         >  5:72078/1‑149 (MQ=255)
                                                                     cTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGt                                                                        <  5:56141/149‑1 (MQ=255)
                                                                     cTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGATCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGt                                                                        >  2:38414/1‑149 (MQ=255)
                                                                      tCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTg                                                                       >  3:259724/1‑149 (MQ=255)
                                                                        gACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGc                                                                      >  1:155134/1‑148 (MQ=255)
                                                                          cccGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCAc                                                                   <  8:24004/149‑1 (MQ=255)
                                                                           cggCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACg                                                                  <  8:106341/147‑1 (MQ=255)
                                                                            cGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTg                                                                                                                           >  1:397482/1‑91 (MQ=255)
                                                                            cGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTg                                                                                                                           <  2:397482/91‑1 (MQ=255)
                                                                            cGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACgg                                                                 >  1:188539/1‑149 (MQ=255)
                                                                            cGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACgg                                                                 <  8:279030/149‑1 (MQ=255)
                                                                            cGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACgg                                                                 <  1:353828/149‑1 (MQ=255)
                                                                                aTTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAAt                                                             >  5:355813/1‑149 (MQ=255)
                                                                                  tcgtcgATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATgc                                                           >  3:8693/1‑149 (MQ=255)
                                                                                    gtcgATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATgcg                                                          <  4:53196/148‑1 (MQ=255)
                                                                                    gtcgATATCGGCGGCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATgcgg                                                         <  7:381807/149‑1 (MQ=255)
                                                                                        atatCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCggcggc                                                     <  1:208614/149‑1 (MQ=255)
                                                                                           tCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGTCAATGCGGCGGCCAg                                                  >  8:67204/1‑149 (MQ=255)
                                                                                           tCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAg                                                  <  3:61513/149‑1 (MQ=255)
                                                                                           tCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAg                                                  <  3:303234/149‑1 (MQ=255)
                                                                                              gCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAAcgcc                                                                                   >  4:250464/1‑113 (MQ=255)
                                                                                              gCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAAcgcc                                                                                   <  3:250464/113‑1 (MQ=255)
                                                                                              gCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCg                                               >  1:376567/1‑149 (MQ=255)
                                                                                                 tCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGttt                                            <  3:304640/149‑1 (MQ=255)
                                                                                                     ccAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTgccg                                        <  4:8693/149‑1 (MQ=255)
                                                                                                      cAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTgccgc                                       <  4:259724/149‑1 (MQ=255)
                                                                                                            ggCCGCGACTGGCGCATCGTGATGGCACCGCCGCTGCTCAGCGGACCACAGGCCGACACGCACGCCGAGTGACAGAAGGTCCTGCACGAGGTCAACGCCGCGCGCAAGCTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGccgcccc                                 >  4:202938/1‑149 (MQ=255)
                                                                                                               cgcgACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGccgccccgc                               >  5:181894/1‑148 (MQ=255)
                                                                                                               cgcgACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGccgccccgc                               >  3:394483/1‑148 (MQ=255)
                                                                                                                gcgACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCa                             <  2:376567/149‑1 (MQ=255)
                                                                                                                    cTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACt                           >  4:108197/1‑147 (MQ=255)
                                                                                                                       gcgcATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGAGGCCCGCCGTTTGCCGCCGCCCCGCCACTAAGCt                      >  1:327454/1‑149 (MQ=255)
                                                                                                                        cgcATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATgcgg                                                         <  2:400576/113‑1 (MQ=255)
                                                                                                                        cgcATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTg                     <  4:386432/149‑1 (MQ=255)
                                                                                                                          cATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGa                   <  1:274026/149‑1 (MQ=255)
                                                                                                                           aTCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAg                  <  2:155134/149‑1 (MQ=255)
                                                                                                                               tGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGca                >  1:233761/1‑147 (MQ=255)
                                                                                                                               tGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGca                >  6:216843/1‑147 (MQ=255)
                                                                                                                                   ggCGCGGCCGCTGCTCAGCGGCCGCCTGGGCCACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCCCCCCGCCACAAAGCTGGAGCACAGTGc          >  4:141155/1‑149 (MQ=255)
                                                                                                                                    aggcGGCCGCTGAGCAGCGGCCGGCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGtt                                             <  7:307256/110‑1 (MQ=255)
                                                                                                                                    aggcGGCCGCTGAGCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGtt                                             >  8:307256/4‑113 (MQ=255)
                                                                                                                                      aCGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTg        >  1:315703/1‑148 (MQ=255)
                                                                                                                                      aCGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTg        <  2:315703/148‑1 (MQ=255)
                                                                                                                                        ggCCGCTGAGCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGcaca              <  7:26596/140‑1 (MQ=255)
                                                                                                                                        ggCCGCTGAGCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGcc     >  6:175717/1‑149 (MQ=255)
                                                                                                                                        ggCCGCTGAGCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCAGTTTGCCGCCGCCCCGCCACTAAGCTGGAGcaca              >  8:26596/1‑140 (MQ=255)
                                                                                                                                         gCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCg    <  6:54931/149‑1 (MQ=255)
                                                                                                                                            cGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCggg  <  2:233761/149‑1 (MQ=255)
                                                                                                                                            cGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGGGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCAGCCGCCCCTCCACTAAGCTGGAGCACAGTTCTGGCCggg  <  5:40763/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
AGCCTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCGGG  >  NC_002947/4042177‑4042464

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: