Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,113,450 A→T 100% I161I (ATT→ATA)  yceI ← secreted protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,113,4500AT97.3% 713.9 / 7.0 220I161I (ATT→ATA) yceIsecreted protein
Reads supporting (aligned to +/- strand):  ref base A (2/4);  new base T (135/79);  total (137/83)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 220 total aligned reads displayed.

CAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTAATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCACCGTCTTT  >  CP009273/1113305‑1113593
                                                                                                                                                 |                                                                                                                                               
cAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTAt                                                                                                                                               <  2:1016386/148‑1 (MQ=255)
 aaCAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTAtt                                                                                                                                              >  2:785512/1‑148 (MQ=255)
  aCAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTAttt                                                                                                                                             >  1:2301371/1‑148 (MQ=255)
  aCAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTAttt                                                                                                                                             >  2:527402/1‑148 (MQ=255)
   cAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAATTCTGTCTTGATCTTGAAGTCCTTGAGCTTTAtttt                                                                                                                                            >  2:2888001/1‑148 (MQ=255)
    aGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTg                                                                                                                                           >  2:3661788/1‑148 (MQ=255)
     gAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGACGTCCGTGAGCTTTATTTTGc                                                                                                                                          >  2:3328354/1‑148 (MQ=255)
      aataaataTCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGcc                                                                                                                                         <  1:213918/148‑1 (MQ=255)
      aataaataTCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGcc                                                                                                                                         >  1:2762964/1‑148 (MQ=255)
       ataaataTCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCt                                                                                                                                        >  2:342867/1‑148 (MQ=255)
        taaataTCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCtt                                                                                                                                       <  2:1572658/148‑1 (MQ=255)
        taaataTCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCtt                                                                                                                                       >  1:114724/1‑148 (MQ=255)
           atatCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCgg                                                                                                                                    <  2:1896902/148‑1 (MQ=255)
                aTTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCg                                                                                                                               >  2:1743945/1‑148 (MQ=255)
                 ttGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGa                                                                                                                              <  1:2702924/148‑1 (MQ=255)
                 ttGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGa                                                                                                                              >  2:3728692/1‑148 (MQ=255)
                 ttGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTAATTTTGCCTTCGGCCTCGa                                                                                                                              <  1:3117749/148‑1 (MQ=255)
                    ccGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGc                                                                                                                           <  1:2265789/148‑1 (MQ=255)
                         tGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCa                                                                                                                      >  1:2876308/1‑148 (MQ=255)
                             tGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTAATTTTGCCTTCGGCCTCGAAGCCTGCACGtt                                                                                                                  >  1:1460669/1‑148 (MQ=255)
                               ccATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTa                                                                                                                >  2:2885610/1‑148 (MQ=255)
                                cATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTac                                                                                                               >  2:746039/1‑148 (MQ=255)
                                   ccGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTaccac                                                                                                            >  2:3277025/1‑148 (MQ=255)
                                    cGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTaccacc                                                                                                           <  1:2018088/148‑1 (MQ=255)
                                       cAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCa                                                                                                        >  1:622611/1‑148 (MQ=255)
                                       cAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCa                                                                                                        >  2:3586047/1‑148 (MQ=255)
                                       cAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCa                                                                                                        >  1:1762231/1‑148 (MQ=255)
                                       cAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCa                                                                                                        >  1:3254729/1‑148 (MQ=255)
                                       cAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCa                                                                                                        >  1:1350459/1‑148 (MQ=255)
                                          cATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATgg                                                                                                     >  1:2737795/1‑148 (MQ=255)
                                           aTTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATggg                                                                                                    >  2:1124087/1‑148 (MQ=255)
                                                aCGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGgtcgtc                                                                                               >  2:1372052/1‑148 (MQ=255)
                                                aCGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGgtcgtc                                                                                               >  2:1372344/1‑148 (MQ=255)
                                                   cTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGcc                                                                                            >  2:2000309/1‑148 (MQ=255)
                                                   cTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGcc                                                                                            >  1:2666027/1‑148 (MQ=255)
                                                   cTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGcc                                                                                            >  2:3020852/1‑148 (MQ=255)
                                                   cTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCAATTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGcc                                                                                            <  1:112876/148‑1 (MQ=255)
                                                       cTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGa                                                                                        <  2:443966/148‑1 (MQ=255)
                                                                 aCGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCaattaatt                                                                              >  2:2124155/1‑148 (MQ=255)
                                                                      ttCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCt                                                                         <  1:3326711/148‑1 (MQ=255)
                                                                      ttCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCt                                                                         >  1:1519155/1‑148 (MQ=255)
                                                                       tCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCtt                                                                        >  1:2417801/1‑148 (MQ=255)
                                                                        ccACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTc                                                                       <  1:608693/148‑1 (MQ=255)
                                                                             gAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCg                                                                  >  2:2837510/1‑148 (MQ=255)
                                                                             gAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCg                                                                  >  2:2560106/1‑148 (MQ=255)
                                                                             gAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTAATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCg                                                                  <  1:1062120/148‑1 (MQ=255)
                                                                              aaaTAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGt                                                                 >  2:1505277/1‑148 (MQ=255)
                                                                               aataatCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTg                                                                <  1:89978/148‑1 (MQ=255)
                                                                                ataatCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGTCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGa                                                               <  2:2161913/148‑1 (MQ=255)
                                                                                   atCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACgg                                                            >  2:2417096/1‑148 (MQ=255)
                                                                                    tCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACggg                                                           <  1:713403/148‑1 (MQ=255)
                                                                                    tCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACggg                                                           >  1:3722476/1‑148 (MQ=255)
                                                                                    tCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACggg                                                           <  1:3593782/148‑1 (MQ=255)
                                                                                     cAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGt                                                          <  1:1268494/148‑1 (MQ=255)
                                                                                      aGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGtt                                                         >  2:3495017/1‑148 (MQ=255)
                                                                                       gATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGttt                                                        >  2:3047427/1‑148 (MQ=255)
                                                                                        aTCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTTAATTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTAGCCCTGACCAATTAATATCGCTTCCAGAGTGACGGGTTTg                                                       >  1:148688/1‑148 (MQ=255)
                                                                                        aTCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTg                                                       >  1:2465575/1‑148 (MQ=255)
                                                                                        aTCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTg                                                       >  1:2472175/1‑148 (MQ=255)
                                                                                        aTCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTg                                                       >  1:2615238/1‑148 (MQ=255)
                                                                                         tCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTgg                                                      >  2:1798807/1‑148 (MQ=255)
                                                                                         tCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTgg                                                      <  1:2825213/148‑1 (MQ=255)
                                                                                          ccACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGt                                                     <  1:3052387/148‑1 (MQ=255)
                                                                                           cACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGtt                                                    >  2:928255/1‑148 (MQ=255)
                                                                                            aCTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCACTTAATTTCGCTTCCAGCGTGACGGGTTTGGTTa                                                   <  2:2310849/148‑1 (MQ=255)
                                                                                               tCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGc                                                >  2:2449986/1‑148 (MQ=255)
                                                                                               tCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGc                                                >  1:3213253/1‑148 (MQ=255)
                                                                                                    agaAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTc                                           <  1:818443/148‑1 (MQ=255)
                                                                                                    agaAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTc                                           >  2:2537912/1‑148 (MQ=255)
                                                                                                        gCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCg                                       >  1:1978774/1‑148 (MQ=255)
                                                                                                        gCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCg                                       >  1:216862/1‑148 (MQ=255)
                                                                                                        gCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCg                                       >  1:217259/1‑148 (MQ=255)
                                                                                                        gCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCAGTCAGCg                                       >  1:1978705/1‑148 (MQ=255)
                                                                                                         cTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGt                                      >  2:2668265/1‑148 (MQ=255)
                                                                                                            gACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAg                                   >  2:571400/1‑148 (MQ=255)
                                                                                                             aCCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGa                                  <  1:1273699/148‑1 (MQ=255)
                                                                                                             aCCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGa                                  <  1:1703595/148‑1 (MQ=255)
                                                                                                              ccTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGAt                                 >  2:1857646/1‑148 (MQ=255)
                                                                                                                tAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCg                               >  2:2798899/1‑148 (MQ=255)
                                                                                                                  aaTCTGTCTTGATATTGAAGTCCTTGAGCTTAATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGcc                             >  1:518002/1‑148 (MQ=255)
                                                                                                                   aTCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCg                            >  1:3505767/1‑148 (MQ=255)
                                                                                                                      tGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTa                         <  1:917650/148‑1 (MQ=255)
                                                                                                                        tCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAat                       >  1:708701/1‑148 (MQ=255)
                                                                                                                        tCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAat                       >  1:2888981/1‑148 (MQ=255)
                                                                                                                        tCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAat                       >  1:2347949/1‑148 (MQ=255)
                                                                                                                             atatTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCa                  >  2:488908/1‑148 (MQ=255)
                                                                                                                              tatTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAg                 <  2:611939/148‑1 (MQ=255)
                                                                                                                              tatTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAg                 <  2:2358028/148‑1 (MQ=255)
                                                                                                                               atTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGTTTACGCCATTCAGAGTCAGATCGCCGTTAATATCCAGt                >  1:181090/1‑148 (MQ=255)
                                                                                                                               atTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGt                >  1:505653/1‑148 (MQ=255)
                                                                                                                               atTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGt                <  1:746045/148‑1 (MQ=255)
                                                                                                                                ttGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGtt               >  2:3162503/1‑148 (MQ=255)
                                                                                                                                ttGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGtt               >  1:2190901/1‑148 (MQ=255)
                                                                                                                                 tGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTc              >  1:837577/1‑148 (MQ=255)
                                                                                                                                     gTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCa          <  2:505056/148‑1 (MQ=255)
                                                                                                                                      tCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCAc         >  1:3125373/1‑148 (MQ=255)
                                                                                                                                          tGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCACCGtc     <  1:2850492/148‑1 (MQ=255)
                                                                                                                                             gCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCACCGtcttt  >  2:3300654/1‑148 (MQ=255)
                                                                                                                                             gCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCACCGtcttt  <  1:670422/148‑1 (MQ=255)
                                                                                                                                             gCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCACCGtcttt  >  2:1400656/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTAATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCACCGTCTTT  >  CP009273/1113305‑1113593

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