Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,741,508 A→T 7.6% intergenic (‑179/+49) PP_5038 ← / ← PP_5039 lipoprotein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,741,5080AT 7.6% 193.6 / 5.3 66intergenic (‑179/+49)PP_5038/PP_5039lipoprotein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/28);  new base T (3/2);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

AGTGGCCTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGAAGCGGGC  >  NC_002947/5741361‑5741655
                                                                                                                                                   |                                                                                                                                                   
aGTGGCCTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAAt                                                                                                                                                    <  3:41201/149‑1 (MQ=255)
         cAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCGCTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCg                                                                                                                                             <  8:40432/147‑1 (MQ=255)
         cAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGc                                                                                                                                           <  7:242185/149‑1 (MQ=255)
          aGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCAATGTACCTACACCAATGGTCGGATGGCCAATACTGTGCGATTAATCACAAAGCAATAATCATGCGACCAATAAAAAAGGCCGCTCCCATTGGGGGCGGcc                                                                                                                                          >  7:208770/1‑149 (MQ=255)
            tGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCtt                                                                                                                                        <  4:2260/149‑1 (MQ=255)
             gACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGTTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTg                                                                                                                                       >  6:174785/1‑149 (MQ=255)
               ccTTCGACGGGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCTCCCATTGGGGGCGGCCttttt                                                                                                                                     <  8:208770/149‑4 (MQ=255)
                      cGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGAt                                                                                                                              <  3:268000/149‑1 (MQ=255)
                            aaTCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGc                                                                                                                        <  1:117821/149‑1 (MQ=255)
                                  tcccAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCt                                                                                                                          >  7:6528/1‑141 (MQ=255)
                                  tcccAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCt                                                                                                                          <  8:6528/141‑1 (MQ=255)
                                  tcccAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCAcag                                                                                                                    >  7:76546/1‑147 (MQ=255)
                                        gTTCCATGCCACACCCCCTATGTCCCAACACCATTGGTCGGATGGCCACTACTGTGCGAATAATCACAAAGCATTAAACATGCGCCCAATAAAAAAGACCGCCCCCAATGGGAGCGGCATAGCGTCGCGATGCATGAACAGCAGCcccc                                                                                                            >  2:272705/1‑149 (MQ=255)
                                         ttCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCAGAACGGATTAATCATGCGCCGAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGccccg                                                                                                           <  7:207337/149‑2 (MQ=255)
                                             cTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTggg                                                                                                       >  5:224021/1‑149 (MQ=255)
                                             cTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTggg                                                                                                       >  1:267739/1‑149 (MQ=255)
                                                ccccACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGc                                                                                                      <  4:209956/147‑1 (MQ=255)
                                                  ccACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCt                                                                                                    >  7:171318/1‑147 (MQ=255)
                                                       cccTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAAcc                                                                                             >  4:304104/1‑149 (MQ=255)
                                                           aTGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGa                                                                                         >  3:264779/1‑149 (MQ=255)
                                                                ccTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAg                                                                                    <  8:171318/149‑1 (MQ=255)
                                                                 cTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAgg                                                                                   >  3:167413/1‑149 (MQ=255)
                                                                       cATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTggg                                                                                                       <  7:55388/123‑1 (MQ=255)
                                                                       cATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTggg                                                                                                       >  8:55388/1‑123 (MQ=255)
                                                                        aTTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAg                                                                             >  1:127034/1‑148 (MQ=255)
                                                                            gTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCaaa                                                                         >  8:239712/1‑148 (MQ=255)
                                                                             tctgATGGCCACTACTGTGCGAATAATCACAAAGCATTAATCATGCGCCCATTAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCaaaaa                                                                       <  1:64205/146‑1 (MQ=255)
                                                                             tCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGc                                                                            <  2:160064/144‑1 (MQ=255)
                                                                             tCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGc                                                                            >  1:160064/1‑144 (MQ=255)
                                                                             tCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCaaaaa                                                                       <  8:76546/149‑1 (MQ=255)
                                                                                 aTGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTc                                                                    >  1:203661/1‑148 (MQ=255)
                                                                                     ccACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCAtct                                                               <  4:167413/149‑1 (MQ=255)
                                                                                      cACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCAtctc                                                              <  3:304104/149‑1 (MQ=255)
                                                                                         tactGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTc                                                           >  2:41489/1‑149 (MQ=255)
                                                                                         tactGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTc                                                           >  1:215975/1‑149 (MQ=255)
                                                                                         tactGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTc                                                           <  2:267739/149‑1 (MQ=255)
                                                                                           ctGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCACCTTGGGCCTCAGAACCGGGACCCAGGCTACAGCAAAAAGTCCAACTCTTCGc                                                         >  8:258957/1‑149 (MQ=255)
                                                                                            tgtgCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCt                                                        >  5:329001/1‑149 (MQ=255)
                                                                                                      aTCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCgg                                               <  6:224021/148‑1 (MQ=255)
                                                                                                      aTCACAAAGCAATAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGGACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCgg                                               <  8:259216/148‑1 (MQ=255)
                                                                                                        cacaAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCg                                            >  7:156066/1‑149 (MQ=255)
                                                                                                         acaAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCg                                            >  3:219687/1‑148 (MQ=255)
                                                                                                         acaAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCg                                            >  4:181525/1‑148 (MQ=255)
                                                                                                          caAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCt                                          <  6:329001/149‑1 (MQ=255)
                                                                                                           aaaGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTc                                         >  7:304168/1‑149 (MQ=255)
                                                                                                               cATTAATCATGCGCCCAATAAAACAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGTTGGTGCTCGGCCGCTCCAGc                                     <  8:270322/149‑1 (MQ=255)
                                                                                                               cATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTc                                                                    >  7:85530/1‑118 (MQ=255)
                                                                                                               cATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTc                                                                    <  8:85530/118‑1 (MQ=255)
                                                                                                               cATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGc                                     >  5:73131/1‑149 (MQ=255)
                                                                                                                   aaTCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGc                                                   >  2:294781/1‑131 (MQ=255)
                                                                                                                   aaTCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGc                                                   <  1:294781/131‑1 (MQ=255)
                                                                                                                    aTCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCa                                 <  4:143385/148‑1 (MQ=255)
                                                                                                                    aTCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAg                                >  5:125590/1‑149 (MQ=255)
                                                                                                                       aaccccccAAAAAAAAACCCCCCCCCCAAACCGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGtt                             <  2:127034/144‑1 (MQ=255)
                                                                                                                       aTGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGtt                             >  5:307513/1‑149 (MQ=255)
                                                                                                                         gcgcCCAATAAAAAAGGCCGCCCCCAATGGGGGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAACAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTtgc                           >  1:36485/1‑149 (MQ=255)
                                                                                                                          cgcCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTtgcg                          <  6:35187/149‑1 (MQ=255)
                                                                                                                               aaTAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGACCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTtgcggtgcg                     >  7:58349/1‑149 (MQ=255)
                                                                                                                                  aaaaaaGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCggg                   >  1:146692/1‑148 (MQ=255)
                                                                                                                                    gcaaGGCCGCTCCCATTGGGGGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAg                >  1:236561/3‑149 (MQ=255)
                                                                                                                                    aaaaGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCATAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAg                >  5:274070/1‑149 (MQ=255)
                                                                                                                                     aaaGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAgc               >  3:281202/1‑149 (MQ=255)
                                                                                                                                       aGGCCGCCCCCAATGGGGGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAACAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAgcac             <  2:36485/149‑3 (MQ=255)
                                                                                                                                         gccgctCCCATTGGGGGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAgcgc             >  6:291153/1‑147 (MQ=255)
                                                                                                                                         gccgctCCCATTGGGGGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAgcgc             <  5:291153/147‑1 (MQ=255)
                                                                                                                                         gccgccCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGccc           <  5:174785/149‑1 (MQ=255)
                                                                                                                                                cccaTGGGCCCGGCCCTGCGGCCCCCTGCCCGCACACCCCGCCCCCTGGGCCCCAGCACCCGGACCCCCGCCCCACCATAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGAAGCgg    <  7:239712/146‑1 (MQ=255)
                                                                                                                                                   aTGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGCACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGAAGCGGGc  >  2:233720/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGTGGCCTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGAAGCGGGC  >  NC_002947/5741361‑5741655

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: