Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,029,042 A→G 7.0% 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,029,0420AG 7.0% 246.2 / 6.8 85intergenic (+71/‑368)PP_4432/PP_4433Xaa‑Pro aminopeptidase/amino acid MFS transporter
Reads supporting (aligned to +/- strand):  ref base A (38/41);  new base G (3/3);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.27e-01

ACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGAAGTAATAC  >  NC_002947/5028894‑5029187
                                                                                                                                                    |                                                                                                                                                 
aCCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAg                                                                                                                                                  <  2:1407929/150‑1 (MQ=255)
 ccGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGc                                                                                                                                                 >  1:664107/1‑150 (MQ=255)
      agCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCTTTAAGGCACCAt                                                                                                                                            <  1:524788/150‑1 (MQ=255)
      agCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGAATAATCGGCTGCTCTTTGTCAAGGACTGCCGCAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCTTTAAGGCACCAt                                                                                                                                            <  1:1869061/150‑1 (MQ=255)
         aTCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATATGATGGTGCCTTAAAGCACCATCAc                                                                                                                                         >  1:1333744/1‑150 (MQ=255)
          tCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCTTTAAGGCACCATCAca                                                                                                                                        <  1:1447241/150‑2 (MQ=255)
            cTGATCACTGATACTGGGGTCGAGGCGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTcc                                                                                                                                      >  1:919437/1‑150 (MQ=255)
             tGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCACGCTTTGCGTCATGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCa                                                                                                                                     <  2:530094/150‑1 (MQ=255)
              gATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAg                                                                                                                                    <  1:1291863/150‑1 (MQ=255)
              gATCACTGATACTGGGGTCGAGACGATCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAg                                                                                                                                    >  1:991015/1‑150 (MQ=255)
                    tGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCAc                                                                                                                                <  2:863347/148‑1 (MQ=255)
                    tGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCAc                                                                                                                                >  1:863347/1‑148 (MQ=255)
                       tACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACCGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCtt                                                                                                                           >  2:1017322/1‑150 (MQ=255)
                        aCTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCttt                                                                                                                          <  2:131501/150‑1 (MQ=255)
                           gtggTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGtt                                                                                                                       <  1:1624074/148‑1 (MQ=255)
                           ggggTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATATGGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGtt                                                                                                                       <  2:699552/150‑1 (MQ=255)
                           ggggTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGGCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACCCCAGGCACTCCTTTGtt                                                                                                                       >  1:1029972/1‑150 (MQ=255)
                                cGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACa                                                                                                                  <  1:907895/150‑1 (MQ=255)
                                 gagaCGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACaa                                                                                                                 >  2:1596398/1‑150 (MQ=255)
                                   gaCGTTCTGCTCGACTCCTCGCAGGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACaaaa                                                                                                               >  2:946837/1‑150 (MQ=255)
                                       ttCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAg                                                                                                           <  2:1379227/150‑1 (MQ=255)
                                            cTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGc                                                                                                      >  1:1846494/1‑150 (MQ=255)
                                              cGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAATGTTTGAACGTCGCCGGCGTCATGGTGCCACGCCTAGCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAg                                                                                                    >  1:514371/1‑150 (MQ=255)
                                               gACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGc                                                                                                   >  2:9765/1‑150 (MQ=255)
                                               gACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGc                                                                                                   >  2:1209666/1‑150 (MQ=255)
                                                aCTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGcc                                                                                                  >  2:802406/1‑150 (MQ=255)
                                                     tcGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGaa                                                                                             <  1:1596398/150‑1 (MQ=255)
                                                     tcGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGaa                                                                                             <  1:9765/150‑1 (MQ=255)
                                                       gCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAaga                                                                                           >  2:769575/1‑150 (MQ=255)
                                                       gCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAaga                                                                                           <  2:514371/150‑1 (MQ=255)
                                                         aaGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAgg                                                                                         <  2:664107/150‑1 (MQ=255)
                                                         aaGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAgg                                                                                         >  1:893757/1‑150 (MQ=255)
                                                          aGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGt                                                                                        <  2:1761342/150‑1 (MQ=255)
                                                             ttttCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGa                                                                                     >  1:1786555/1‑150 (MQ=255)
                                                               ttCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGaaa                                                                                   <  1:310709/150‑1 (MQ=255)
                                                               ttCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGaaa                                                                                   <  1:1209666/150‑1 (MQ=255)
                                                                    cAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGa                                                                              >  2:314600/1‑150 (MQ=255)
                                                                      aGGACTGACGAGAGGTACGCAAGTAGCGGGAGTTTTAGTGCCACGGCTATCTTCAGAGACAATGTGATGGTGCATCAAAGAACCATCACGCTAGCCCCTCCTTTATTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGAtt                                                                            <  2:216055/150‑1 (MQ=255)
                                                                            gACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGTGACAAGGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTccc                                                                      <  2:1296609/150‑1 (MQ=255)
                                                                              cGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCa                                                                    <  1:769575/150‑1 (MQ=255)
                                                                               gAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAg                                                                   >  2:81652/1‑150 (MQ=255)
                                                                                aaaaGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCTCCAgg                                                                  >  2:1624210/1‑150 (MQ=255)
                                                                                aaaaGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAgg                                                                  <  1:1372482/150‑1 (MQ=255)
                                                                                 aaaGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAggg                                                                 <  1:953640/150‑1 (MQ=255)
                                                                                 aaaGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAggg                                                                 >  2:106484/1‑150 (MQ=255)
                                                                                    gTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGc                                                              <  2:1640446/150‑1 (MQ=255)
                                                                                    gTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGc                                                              >  1:1640446/1‑150 (MQ=255)
                                                                                     tttGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTg                                                                                                                         >  2:791863/1‑90 (MQ=255)
                                                                                     tttGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTg                                                                                                                         <  1:791863/90‑1 (MQ=255)
                                                                                     tttGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCg                                                             >  2:137661/1‑150 (MQ=255)
                                                                                      ttGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGCAGTTGAAAGTTGATTGTTCCCCAGGGAGCGa                                                            <  1:816138/150‑1 (MQ=255)
                                                                                      ttGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGa                                                            <  2:1149044/150‑1 (MQ=255)
                                                                                       tggaaGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAAGGTGATGGTGCCTTAAAGCAACATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAg                                                           <  2:38831/147‑1 (MQ=255)
                                                                                            gTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGc                                                                                                          >  2:1259035/1‑98 (MQ=255)
                                                                                            gTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGc                                                                                                          <  1:1259035/98‑1 (MQ=255)
                                                                                               gCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGt                                                          <  1:1628833/143‑1 (MQ=255)
                                                                                               gCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGt                                                          >  2:1628833/1‑143 (MQ=255)
                                                                                                cTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCtgt                                                  <  2:1846494/150‑1 (MQ=255)
                                                                                                 tGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCtgtg                                                 >  2:1558671/1‑150 (MQ=255)
                                                                                                    aGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCt                                              >  1:753609/1‑150 (MQ=255)
                                                                                                     gTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCt                                                                                                                            <  1:1140320/71‑1 (MQ=255)
                                                                                                     gTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCt                                                                                                                            >  2:1140320/1‑71 (MQ=255)
                                                                                                     gTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTa                                             >  2:1502960/1‑150 (MQ=255)
                                                                                                      ttATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTaa                                            >  1:1288276/1‑150 (MQ=255)
                                                                                                       tataGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAAc                                           <  1:240877/150‑1 (MQ=255)
                                                                                                       tataGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAAc                                           <  1:1558671/150‑1 (MQ=255)
                                                                                                        ataGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACg                                          <  2:1288276/150‑1 (MQ=255)
                                                                                                              ccACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGaa                                    >  1:153274/1‑150 (MQ=255)
                                                                                                              ccACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGa                                     <  2:170746/149‑1 (MQ=255)
                                                                                                              ccACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGa                                     >  1:170746/1‑149 (MQ=255)
                                                                                                                aCGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAAc                                  >  2:1863102/1‑150 (MQ=255)
                                                                                                                    cTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTAt                              <  1:1502960/150‑1 (MQ=255)
                                                                                                                    cTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTAt                              <  1:1017322/150‑1 (MQ=255)
                                                                                                                     tATCTTCAGCGACAATGTGATGGTGCCTTAAAGCCCCCTCCCGCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATc                             <  1:959165/150‑1 (MQ=255)
                                                                                                                     tATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATc                             <  1:802406/150‑1 (MQ=255)
                                                                                                                               gACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGacgacg                   >  1:272938/1‑150 (MQ=255)
                                                                                                                                    tgtgATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGaa              <  2:1247104/150‑1 (MQ=255)
                                                                                                                                     gtgATGGTGCTTTAAGGCACCATCACTCCAGCCACTCCTTTGTTTTAGAAAAAAAGCACGCAGCCTTGAGGAGGTTGAAAGTTGATTGTTCCCCAGGGAGTGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGATGCTGAAt             >  2:977619/1‑150 (MQ=255)
                                                                                                                                     gtgATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAAt             <  1:106484/150‑1 (MQ=255)
                                                                                                                                       gATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAg           <  2:893757/150‑1 (MQ=255)
                                                                                                                                        aTGGTGCTTTAAGGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGa          >  2:186655/1‑150 (MQ=255)
                                                                                                                                           gTGCTTTAAGGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGAAGt       >  2:1131071/1‑150 (MQ=255)
                                                                                                                                            tGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGAAGTa      >  1:1264164/1‑150 (MQ=255)
                                                                                                                                                ttAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGAAGTAATAc  >  2:236265/1‑150 (MQ=255)
                                                                                                                                                ttAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGAAGTAATAc  >  2:1501139/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGAAGTAATAC  >  NC_002947/5028894‑5029187

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