Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 337,825 C→G 17.0% intergenic (‑72/+57) PP_0279 ← / ← PP_0280 hypothetical protein/amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947337,8250CG17.0% 182.2 / 25.8 66intergenic (‑72/+57)PP_0279/PP_0280hypothetical protein/amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (24/30);  new base G (5/6);  total (29/37)
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.45e-01

GCTGGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACA  >  NC_002947/337678‑337970
                                                                                                                                                    |                                                                                                                                                 
gctgGTGCTCTATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCaa                                                                                                                                                 <  1:1447144/150‑1 (MQ=255)
gctgGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCaa                                                                                                                                                 <  1:1734344/150‑1 (MQ=255)
 ctgGTGCTTGTTCTGGAACTGGTTAAGGGCATGGAAGCGGGTGCTTAGGTGGGTGCTGTGGATGTGGTTGGGCATGGCAAGGAGGTGTGGTGGTTTGGCGGGCCTATTCGCGGGCACGCCCGCTGTCACAAGTATAGCGCCACATTTAAg                                                                                                                                                <  2:121999/150‑1 (MQ=255)
       gtCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGgc                                                                                                                                          <  1:1497189/149‑1 (MQ=255)
       cTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGGGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATTGCAAGGATGCGTGGTGGTTTTGAGGGCCTCTTCGCGAGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGgc                                                                                                                                          <  1:587237/150‑1 (MQ=255)
           aTCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCGGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTa                                                                                                                                      <  1:750875/150‑1 (MQ=255)
            tCAGGAATTGGTCGAGGGCTTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGTGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGTGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGTCACATTCAAGGTTGGCGCTac                                                                                                                                     <  2:891904/150‑1 (MQ=255)
              aGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACCTGTATAGCGCCACATTCAAGGTTGGCGCTacac                                                                                                                                   <  1:912829/150‑1 (MQ=255)
              aGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTacac                                                                                                                                   <  1:1641595/150‑1 (MQ=255)
                    tGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTgg                                                                                                                             <  1:1437200/150‑1 (MQ=255)
                    tGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTgg                                                                                                                             <  1:2091677/150‑1 (MQ=255)
                    tGGGCGAGGGCCTGGAATTGGGTGCCTGGCTTGGGGCTTTGGAGGTGACTGGGCATGGCAAGGACGCGTTGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTgg                                                                                                                             <  1:686196/150‑1 (MQ=255)
                      gTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGa                                                                                                                           <  1:1299517/150‑1 (MQ=21)
                        cGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCg                                                                                                                        <  1:2020815/150‑1 (MQ=21)
                          aGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCgg                                                                                                                       <  1:2315872/150‑1 (MQ=255)
                              ccTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGt                                                                                                                   <  1:1134793/150‑1 (MQ=21)
                                tGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGt                                                                                                                   >  1:202733/1‑148 (MQ=255)
                                tGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGt                                                                                                                   <  2:202733/148‑1 (MQ=255)
                                             ttAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGc                                                                                                                         >  2:595318/1‑129 (MQ=255)
                                             ttAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGc                                                                                                                         <  1:595318/129‑1 (MQ=255)
                                             ttAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCggg                                                                                                                      <  1:692831/132‑1 (MQ=255)
                                             ttAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCggg                                                                                                                      >  2:692831/1‑132 (MQ=255)
                                                         tgtgGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCAACCTTGAATGTGGCGCTATACCtgtg                                                                                                                              <  2:2178514/112‑1 (MQ=18)
                                                         tgtgGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCAACCTTGAATGTGGCGCTATACCtgtg                                                                                                                              >  1:2178514/1‑112 (MQ=18)
                                                                gtggtTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCt                                                                                 >  1:1980798/1‑150 (MQ=255)
                                                                       ggCATGGCAAGGATGTGTGGTGGTTTTGATGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTggcagg                                                                          <  2:398727/150‑1 (MQ=255)
                                                                                              ttttGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTcc                                                   >  1:1863579/1‑150 (MQ=255)
                                                                                                                gggCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAg                                       >  2:435847/1‑144 (MQ=18)
                                                                                                                gggCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAg                                       <  1:435847/144‑1 (MQ=18)
                                                                                                                   cACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGaa                                        >  2:2204622/1‑140 (MQ=255)
                                                                                                                   cACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGaa                                        <  1:2204622/140‑1 (MQ=255)
                                                                                                                    aCGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTggcag                                                                           <  1:2235208/104‑1 (MQ=255)
                                                                                                                    aCGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTggcag                                                                           >  2:2235208/1‑104 (MQ=255)
                                                                                                                      gCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTg                                  >  1:2072102/1‑143 (MQ=18)
                                                                                                                      gCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTg                                  <  2:2072102/143‑1 (MQ=18)
                                                                                                                        ccGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCtt                                   <  1:410643/140‑1 (MQ=255)
                                                                                                                        ccGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCtt                                   >  2:410643/1‑140 (MQ=255)
                                                                                                                          gCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTggca                                                                            >  1:426542/1‑97 (MQ=255)
                                                                                                                          gCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTggca                                                                            <  2:426542/97‑1 (MQ=255)
                                                                                                                          gCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGa                       >  1:1707662/1‑150 (MQ=255)
                                                                                                                            tCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                     <  2:1980798/150‑1 (MQ=255)
                                                                                                                              ccACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACga                   >  2:372534/1‑150 (MQ=21)
                                                                                                                              ccACAGGTGTAGCGCCAACCTTGAATGTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACga                   >  2:1135081/1‑150 (MQ=25)
                                                                                                                              ccACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACga                   >  2:215973/1‑150 (MQ=255)
                                                                                                                               cacaGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCGGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACgag                  >  2:1649056/1‑150 (MQ=255)
                                                                                                                                acaGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGt                                                                                            >  2:1682117/1‑75 (MQ=255)
                                                                                                                                acaGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGt                                                                                            <  1:1682117/75‑1 (MQ=255)
                                                                                                                                acaGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCtgtg                                         <  2:911441/126‑4 (MQ=255)
                                                                                                                                acaGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCtgtg                                         >  1:911441/1‑123 (MQ=255)
                                                                                                                                acaGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACgaga                 >  1:272122/1‑150 (MQ=255)
                                                                                                                                  aGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTg                                  >  2:1638688/1‑131 (MQ=255)
                                                                                                                                  aGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTg                                  <  1:1638688/131‑1 (MQ=255)
                                                                                                                                  aGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATc               >  2:290422/1‑150 (MQ=255)
                                                                                                                                    gTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCAc             >  1:1146951/1‑150 (MQ=255)
                                                                                                                                      ataGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGcag                                                                                                      >  2:154896/1‑59 (MQ=255)
                                                                                                                                      ataGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGcag                                                                                                      <  1:154896/59‑1 (MQ=255)
                                                                                                                                       taGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCa                                                          <  1:1197706/102‑1 (MQ=255)
                                                                                                                                       taGCGCCACATTCAAGGATGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCa                                                          >  2:1197706/1‑102 (MQ=255)
                                                                                                                                            ccACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCa                                                          >  2:2097353/1‑97 (MQ=255)
                                                                                                                                            ccACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCa                                                          <  1:2097353/97‑1 (MQ=255)
                                                                                                                                            ccACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTAAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                     >  1:1963530/1‑134 (MQ=255)
                                                                                                                                            ccACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTAAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                     <  2:1963530/134‑1 (MQ=255)
                                                                                                                                             cacaTTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGTCAGCGAAGCATGCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGaa    >  1:408169/1‑150 (MQ=255)
                                                                                                                                                caTTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAaca  <  1:215973/150‑1 (MQ=255)
                                                                                                                                                  ttCAAGGTTGGCGCTACACCTGTGGGAGCGGGAGTGCCCGCGAAGCAGTCAGTGa                                                                                               >  1:403178/1‑55 (MQ=255)
                                                                                                                                                  ttCAAGGTTGGCGCTACACCTGTGGGAGCGGGAGTGCCCGCGAAGCAGTCAGTGa                                                                                               <  2:403178/55‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GCTGGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACA  >  NC_002947/337678‑337970

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