Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,097,943 T→A 100% N19Y (AAT→TAT)  csgE ← curlin secretion specificity factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,097,9430TA95.3% 745.0 / 18.5 235N19Y (AAT→TAT) csgEcurlin secretion specificity factor
Reads supporting (aligned to +/- strand):  ref base T (5/6);  new base A (122/102);  total (127/108)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 236 total aligned reads displayed.

GCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGT  >  CP009273/1097798‑1098088
                                                                                                                                                 |                                                                                                                                                 
gCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATAcc                                                                                                                                                 <  1:1191891/148‑1 (MQ=255)
     gggCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTgcg                                                                                                                                            >  2:2881839/1‑148 (MQ=255)
      ggCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTgcgg                                                                                                                                           >  2:2286798/1‑148 (MQ=255)
       gCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTgcggc                                                                                                                                          >  2:3044833/1‑148 (MQ=255)
        ccTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTgcggcg                                                                                                                                         >  1:1019134/1‑148 (MQ=255)
           ttCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTCTGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGacc                                                                                                                                      >  2:358419/1‑146 (MQ=255)
           ttCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAAc                                                                                                                                      >  1:251233/1‑148 (MQ=255)
           ttCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAAc                                                                                                                                      >  2:1896874/1‑148 (MQ=255)
            tCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACa                                                                                                                                     >  1:2380004/1‑148 (MQ=255)
             cATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAg                                                                                                                                    >  1:2269043/1‑148 (MQ=255)
              aTTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGa                                                                                                                                   >  1:538114/1‑148 (MQ=255)
               ttAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGaa                                                                                                                                  >  1:1103109/1‑148 (MQ=255)
               ttAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGaa                                                                                                                                  <  1:2872033/148‑1 (MQ=255)
                  aTCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAAtt                                                                                                                               <  1:1715160/148‑1 (MQ=255)
                   tCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTc                                                                                                                              >  2:2578870/1‑148 (MQ=255)
                     gTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTg                                                                                                                            <  1:2934916/148‑1 (MQ=255)
                             aCCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCcaca                                                                                                                    <  1:1977407/148‑1 (MQ=255)
                               ccGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAAt                                                                                                                  >  1:2887071/1‑148 (MQ=255)
                                   ataGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAg                                                                                                              <  2:2293209/148‑1 (MQ=255)
                                       tCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTa                                                                                                          <  2:751257/148‑1 (MQ=255)
                                        cACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTaa                                                                                                         <  2:1382880/148‑1 (MQ=255)
                                        cACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTaa                                                                                                         <  1:125161/148‑1 (MQ=255)
                                        cACTTTCCCATTTAACACTAAAGGCTCGGTAAAAATCAACTCCAAAATAGCAAACAGAAAGGTCACTTAGAAATCCCGGGACTTCTACCTCAACGGCGTTAAGATACCCTCCGGCGAACACAAATTCAGCCGCCACAATCCAGccaaa                                                                                                         >  2:1687668/1‑144 (MQ=255)
                                          cTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAAt                                                                                                       >  1:602832/1‑148 (MQ=255)
                                           tttCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATa                                                                                                      >  2:1700021/1‑148 (MQ=255)
                                            ttCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATaa                                                                                                     >  1:1894064/1‑148 (MQ=255)
                                            ttCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATaa                                                                                                     >  2:1037715/1‑148 (MQ=255)
                                              cccATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACg                                                                                                   <  2:3016170/148‑1 (MQ=255)
                                               ccATTTATCACTAAAGGCTCGGTAAAAATCATGGACAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATACTGCCGCCACAATCCAGCGTAAATAACGt                                                                                                  >  1:3263952/1‑148 (MQ=255)
                                                  tttATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTc                                                                                               >  1:2463511/1‑148 (MQ=255)
                                                  tttATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTc                                                                                               <  1:66177/148‑1 (MQ=255)
                                                  tttATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTc                                                                                               >  2:749681/1‑148 (MQ=255)
                                                   ttATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCa                                                                                              >  2:1540981/1‑148 (MQ=255)
                                                    tATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCAt                                                                                             <  2:2140165/148‑1 (MQ=255)
                                                      tCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATgg                                                                                           <  2:1215237/148‑1 (MQ=255)
                                                      tCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATgg                                                                                           <  1:461143/148‑1 (MQ=255)
                                                       cACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGc                                                                                          <  2:96940/148‑1 (MQ=255)
                                                        aCTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCt                                                                                         <  2:1666850/148‑1 (MQ=255)
                                                         cTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCtt                                                                                        >  1:2139508/1‑148 (MQ=255)
                                                              ggCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCg                                                                                   >  2:965985/1‑148 (MQ=255)
                                                              ggCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCg                                                                                   >  2:2986847/1‑148 (MQ=255)
                                                              ggCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCg                                                                                   >  2:3265488/1‑148 (MQ=255)
                                                                       aaaaTCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGtt                                                                          >  1:1359115/1‑148 (MQ=255)
                                                                       aaaaTCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGtt                                                                          >  1:2050582/1‑148 (MQ=255)
                                                                       aaaaTCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGtt                                                                          >  1:3106404/1‑148 (MQ=255)
                                                                        aaaTCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTa                                                                         <  2:2570745/148‑1 (MQ=255)
                                                                         aaTCATGGCCAATAGAAGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGCTACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTAt                                                                        <  2:1965014/148‑1 (MQ=255)
                                                                          aTCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATc                                                                       >  1:2509785/1‑148 (MQ=255)
                                                                          aTCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATc                                                                       >  2:3442884/1‑148 (MQ=255)
                                                                            cATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGt                                                                     <  1:408274/148‑1 (MQ=255)
                                                                            cATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGt                                                                     >  1:366924/1‑148 (MQ=255)
                                                                                gCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGc                                                                 >  2:656416/1‑148 (MQ=255)
                                                                                   aaTAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCa                                                              <  1:2014351/148‑1 (MQ=255)
                                                                                    aTAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAg                                                             >  2:2871294/1‑148 (MQ=255)
                                                                                    aTAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAg                                                             <  2:1648231/148‑1 (MQ=255)
                                                                                        aTGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGaaa                                                         <  2:1548050/148‑1 (MQ=255)
                                                                                         tGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAAc                                                        >  2:3383226/1‑148 (MQ=255)
                                                                                          gAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAAcc                                                       >  1:2099544/1‑148 (MQ=255)
                                                                                            aaCAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGc                                                     <  1:1932818/148‑1 (MQ=255)
                                                                                            aaCAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGc                                                     >  1:2507852/1‑148 (MQ=255)
                                                                                             aCAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCt                                                    <  2:773496/148‑1 (MQ=255)
                                                                                             aCAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCt                                                    >  2:2221214/1‑148 (MQ=255)
                                                                                                 tATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTtgtg                                                >  2:396000/1‑148 (MQ=255)
                                                                                                  aTGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTgtgt                                               <  2:1833106/148‑1 (MQ=255)
                                                                                                  aTGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTgtgt                                               >  1:835484/1‑148 (MQ=255)
                                                                                                       cAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGt                                          >  1:2979509/1‑148 (MQ=255)
                                                                                                        aGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGtt                                         <  1:406419/148‑1 (MQ=255)
                                                                                                        aGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGtt                                         <  1:406445/148‑1 (MQ=255)
                                                                                                          ttAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGtttt                                       <  1:3035073/148‑1 (MQ=255)
                                                                                                          ttAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGtttt                                       <  2:2380004/148‑1 (MQ=255)
                                                                                                            aGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGtttttt                                     <  1:750631/148‑1 (MQ=255)
                                                                                                            aGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGtttttt                                     <  1:965985/148‑1 (MQ=255)
                                                                                                            aGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGtttttt                                     <  2:2726289/148‑1 (MQ=255)
                                                                                                             gCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTa                                    <  2:835104/148‑1 (MQ=255)
                                                                                                              cAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTAc                                   >  1:1676167/1‑148 (MQ=255)
                                                                                                                aTCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACgg                                 <  1:1428749/148‑1 (MQ=255)
                                                                                                                 tCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTGGTGTCCGGTTTTTTACGGc                                >  1:1488948/1‑148 (MQ=255)
                                                                                                                 tCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGc                                >  1:353587/1‑148 (MQ=255)
                                                                                                                 tCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGc                                >  1:317589/1‑148 (MQ=255)
                                                                                                                  cccGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCt                               <  2:2233715/148‑1 (MQ=255)
                                                                                                                   ccGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTa                              <  1:2820366/148‑1 (MQ=255)
                                                                                                                   ccGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTa                              <  1:677392/148‑1 (MQ=255)
                                                                                                                   ccGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTa                              <  2:2463511/148‑1 (MQ=255)
                                                                                                                     gggACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTAtc                            >  1:2804447/1‑148 (MQ=255)
                                                                                                                       gACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTAtctt                          <  1:3148644/148‑1 (MQ=255)
                                                                                                                        aCTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCAGGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTAtcttc                         <  2:2973889/148‑1 (MQ=255)
                                                                                                                           tCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTg                      <  1:154568/148‑1 (MQ=255)
                                                                                                                           tCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTg                      <  2:2084527/148‑1 (MQ=255)
                                                                                                                           tCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTg                      <  2:1296261/148‑1 (MQ=255)
                                                                                                                              aCCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGaaa                   <  2:1103109/148‑1 (MQ=255)
                                                                                                                               ccTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGaaaa                  >  2:899728/1‑148 (MQ=255)
                                                                                                                               ccTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGaaaa                  >  1:1798450/1‑148 (MQ=255)
                                                                                                                                  cAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGAt               <  2:418285/148‑1 (MQ=255)
                                                                                                                                   aaCGGAGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGAtt              <  2:2826720/148‑1 (MQ=255)
                                                                                                                                    aCGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATta             >  2:1620644/1‑148 (MQ=255)
                                                                                                                                       gCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATaa          <  2:2511734/148‑1 (MQ=255)
                                                                                                                                       gCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATaa          >  2:3164343/1‑148 (MQ=255)
                                                                                                                                          tGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGa       >  1:1377917/1‑148 (MQ=255)
                                                                                                                                          tGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGa       >  1:2662639/1‑148 (MQ=255)
                                                                                                                                          tGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGa       >  1:2659701/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGT  >  CP009273/1097798‑1098088

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: