Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,029,037 C→T 7.2% intergenic (+66/‑373) 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,0370CT 7.2% 235.8 / 6.9 84intergenic (+66/‑373)PP_4432/PP_4433Xaa‑Pro aminopeptidase/amino acid MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (36/41);  new base T (3/3);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

GATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGAAGTA  >  NC_002947/5028890‑5029183
                                                                                                                                                   |                                                                                                                                                  
gaTTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGCTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCtt                                                                                                                                                  <  1:103846/150‑1 (MQ=255)
    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)
                                      agaCGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGcc                                                                                                                                                    <  1:1214365/110‑1 (MQ=255)
                                      agaCGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGcc                                                                                                                                                    >  2:1214365/1‑110 (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:514371/150‑1 (MQ=255)
                                                           gCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAaga                                                                                       >  2:769575/1‑150 (MQ=255)
                                                             aaGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAgg                                                                                     >  1:893757/1‑150 (MQ=255)
                                                             aaGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAgg                                                                                     <  2:664107/150‑1 (MQ=255)
                                                              aGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGt                                                                                    <  2:1761342/150‑1 (MQ=255)
                                                                 ttttCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGa                                                                                 >  1:1786555/1‑150 (MQ=255)
                                                                   ttCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGaaa                                                                               <  1:1209666/150‑1 (MQ=255)
                                                                   ttCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGaaa                                                                               <  1:310709/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                                                             >  2:106484/1‑150 (MQ=255)
                                                                                     aaaGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAggg                                                             <  1:953640/150‑1 (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                                                                                                      <  1:1259035/98‑1 (MQ=255)
                                                                                                gTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGc                                                                                                      >  2:1259035/1‑98 (MQ=255)
                                                                                                   gCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGt                                                      >  2:1628833/1‑143 (MQ=255)
                                                                                                   gCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGt                                                      <  1:1628833/143‑1 (MQ=255)
                                                                                                    cTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCtgt                                              <  2:1846494/150‑1 (MQ=255)
                                                                                                     tGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCtgtg                                             >  2:1558671/1‑150 (MQ=255)
                                                                                                        aGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCt                                          >  1:753609/1‑150 (MQ=255)
                                                                                                         gTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCt                                                                                                                        >  2:1140320/1‑71 (MQ=255)
                                                                                                         gTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCt                                                                                                                        <  1:1140320/71‑1 (MQ=255)
                                                                                                         gTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTa                                         >  2:1502960/1‑150 (MQ=255)
                                                                                                          ttATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTaa                                        >  1:1288276/1‑150 (MQ=255)
                                                                                                           tataGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAAc                                       <  1:1558671/150‑1 (MQ=255)
                                                                                                           tataGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAAc                                       <  1:240877/150‑1 (MQ=255)
                                                                                                            ataGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACg                                      <  2:1288276/150‑1 (MQ=255)
                                                                                                                  ccACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGaa                                >  1:153274/1‑150 (MQ=255)
                                                                                                                  ccACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGa                                 >  1:170746/1‑149 (MQ=255)
                                                                                                                  ccACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGa                                 <  2:170746/149‑1 (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)
                                                                                                                                                   |                                                                                                                                                  
GATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGAAGTA  >  NC_002947/5028890‑5029183

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: