Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,028,994 A→G 6.8% intergenic (+71/‑368) PP_4432 → / → PP_4433 Xaa‑Pro aminopeptidase/amino acid MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,028,9940AG 6.8% 164.8 / 2.3 59intergenic (+71/‑368)PP_4432/PP_4433Xaa‑Pro aminopeptidase/amino acid MFS transporter
Reads supporting (aligned to +/- strand):  ref base A (25/30);  new base G (2/2);  total (27/32)
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.40e-01

CCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGAAGTAATACCC  >  NC_002947/5028847‑5029141
                                                                                                                                                              |                                                                                                                                                    
ccGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGACGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTaaa                                                                                                                                                      <  3:194863/148‑1 (MQ=255)
       cATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGGCTTTTCTTCAAGGACTGACGAAAAGTTTGAAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCAACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGAACCAt                                                                                                                                               >  1:248774/1‑149 (MQ=255)
         tCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCTTTAAGGCACCATCAc                                                                                                                                            <  4:189807/149‑1 (MQ=255)
         tCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAc                                                                                                                                            >  4:215439/1‑149 (MQ=255)
                   tGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACt                                                                                                                                 >  6:102999/1‑149 (MQ=255)
                   tGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACt                                                                                                                                 >  4:94445/1‑149 (MQ=255)
                         tGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTg                                                                                                                           >  5:266793/1‑149 (MQ=255)
                           gggTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGtt                                                                                                                         <  1:37459/149‑1 (MQ=255)
                           gggTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGtt                                                                                                                         <  3:94445/149‑1 (MQ=255)
                                gagaCGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTAc                                                                                                                     >  5:174132/1‑148 (MQ=255)
                                gagaCGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTAc                                                                                                                     >  4:116902/1‑148 (MQ=255)
                                  gaCGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCCTCTCCCTTGCCCCCGCCTAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCttt                                                                                                                            <  1:268680/148‑1 (MQ=255)
                                  gaCGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCa                                                                                                                                       >  8:132806/1‑128 (MQ=255)
                                  gaCGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCa                                                                                                                                       <  7:132806/128‑1 (MQ=255)
                                       tCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACaaaaaaa                                                                                                              <  3:184051/148‑1 (MQ=255)
                                        cTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATTGTGCCTTAAAACACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAg                                                                                                             <  3:76482/148‑1 (MQ=255)
                                        cTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAACACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAg                                                                                                             <  3:72085/148‑1 (MQ=255)
                                              gACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAGAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCCCGGCTATGTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACg                                                                                                         >  1:53113/1‑147 (MQ=255)
                                                    ttgCAAACTTTTCGTCAGTCCTTGACGAAAAGCTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTg                                                                                                 >  2:149420/3‑148 (MQ=38)
                                                    ttgCAAACTTTTCGTCAGTCCTTGACGAAAAGCTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTg                                                                                                 >  1:159111/3‑148 (MQ=38)
                                                      gCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTa                                                                                                                      >  7:193758/1‑125 (MQ=255)
                                                      gCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTa                                                                                                                      <  8:193758/125‑1 (MQ=255)
                                                      gCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAag                                                                                              <  3:116902/149‑1 (MQ=255)
                                                        aaGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAagag                                                                                            <  6:174132/149‑1 (MQ=255)
                                                        aaGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAagag                                                                                            >  3:284781/1‑149 (MQ=255)
                                                        aaGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAagag                                                                                            <  3:317892/149‑1 (MQ=255)
                                                           gCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAgg                                                                                           <  1:75642/149‑1 (MQ=255)
                                                                tCGTCAAGGACTGACGAAAAGGTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGGGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGaa                                                                                      <  6:236547/148‑1 (MQ=255)
                                                                     aaGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGa                                                                                <  2:290172/149‑1 (MQ=255)
                                                                       ggACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGAtt                                                                              <  2:53113/149‑1 (MQ=255)
                                                                        gACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTg                                                                             >  7:177147/1‑149 (MQ=255)
                                                                           tGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTc                                                                          >  7:258717/1‑149 (MQ=255)
                                                                              cGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTccc                                                                        <  4:284781/148‑1 (MQ=255)
                                                                                         gCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGa                                                              <  3:215439/148‑1 (MQ=255)
                                                                                           aaGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGt                                                            >  7:197873/1‑148 (MQ=255)
                                                                                             gTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCaa                                                         >  2:132986/1‑149 (MQ=255)
                                                                                              tCGCTGG‑‑‑‑‑‑‑‑‑TGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAAc                                                        >  8:34165/1‑149 (MQ=255)
                                                                                                             taGTTATATTTCCACTGCTATATTCAGCGACAATGTTATTGTGCCTTAAATCACCATCACTCCAGCCACTACTTTGTTTTACAAAAAAAGCACGCAGACTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGTAGCGAGTCAACCTGTGc                                                  <  2:248774/148‑1 (MQ=255)
                                                                                                                ttATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCt                                                >  2:214328/1‑148 (MQ=255)
                                                                                                                        ccACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTg                                        >  5:244367/1‑148 (MQ=255)
                                                                                                                        ccACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTg                                        >  3:288382/1‑148 (MQ=255)
                                                                                                                          aCGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGaaa                                     <  8:197873/149‑1 (MQ=255)
                                                                                                                          aCGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGaaa                                     <  6:248042/149‑1 (MQ=255)
                                                                                                                             gCTATCTTCAGCGACAATGTGATTGTTCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACtt                                  <  2:163234/149‑1 (MQ=255)
                                                                                                                                 tCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTAt                                <  4:288382/147‑1 (MQ=255)
                                                                                                                                 tCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTAt                                <  6:266793/147‑1 (MQ=255)
                                                                                                                                     cAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCaa                             <  1:149420/146‑1 (MQ=255)
                                                                                                                                      aGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGgac                         <  1:214328/149‑1 (MQ=255)
                                                                                                                                      aGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGga                          <  1:100318/148‑1 (MQ=255)
                                                                                                                                       gCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGgac                         >  6:257194/1‑148 (MQ=255)
                                                                                                                                                tgATGGTGCTTTAAGGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGaaa                                     <  3:150513/127‑1 (MQ=255)
                                                                                                                                                tgATGGTGCTTTAAGGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGaaa                                     >  4:150513/1‑127 (MQ=255)
                                                                                                                                                tgATGGTGCTTTAAGGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGaa                >  8:14675/1‑148 (MQ=255)
                                                                                                                                                tgATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGaa                <  8:177147/148‑1 (MQ=255)
                                                                                                                                                     gTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGaa           <  1:132986/148‑1 (MQ=255)
                                                                                                                                                     gTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGaa           <  4:169492/148‑1 (MQ=255)
                                                                                                                                                      tGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGAAg          >  7:263191/1‑148 (MQ=255)
                                                                                                                                                        ccTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACg                                            >  2:216891/1‑112 (MQ=255)
                                                                                                                                                        ccTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACg                                            <  1:216891/112‑1 (MQ=255)
                                                                                                                                                        ccTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGAAGTa        >  2:183877/1‑148 (MQ=255)
                                                                                                                                                            aaaGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGAAGTAATa     <  5:102999/147‑1 (MQ=255)
                                                                                                                                                             aaGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGCATAGAAGTAATAccc  <  5:257194/149‑1 (MQ=255)
                                                                                                                                                              |                                                                                                                                                    
CCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGAAGTAATACCC  >  NC_002947/5028847‑5029141

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: