Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,430,873 G→T 100% S93* (TCA→TAA)  ompN ← outer membrane pore protein N, non‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,430,8730GT95.7% 619.8 / 17.6 208S93* (TCA→TAA) ompNouter membrane pore protein N, non‑specific
Reads supporting (aligned to +/- strand):  ref base G (4/5);  new base T (116/83);  total (120/88)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 208 total aligned reads displayed.

GTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTATC  >  CP009273/1430728‑1431015
                                                                                                                                                 |                                                                                                                                              
gTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTAc                                                                                                                                              <  1:1395018/148‑1 (MQ=255)
  cTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTg                                                                                                                                            <  1:150423/148‑1 (MQ=255)
  cTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTg                                                                                                                                            >  2:1520150/1‑148 (MQ=255)
    gCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGt                                                                                                                                          >  1:1644606/1‑148 (MQ=255)
     cATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGtt                                                                                                                                         <  2:2001655/148‑1 (MQ=255)
        tGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGttttt                                                                                                                                      <  1:1624087/148‑1 (MQ=255)
          gTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCAATAATTACGTCCTTAATCGAAAGAACCTTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTg                                                                                                                                    <  1:3050552/148‑1 (MQ=255)
           tataAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGa                                                                                                                                   >  2:289579/1‑148 (MQ=255)
           tataAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGa                                                                                                                                   >  2:2051464/1‑148 (MQ=255)
            ataAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGaa                                                                                                                                  >  2:1239443/1‑148 (MQ=255)
            ataAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGaa                                                                                                                                  <  2:1644606/148‑1 (MQ=255)
             taAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAg                                                                                                                                 <  2:1450764/148‑1 (MQ=255)
               agagTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGAt                                                                                                                               <  2:582232/148‑1 (MQ=255)
               agagTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGAt                                                                                                                               <  1:3142689/148‑1 (MQ=255)
               agagTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGAt                                                                                                                               <  1:1184372/148‑1 (MQ=255)
                 agTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTc                                                                                                                             >  1:748704/1‑148 (MQ=255)
                 agTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTc                                                                                                                             >  1:2219330/1‑148 (MQ=255)
                   tCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAg                                                                                                                           <  2:3353000/148‑1 (MQ=255)
                      ccGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAgtgt                                                                                                                        <  2:946632/148‑1 (MQ=255)
                       cGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGtgtt                                                                                                                       >  1:36650/1‑148 (MQ=255)
                       cGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGtgtt                                                                                                                       <  1:2210778/148‑1 (MQ=255)
                            aaTTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTg                                                                                                                  >  1:3184840/1‑148 (MQ=255)
                            aaTTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTg                                                                                                                  >  1:3139025/1‑148 (MQ=255)
                            aaTTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTg                                                                                                                  <  1:3005220/148‑1 (MQ=255)
                              ttCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGcc                                                                                                                >  1:3387790/1‑148 (MQ=255)
                              ttCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGcc                                                                                                                >  1:69379/1‑148 (MQ=255)
                                 aGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGa                                                                                                             >  1:2390442/1‑148 (MQ=255)
                                   gcagcaTATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAat                                                                                                           >  2:3377310/1‑148 (MQ=255)
                                        atatCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGt                                                                                                      >  1:2646646/1‑148 (MQ=255)
                                         tatCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTa                                                                                                     <  1:256030/148‑1 (MQ=255)
                                          atCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTAt                                                                                                    >  2:2942078/1‑148 (MQ=255)
                                           tCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTAtt                                                                                                   <  2:2861738/148‑1 (MQ=255)
                                             ggTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTcc                                                                                                 <  2:3507556/148‑1 (MQ=255)
                                               tCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCa                                                                                               >  2:114191/1‑148 (MQ=255)
                                               tCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCa                                                                                               >  2:1937606/1‑148 (MQ=255)
                                               tCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCa                                                                                               >  1:2699914/1‑148 (MQ=255)
                                                ccATCCTTCGATGTAGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGAAATCTGCAATTTTCAGCCCGGCAAATGACAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGTATATTGTATTGCCAt                                                                                              <  2:1762489/148‑1 (MQ=255)
                                                ccAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCAt                                                                                              <  2:1218049/148‑1 (MQ=255)
                                                  aGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTg                                                                                            <  2:2322306/148‑1 (MQ=255)
                                                  aGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTg                                                                                            >  1:283293/1‑148 (MQ=255)
                                                  aGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTg                                                                                            >  1:1083245/1‑148 (MQ=255)
                                                     cTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGAcc                                                                                         <  1:2100136/148‑1 (MQ=255)
                                                     cTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGAcc                                                                                         >  2:1003214/1‑148 (MQ=255)
                                                      ttCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCg                                                                                        >  1:1674317/1‑148 (MQ=255)
                                                            gTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCa                                                                                  >  1:3526219/1‑148 (MQ=255)
                                                              cGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGt                                                                                >  2:549227/1‑148 (MQ=255)
                                                                tACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGa                                                                              >  2:2002159/1‑148 (MQ=255)
                                                                tACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGa                                                                              <  1:1998644/148‑1 (MQ=255)
                                                                 aCATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAg                                                                             >  2:3209109/1‑148 (MQ=255)
                                                                     tACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGa                                                                         >  1:237889/1‑148 (MQ=255)
                                                                     tACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGa                                                                         <  1:2876138/148‑1 (MQ=255)
                                                                         ccATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGt                                                                     <  2:3393654/148‑1 (MQ=255)
                                                                          cATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGtt                                                                    >  2:83333/1‑148 (MQ=255)
                                                                          cATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGtt                                                                    >  2:942753/1‑148 (MQ=255)
                                                                            tAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTaa                                                                  <  1:1031661/148‑1 (MQ=255)
                                                                            tAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTaa                                                                  <  1:306410/148‑1 (MQ=255)
                                                                            tAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTaa                                                                  >  1:2957767/1‑148 (MQ=255)
                                                                             aaTTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAAt                                                                 >  2:86340/1‑148 (MQ=255)
                                                                              aTTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAAtt                                                                >  2:3205049/1‑148 (MQ=255)
                                                                                tACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTg                                                              >  1:2667662/1‑148 (MQ=255)
                                                                                    tCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGtt                                                          >  2:2671190/1‑148 (MQ=255)
                                                                                     ccGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGttt                                                         >  1:1085421/1‑148 (MQ=255)
                                                                                        tAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGc                                                      >  2:489362/1‑148 (MQ=255)
                                                                                        tAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGc                                                      >  2:800553/1‑148 (MQ=255)
                                                                                         aaTCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTACCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGcc                                                     >  1:2142083/1‑148 (MQ=255)
                                                                                           tCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTATCCAGTGAGTTTATCTTTAATTTGGGTTTCGCCtt                                                   >  1:1551734/1‑148 (MQ=255)
                                                                                           tCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCtt                                                   >  2:2652981/1‑148 (MQ=255)
                                                                                            cGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCttt                                                  >  2:934822/1‑148 (MQ=255)
                                                                                            cGAAAGAACCGAAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTATGAATATTCAGTGTTGTTTGACAGAATATTGTATTCCCATTTACCGTAGCCAGTGAGTTGAACGTTAATTTGGGTTTCGCCttt                                                  >  2:563401/1‑148 (MQ=255)
                                                                                                aGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTaaaa                                              <  1:2077256/148‑1 (MQ=255)
                                                                                                    ccGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCa                                          <  1:505631/148‑1 (MQ=255)
                                                                                                     cGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAg                                         <  1:1906143/148‑1 (MQ=255)
                                                                                                     cGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAg                                         <  1:1239443/148‑1 (MQ=255)
                                                                                                      gTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGa                                        >  2:246025/1‑148 (MQ=255)
                                                                                                      gTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGa                                        <  2:2814901/148‑1 (MQ=255)
                                                                                                        aaTCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGAcg                                      <  2:356583/148‑1 (MQ=255)
                                                                                                          tCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGAcgcg                                    >  1:2355116/1‑148 (MQ=255)
                                                                                                            tGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCa                                  <  1:3219162/148‑1 (MQ=255)
                                                                                                             gCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCAt                                 >  1:502775/1‑148 (MQ=255)
                                                                                                                aaTTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGc                              >  2:2558146/1‑148 (MQ=255)
                                                                                                                 aTTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGct                             >  2:3199865/1‑148 (MQ=255)
                                                                                                                 aTTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGct                             <  2:71180/148‑1 (MQ=255)
                                                                                                                           ggCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCa                   >  2:841249/1‑148 (MQ=255)
                                                                                                                            gCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAt                  <  2:3139025/148‑1 (MQ=255)
                                                                                                                            gCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAt                  >  2:182024/1‑148 (MQ=255)
                                                                                                                             cAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATc                 >  2:1685464/1‑148 (MQ=255)
                                                                                                                             cAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATc                 >  2:3446346/1‑148 (MQ=255)
                                                                                                                              aaaTGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCt                <  2:2219330/148‑1 (MQ=255)
                                                                                                                               aaTGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCtt               >  1:810028/1‑148 (MQ=255)
                                                                                                                               aaTGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCtt               <  2:1486480/148‑1 (MQ=255)
                                                                                                                                aTGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCttt              >  2:3510346/1‑148 (MQ=255)
                                                                                                                                  gCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCg            >  2:2545685/1‑148 (MQ=255)
                                                                                                                                   ccAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGc           >  1:573257/1‑148 (MQ=255)
                                                                                                                                   ccAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGc           >  2:1761973/1‑148 (MQ=255)
                                                                                                                                     aGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTg         <  1:2431393/148‑1 (MQ=255)
                                                                                                                                     aGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTg         >  2:746821/1‑148 (MQ=255)
                                                                                                                                       aCGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGaa       >  2:3072285/1‑148 (MQ=255)
                                                                                                                                         gggTCCATTACTGGTTTTTTGAAGTTTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAAtt     <  1:1776277/148‑1 (MQ=255)
                                                                                                                                          ggTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTa    <  2:210711/148‑1 (MQ=255)
                                                                                                                                            tCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTATc  >  1:2544491/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTATC  >  CP009273/1430728‑1431015

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