Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,181,504 G→T 100% D69Y (GAT→TAT)  pepT → peptidase T

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,181,5040GT99.1% 739.3 / ‑2.1 225D69Y (GAT→TAT) pepTpeptidase T
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base T (110/113);  total (112/113)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 225 total aligned reads displayed.

GGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGTTATGT  >  CP009273/1181362‑1181648
                                                                                                                                              |                                                                                                                                                
 gTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATcc                                                                                                                                            <  1:1077628/148‑1 (MQ=255)
  tGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATccc                                                                                                                                           <  2:669462/148‑1 (MQ=255)
         ggTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGAtt                                                                                                                                    <  2:2994385/148‑1 (MQ=255)
           ttCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTgg                                                                                                                                  <  2:750289/148‑1 (MQ=255)
           ttCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTgg                                                                                                                                  <  2:3444984/148‑1 (MQ=255)
            tCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGc                                                                                                                                 >  1:427946/1‑148 (MQ=255)
            tCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGc                                                                                                                                 <  1:75325/148‑1 (MQ=255)
             cccAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCt                                                                                                                                >  1:1716019/1‑148 (MQ=255)
              ccAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCtt                                                                                                                               >  1:2326831/1‑148 (MQ=255)
              ccAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCtt                                                                                                                               <  2:585949/148‑1 (MQ=255)
               cAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCttt                                                                                                                              <  2:662276/148‑1 (MQ=255)
                 gCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTAt                                                                                                                            >  2:1649146/1‑148 (MQ=255)
                  cACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTAtt                                                                                                                           <  1:205351/148‑1 (MQ=255)
                   aCGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTAttt                                                                                                                          <  2:237594/148‑1 (MQ=255)
                    cGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTtc                                                                                                                         >  2:2599825/1‑148 (MQ=255)
                    cGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTtc                                                                                                                         >  2:2332401/1‑148 (MQ=255)
                    cGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTtc                                                                                                                         <  2:2210972/148‑1 (MQ=255)
                    cGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTtc                                                                                                                         >  1:353479/1‑148 (MQ=255)
                    cGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTtc                                                                                                                         >  2:845124/1‑148 (MQ=255)
                         ggCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCAtg                                                                                                                    <  1:2222208/148‑1 (MQ=255)
                             aaTGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGa                                                                                                                <  1:2979038/148‑1 (MQ=255)
                               tGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATa                                                                                                              <  1:1710377/148‑1 (MQ=255)
                                    gTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCa                                                                                                         <  1:396394/148‑1 (MQ=255)
                                      tATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCAcc                                                                                                       >  1:549017/1‑148 (MQ=255)
                                      tATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCAcc                                                                                                       >  1:1748229/1‑148 (MQ=255)
                                         tGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGa                                                                                                    <  2:3188606/148‑1 (MQ=255)
                                         tGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGa                                                                                                    <  2:3188883/148‑1 (MQ=255)
                                           cATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGAtt                                                                                                  <  2:1518327/148‑1 (MQ=255)
                                            aTCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTg                                                                                                 >  2:332003/1‑148 (MQ=255)
                                              ctgctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCa                                                                                               >  2:1811581/1‑148 (MQ=255)
                                              ctgctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCa                                                                                               <  1:2591902/148‑1 (MQ=255)
                                              ctgctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCa                                                                                               <  1:894613/148‑1 (MQ=255)
                                               tgctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAg                                                                                              <  1:3301515/148‑1 (MQ=255)
                                                   gAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGc                                                                                          <  2:1664897/148‑1 (MQ=255)
                                                    aaaGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCa                                                                                         >  1:3267095/1‑148 (MQ=255)
                                                     aaGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCaa                                                                                        >  1:1996722/1‑148 (MQ=255)
                                                       gagCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCaaaa                                                                                      >  1:3180584/1‑148 (MQ=255)
                                                       gagCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCaaaa                                                                                      <  2:2570367/148‑1 (MQ=255)
                                                           agcTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAAtgt                                                                                  <  1:3093887/148‑1 (MQ=255)
                                                             cTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGa                                                                                >  2:430851/1‑148 (MQ=255)
                                                               cGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAAt                                                                              <  2:1518050/148‑1 (MQ=255)
                                                                 aaGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATcc                                                                            <  1:1125072/148‑1 (MQ=255)
                                                                  agagATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTTCAGCGGCAAAAATGTGAATCCg                                                                           >  2:1247063/1‑148 (MQ=255)
                                                                   gagaTGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGc                                                                          >  2:787862/1‑148 (MQ=255)
                                                                   gagaTGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGc                                                                          >  2:423660/1‑148 (MQ=255)
                                                                     gaTGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCaa                                                                        <  2:3429441/148‑1 (MQ=255)
                                                                       tGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAAt                                                                      <  2:2799331/148‑1 (MQ=255)
                                                                          ggCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAAttgt                                                                   <  1:2989750/148‑1 (MQ=255)
                                                                              tATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGaa                                                               >  1:2630174/1‑148 (MQ=255)
                                                                               aTCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCACATTGTTGaaa                                                              >  2:3115910/1‑148 (MQ=255)
                                                                               aTCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGaaa                                                              >  1:2451334/1‑148 (MQ=255)
                                                                                 cAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAAc                                                            >  1:1453380/1‑148 (MQ=255)
                                                                                 cAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAAc                                                            >  1:1453595/1‑148 (MQ=255)
                                                                                   aTGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTa                                                          >  2:3269992/1‑148 (MQ=255)
                                                                                   aTGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTa                                                          <  1:2398012/148‑1 (MQ=255)
                                                                                    tgtgACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTAt                                                         >  2:169020/1‑148 (MQ=255)
                                                                                    tgtgACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTAt                                                         >  2:423347/1‑148 (MQ=255)
                                                                                         ccTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCgg                                                    <  1:3485677/148‑1 (MQ=255)
                                                                                          ctaaaGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGt                                                   >  1:3214807/4‑148 (MQ=255)
                                                                                          cTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTTTTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGt                                                   >  2:1495209/1‑148 (MQ=255)
                                                                                          cTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGt                                                   <  2:767225/148‑1 (MQ=255)
                                                                                               gTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGa                                              <  1:3358500/148‑1 (MQ=255)
                                                                                                 gagaAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGata                                            <  1:274256/148‑1 (MQ=255)
                                                                                                 gagaAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGata                                            <  2:2998415/148‑1 (MQ=255)
                                                                                                  agaAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGTAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGatat                                           >  1:780367/1‑148 (MQ=255)
                                                                                                  agaAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGatat                                           <  1:1264303/148‑1 (MQ=255)
                                                                                                    aaGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTg                                         >  2:729406/1‑148 (MQ=255)
                                                                                                      gggCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTgcg                                       <  1:1649146/148‑1 (MQ=255)
                                                                                                        gCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCt                                     >  2:828647/1‑148 (MQ=255)
                                                                                                        gCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCt                                     >  2:2587941/1‑148 (MQ=255)
                                                                                                          aCTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTgg                                   <  2:3180584/148‑1 (MQ=255)
                                                                                                           cTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTggg                                  >  1:1680694/1‑148 (MQ=255)
                                                                                                             ttGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTa                                >  2:247193/1‑148 (MQ=255)
                                                                                                               gATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATc                              <  1:2355626/148‑1 (MQ=255)
                                                                                                                 tGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCgg                            <  2:2024637/148‑1 (MQ=255)
                                                                                                                 tGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCgg                            >  1:2018149/1‑148 (MQ=255)
                                                                                                                    cGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGa                         >  2:2393231/1‑148 (MQ=255)
                                                                                                                    cGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGa                         >  2:1857663/1‑148 (MQ=255)
                                                                                                                       cGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGAGGa                      >  2:2047384/1‑148 (MQ=255)
                                                                                                                        gTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGaa                     <  2:1398285/148‑1 (MQ=255)
                                                                                                                         ttACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATAACTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGTAg                    >  1:192753/1‑148 (MQ=255)
                                                                                                                           aCCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGtt                  <  1:587454/148‑1 (MQ=255)
                                                                                                                           aCCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATACGCAAATTGTTGAAAACTATCGCGGTGGGGATATTGCGCTGGGTATCGGCGATGAAGtt                  >  1:1290651/1‑148 (MQ=255)
                                                                                                                           aCCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCACACCGGATTGCAGCGGCAAAAATGTGAAACCGCCAATTTTTGAAAAATATCGCGGGGGCGATATTTCGCTGGGTATCGGCGATGACGtt                  >  1:2522097/1‑148 (MQ=255)
                                                                                                                             cGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGtttt                <  1:1002133/148‑1 (MQ=255)
                                                                                                                             cGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGtttt                <  1:423628/148‑1 (MQ=255)
                                                                                                                              ggCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTa               <  1:3107612/148‑1 (MQ=255)
                                                                                                                                  aaCGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTAGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGAAGAAGTTTTATAAc           >  2:2774106/1‑148 (MQ=255)
                                                                                                                                   aCGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCAcc          >  1:2241363/1‑148 (MQ=255)
                                                                                                                                     gTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCgg        <  2:2842782/148‑1 (MQ=255)
                                                                                                                                       cccTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGtt      <  1:840825/148‑1 (MQ=255)
                                                                                                                                       cccTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGtt      <  1:2599825/148‑1 (MQ=255)
                                                                                                                                       cccTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGAGGAAGTTTTATCAACGGtg      >  1:2812601/1‑147 (MQ=255)
                                                                                                                                       cccTGGCTATATCCCGGCGATTGGCTTTAGTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGtt      <  2:2677454/148‑1 (MQ=255)
                                                                                                                                        ccTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGTTa     <  2:2986531/148‑1 (MQ=255)
                                                                                                                                        ccTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGTTa     <  1:368893/148‑1 (MQ=255)
                                                                                                                                          tGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGTTATg   <  2:2966856/148‑1 (MQ=255)
                                                                                                                                          tGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGTTATg   >  1:288045/1‑148 (MQ=255)
                                                                                                                                           ggCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGTTATGt  <  1:3146835/148‑1 (MQ=255)
                                                                                                                                           ggCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGAGATTGCGCTGGGTATCGGCGATGAAGTTTTATCAACGGTTATGt  >  1:2162327/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                
GGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGTTATGT  >  CP009273/1181362‑1181648

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