Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,904,996 A→C 6.5% intergenic (‑163/+72) PP_4312 ← / ← PP_4313 hypothetical protein/peptidylprolyl isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,904,9960AC 6.5% 194.6 / 4.4 62intergenic (‑163/+72)PP_4312/PP_4313hypothetical protein/peptidylprolyl isomerase
Reads supporting (aligned to +/- strand):  ref base A (25/33);  new base C (2/2);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

GAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGCCCCTACGC  >  NC_002947/4904854‑4905144
                                                                                                                                              |                                                                                                                                                    
gagaATAACGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGcc                                                                                                                                                 >  6:110905/1‑148 (MQ=255)
    atagCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGt                                                                                                                                            <  7:145537/149‑1 (MQ=255)
       tCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAg                                                                                                                                         >  5:141296/1‑149 (MQ=255)
         gggTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCCTTGCCTAgggaacg                                                                                                                                        <  5:107956/148‑6 (MQ=255)
         gggTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAgg                                                                                                                                        <  7:61267/148‑1 (MQ=255)
           gTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGc                                                                                                                                     >  5:312407/1‑149 (MQ=255)
             ttACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCt                                                                                                                                    <  5:382943/148‑1 (MQ=255)
                cccGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTggg                                                                                                                                 <  2:261060/148‑1 (MQ=255)
                 ccGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGgtg                                                                                                                               >  2:8253/1‑149 (MQ=255)
                   gCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGtt                                                                                                                             <  5:110905/149‑1 (MQ=255)
                   gCCACATCGCCCGGCTAGGTGCAGGGCGCCCGTGCAGTGAATGCCAGCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGtt                                                                                                                             <  8:350651/149‑1 (MQ=255)
                     cacaTCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTgc                                                                                                                           <  5:287482/149‑1 (MQ=255)
                           gCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAAtt                                                                                                                     >  1:257955/1‑149 (MQ=255)
                           gCCCGGCTAGTGGCAGGGCGCCCGGGCAGGGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAAtt                                                                                                                     <  4:5487/149‑1 (MQ=255)
                             ccGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATtgt                                                                                                                   >  5:148078/1‑149 (MQ=255)
                             ccGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATtgt                                                                                                                   >  4:327281/1‑149 (MQ=255)
                                    gTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGt                                                                                                            <  3:182183/149‑1 (MQ=255)
                                       gCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGc                                                                                                          >  6:305566/1‑148 (MQ=255)
                                            gcgcCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGGCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAAc                                                                                                    <  2:56285/149‑1 (MQ=255)
                                                ccGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCg                                                                                                <  2:2960/149‑1 (MQ=255)
                                                 cGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCgg                                                                                               <  5:305566/149‑1 (MQ=255)
                                                   tGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGt                                                                                             <  6:148078/149‑1 (MQ=255)
                                                    gCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAAGTGTGTAAGGTGCAAGACCCTCGGGtt                                                                                            >  6:201142/1‑149 (MQ=255)
                                                      aGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTg                                                                                           >  6:219235/1‑148 (MQ=255)
                                                      aGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGc                                                                                          <  8:292899/149‑1 (MQ=255)
                                                       gTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCa                                                                                         <  6:141296/149‑1 (MQ=255)
                                                               cATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATtgt                                                                                                                   <  7:185040/115‑1 (MQ=255)
                                                               cATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATtgt                                                                                                                   >  8:185040/1‑115 (MQ=255)
                                                               cATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGaa                                                                                  >  7:240772/1‑148 (MQ=255)
                                                               cATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAac                                                                                 <  8:240772/149‑1 (MQ=255)
                                                                aTCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAaca                                                                                >  2:41800/1‑149 (MQ=255)
                                                                 tCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCggg                                                                                              <  2:5850/134‑1 (MQ=255)
                                                                 tCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCggg                                                                                              >  1:5850/1‑134 (MQ=255)
                                                                 tCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAacac                                                                               >  5:331538/1‑149 (MQ=255)
                                                                      ggCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCa                                                                            <  3:57443/147‑1 (MQ=255)
                                                                       gCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGc                                                                                                          >  7:328379/1‑116 (MQ=255)
                                                                       gCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGc                                                                                                          <  8:328379/116‑1 (MQ=255)
                                                                         aaCGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTACGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAatcat                                                                       <  2:340713/149‑1 (MQ=255)
                                                                            gcgTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATccc                                                                    <  6:312407/149‑1 (MQ=255)
                                                                                  gggTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGa                                                              <  5:219235/149‑1 (MQ=255)
                                                                                        gagTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACtt                                                        <  8:185678/149‑1 (MQ=255)
                                                                                          gTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTc                                                       >  5:143402/1‑148 (MQ=255)
                                                                                             ccccGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCcag                                                    >  7:189270/1‑148 (MQ=255)
                                                                                              cccGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCcagca                                                  <  6:331538/149‑1 (MQ=255)
                                                                                                  aaCCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTc                                              <  8:149483/149‑1 (MQ=255)
                                                                                                     ccGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGa                                            >  6:198060/1‑148 (MQ=255)
                                                                                                     ccGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGa                                            >  7:135331/1‑148 (MQ=255)
                                                                                                            gTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGaa                                     <  1:41800/148‑1 (MQ=255)
                                                                                                            gTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGaa                                     <  1:8253/148‑1 (MQ=255)
                                                                                                                ctccCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAaggtc                                 <  5:198060/148‑1 (MQ=255)
                                                                                                                   ccTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGGAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCGCGCGTGAGCACTTCCAGCAGTTCGATCTCGAaggtcagg                              <  6:204778/148‑1 (MQ=255)
                                                                                                                    cTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAaggtcaggtc                            >  3:72067/1‑149 (MQ=255)
                                                                                                                               tGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTcggcg                 <  3:327281/149‑1 (MQ=255)
                                                                                                                                 ccGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGaa                                     <  2:13031/127‑1 (MQ=255)
                                                                                                                                 ccGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGaa                                     >  1:13031/1‑127 (MQ=255)
                                                                                                                                  cctaaccTAGGCCTTGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGaa                                     >  6:308018/6‑126 (MQ=255)
                                                                                                                                  cctaaccTAGGCCTTGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGaa                                     <  5:308018/121‑1 (MQ=255)
                                                                                                                                  cctaaccTAGGCCTTGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGa              >  3:68728/6‑149 (MQ=255)
                                                                                                                                  cGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATACCGCGTGAGCACTTCCAGCAGTTCGATCTAGAAGGTAAGGTCCGAATTCGGCGGGa              >  5:43223/1‑149 (MQ=255)
                                                                                                                                    ttCCCTAGGCAAGGCGTAGGTTAGGGGGTGGGTGTTGCGCATTTGTGTATGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGGGAGCACTTCCAGCAGTTCGATGTCGAAGGTCAGGTCCGAATTCGGCGGGATc            <  8:135331/149‑1 (MQ=255)
                                                                                                                                       ccTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGcc         >  3:210705/1‑149 (MQ=255)
                                                                                                                                          aGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGcccc       >  8:298119/1‑148 (MQ=255)
                                                                                                                                              aaGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGCCCCTACGc  >  7:185335/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGCCCCTACGC  >  NC_002947/4904854‑4905144

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: ‑