Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,902,949 G→C 9.5% L30V (CTG→GTG)  PP_5175 ← HlyD family secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,902,9490GC 9.5% 210.1 / 9.5 72L30V (CTG→GTG) PP_5175HlyD family secretion protein
Reads supporting (aligned to +/- strand):  ref base G (30/35);  new base C (5/2);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

TGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGG  >  NC_002947/5902801‑5903096
                                                                                                                                                    |                                                                                                                                                   
tGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAg                                                                                                                                                     <  5:37568/149‑1 (MQ=255)
 gATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAgt                                                                                                                                                    >  4:303834/1‑149 (MQ=255)
   tCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAgtgt                                                                                                                                                  <  5:67963/149‑1 (MQ=255)
      cccACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGtctc                                                                                                                                               <  6:87186/149‑1 (MQ=255)
      ccaaCGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCACCGAAATCGGTGGGTTGTACCTGCAGCACACCCACCCGAGGCACTGtctc                                                                                                                                               >  3:304456/4‑149 (MQ=255)
       ccACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCg                                                                                                                                              <  3:303834/149‑1 (MQ=255)
       ccACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCg                                                                                                                                              >  7:119485/1‑149 (MQ=255)
              gACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGACACTGCCTCgggtg                                                                                                                                          >  4:200990/1‑144 (MQ=255)
              gACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGACACTGCCTCgggtg                                                                                                                                          <  3:200990/146‑3 (MQ=255)
               aCGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAATTCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACctc                                                                                                                                      <  1:258183/149‑1 (MQ=255)
                 gCTACTCGAAATGTTCTTGCCGCCGACGCTGAATGACAAATCGGTCTGCACGCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGGGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACctct                                                                                                                                     <  7:154814/148‑1 (MQ=255)
                   tACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTcc                                                                                                                                  >  7:91567/1‑149 (MQ=255)
                    aCGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCt                                                                                                                                 >  4:124074/1‑149 (MQ=255)
                    aCGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCt                                                                                                                                 >  6:212195/1‑149 (MQ=255)
                      gcgAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTg                                                                                                                               <  4:320158/149‑1 (MQ=255)
                      gcgAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGACACTGCCTCgggt                                                                                                                                           <  1:37187/137‑2 (MQ=255)
                      gcgAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGACACTGCCTCgggt                                                                                                                                           >  2:37187/1‑136 (MQ=255)
                       cgAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGc                                                                                                                              <  6:1031/149‑1 (MQ=255)
                       cgAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGc                                                                                                                              <  5:212195/149‑1 (MQ=255)
                         aaaTGATCTTGCCGCCGACGCGGTATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCgggg                                                                                                                                           >  1:77557/1‑134 (MQ=255)
                         aaaTGATCTTGCCGCCGACGCGGTATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCgggg                                                                                                                                           <  2:77557/134‑1 (MQ=255)
                          aaTGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTa                                                                                                                            >  5:307444/1‑148 (MQ=255)
                               tCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCg                                                                                                                      >  7:191579/1‑149 (MQ=255)
                                  tGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGt                                                                                                                    >  2:217296/1‑148 (MQ=255)
                                   gccgccGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTg                                                                                                                   <  1:294167/148‑1 (MQ=255)
                                    ccgccgACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGa                                                                                                                  <  6:129641/148‑1 (MQ=255)
                                    ccgccgACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAg                                                                                                                 <  2:96811/149‑1 (MQ=255)
                                    ccgccgACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAg                                                                                                                 <  2:3131/149‑1 (MQ=255)
                                            gcgGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTac                                                                                                          >  6:59904/1‑148 (MQ=255)
                                              ggAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTacaa                                                                                                        <  7:37349/148‑1 (MQ=255)
                                              ggAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTacaa                                                                                                        >  8:37349/1‑148 (MQ=255)
                                               gAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTacaaca                                                                                                      <  5:40364/149‑1 (MQ=255)
                                                 aTGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGc                                                                                                    <  8:91567/149‑1 (MQ=255)
                                                 aTGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGc                                                                                                    <  2:17200/149‑1 (MQ=255)
                                                     cAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCa                                                                                                <  8:191579/149‑1 (MQ=255)
                                                      aaaTCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACCCCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAg                                                                                               >  7:239555/1‑149 (MQ=255)
                                                      aaaTCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTAc                                                                                                                           >  5:260539/1‑121 (MQ=255)
                                                      aaaTCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTAc                                                                                                                           <  6:260539/121‑1 (MQ=255)
                                                       aaTCGGTCTGCACCCGCGCCTGCACGTCGCCGGACAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAgg                                                                                              <  7:131651/149‑1 (MQ=255)
                                                         tCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAg                                                                                                                 <  1:132348/128‑1 (MQ=255)
                                                         tCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAg                                                                                                                 >  2:132348/1‑128 (MQ=255)
                                                                 cACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAAc                                                                                      <  6:307444/147‑1 (MQ=255)
                                                                    ccGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGAACTCTTCCTTGCTACAGCCGCTGAGCATGTACAACATCCCCAGGCAGACAACCGACa                                                                                 >  7:292332/1‑149 (MQ=255)
                                                                     cgcgCCTGCACGTGGCCTGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCAGACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAg                                                                                <  4:304456/149‑1 (MQ=255)
                                                                        gcCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGa                                                                                          <  1:211054/136‑1 (MQ=255)
                                                                        gcCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGa                                                                                          >  2:211054/1‑136 (MQ=255)
                                                                              aCGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCAt                                                                       >  4:106940/1‑149 (MQ=255)
                                                                                gTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATc                                                                      >  7:215699/1‑148 (MQ=255)
                                                                                        tCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATg                                                             >  6:33245/1‑149 (MQ=255)
                                                                                            aGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTgg                                                         <  3:124074/149‑1 (MQ=255)
                                                                                             gTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGc                                                        <  5:59904/149‑1 (MQ=255)
                                                                                             gTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGc                                                        <  8:239555/149‑1 (MQ=255)
                                                                                                aCCCTGGCTGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAgg                                                                                              <  7:114782/108‑1 (MQ=255)
                                                                                                aCCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAgg                                                                                              >  8:114782/1‑108 (MQ=255)
                                                                                                aCCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGCTc                                                      <  6:17350/148‑1 (MQ=255)
                                                                                                  ccTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGCTCctt                                                   >  5:142332/1‑149 (MQ=255)
                                                                                                  ccTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGCTCctt                                                   <  5:33245/149‑1 (MQ=255)
                                                                                                        aGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGCTCCTTGCTtgt                                             >  5:4735/1‑149 (MQ=255)
                                                                                                          cGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGCTCCTTGCTtgtg                                            >  6:222507/1‑148 (MQ=255)
                                                                                                                     ggTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGAttt                                <  8:119485/149‑1 (MQ=255)
                                                                                                                     ggTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGATACAGCCGGTGAGCATGTACAAAAGACCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGATGGCTCCTTGCTTGTGCACAAATGAttt                                >  7:215707/1‑149 (MQ=255)
                                                                                                                         gTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTtgtg                             <  3:106940/148‑1 (MQ=255)
                                                                                                                          taactgcttCACACACACACGAGTCTGTGTCTCGGTGACTTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTTTCTCCTTGCTTTTGCACAAATGATTtgttg                            <  8:292332/145‑3 (MQ=255)
                                                                                                                               gctgcACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGc                       >  1:173631/1‑148 (MQ=255)
                                                                                                                                ctgcACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCGTTGCTACAGCCGGTGAGCATGTACAACAGGCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCg                      <  1:217296/148‑1 (MQ=255)
                                                                                                                                   cacaCCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACtg                  >  6:332018/1‑149 (MQ=255)
                                                                                                                                      acccacccGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGa               <  5:332018/149‑1 (MQ=255)
                                                                                                                                       ggtccccGAGACACTGCCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAg              >  4:17835/5‑149 (MQ=255)
                                                                                                                                        ccacccGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGc             >  7:202981/1‑149 (MQ=255)
                                                                                                                                        ccacccGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGc             <  8:215699/149‑1 (MQ=255)
                                                                                                                                        ccacccGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGc             >  5:63241/1‑149 (MQ=255)
                                                                                                                                         tccccGAGACACTGCCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCg            >  1:333480/3‑149 (MQ=255)
                                                                                                                                          acccGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGt           >  7:69922/1‑149 (MQ=255)
                                                                                                                                           cccGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTa          >  7:287281/1‑149 (MQ=255)
                                                                                                                                                   aGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAggg  >  3:111607/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCTGCACGTCGCCGGTCAAAGTGACCCTGGCAGCGAAATCGGTGGGTTGTACCTGCTGCACACCCACCCGAGGCAGTGTCTCGGGGACCTCTTCCTTGCTACAGCCGGTGAGCATGTACAACAGCCCCAGGCAGACAACCGACAGCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGG  >  NC_002947/5902801‑5903096

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