Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,097,943 T→A 66.5% 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,9430TA66.5% 415.7 / inf 346N19Y (AAT→TAT) csgEcurlin secretion specificity factor
Reads supporting (aligned to +/- strand):  ref base T (61/55);  new base A (118/111);  total (180/166)
Fisher's exact test for biased strand distribution p-value = 9.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

Only 100 of 347 total aligned reads displayed.

CGTGCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGTTTT  >  CP009273/1097795‑1098091
                                                                                                                                                    |                                                                                                                                                    
  tGCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATAccc                                                                                                                                                   >  1:2108379/1‑150 (MQ=255)
   gCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCt                                                                                                                                                  >  2:1916581/1‑150 (MQ=255)
     aCTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTgc                                                                                                                                                >  1:1712043/1‑150 (MQ=255)
        gggCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTgcggc                                                                                                                                             >  1:3544700/1‑150 (MQ=255)
         ggCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTgcggcg                                                                                                                                            >  2:1628228/1‑150 (MQ=255)
         ggCCTTTCATTAATCGTTAAGATACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGAAAAAATCAAGGCCAAAAGAAGAAACAGTATGGTCAGTTATCAATCCCGGGAATACTAACAAAAATGCGTGAAGATTCCCTgcggcg                                                                                                                                            >  2:734305/1‑150 (MQ=255)
          gCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGa                                                                                                                                           >  1:1444891/1‑150 (MQ=255)
                 aTTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGaaa                                                                                                                                    <  2:171692/150‑1 (MQ=255)
                     aTCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCt                                                                                                                                >  2:3069770/1‑150 (MQ=255)
                     aTCGTTAAGTTACCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTg                                                                                                                               <  1:17475/150‑1 (MQ=255)
                      tCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTg                                                                                                                               >  2:3794254/1‑150 (MQ=255)
                        gTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTgcc                                                                                                                             <  1:2938170/150‑1 (MQ=255)
                         ttAAGTTACCCGTATAGTCACTGTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTgccg                                                                                                                            <  1:3760284/150‑1 (MQ=255)
                           aaGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTgccgcc                                                                                                                          >  2:2561801/1‑150 (MQ=255)
                             gTTACCCGTATATTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGACATTCTGCCGCcac                                                                                                                        >  2:1573706/1‑150 (MQ=255)
                                aCCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAAt                                                                                                                     >  2:3396690/1‑150 (MQ=255)
                                 cccGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATc                                                                                                                    <  1:660042/150‑1 (MQ=255)
                                  ccGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATcc                                                                                                                   <  1:2672360/150‑1 (MQ=255)
                                  ccGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATcc                                                                                                                   <  1:4073457/150‑1 (MQ=255)
                                   cGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCa                                                                                                                  <  2:1135102/150‑1 (MQ=255)
                                    gTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAg                                                                                                                 >  1:1915607/1‑150 (MQ=255)
                                    gTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAg                                                                                                                 >  2:4586111/1‑150 (MQ=255)
                                        aGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTa                                                                                                             >  1:4290164/1‑150 (MQ=255)
                                            aCTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATATATGAAACAGAATGGTCAGTTAGCAATCCCGGTACTTCTACATCAACGGCGTGAAGAAACCCTGCGGCGAACAGAAATTCTTCATCCACAAACCAGCGaaaaaa                                                                                                         >  2:182013/1‑148 (MQ=255)
                                              tttCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAAc                                                                                                       <  2:4053401/150‑1 (MQ=255)
                                              tttCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAAc                                                                                                       <  1:1563264/150‑1 (MQ=255)
                                               ttCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACg                                                                                                      >  2:2300780/1‑150 (MQ=255)
                                               ttCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACg                                                                                                      >  1:5003336/1‑150 (MQ=255)
                                                tCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGt                                                                                                     >  1:4730942/1‑150 (MQ=255)
                                                 cccATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGtt                                                                                                    <  2:3890107/150‑1 (MQ=255)
                                                  ccATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGttt                                                                                                   >  2:2420278/1‑150 (MQ=255)
                                                  ccATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGttt                                                                                                   <  2:1518322/150‑1 (MQ=255)
                                                   cATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTc                                                                                                  <  1:221159/150‑1 (MQ=255)
                                                   cATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTc                                                                                                  >  2:5158478/1‑150 (MQ=255)
                                                    aTTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCa                                                                                                 <  1:2247403/150‑1 (MQ=255)
                                                    aTTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCa                                                                                                 >  2:1081221/1‑150 (MQ=255)
                                                      ttATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATg                                                                                               >  1:63624/1‑150 (MQ=255)
                                                      ttATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATg                                                                                               <  2:4296711/150‑1 (MQ=255)
                                                       tATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATgg                                                                                              >  2:2479157/1‑150 (MQ=255)
                                                       tATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATgg                                                                                              <  2:3472080/150‑1 (MQ=255)
                                                         tCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCt                                                                                            >  2:1491588/1‑150 (MQ=255)
                                                          cACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCtt                                                                                           >  1:2502239/1‑150 (MQ=255)
                                                          cACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCtt                                                                                           >  1:3871185/1‑150 (MQ=255)
                                                          cACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCtt                                                                                           >  1:2891333/1‑150 (MQ=255)
                                                           aCTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCttt                                                                                          <  1:3794254/150‑1 (MQ=255)
                                                           aCTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCttt                                                                                          <  2:2881296/150‑1 (MQ=255)
                                                            cTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATACCGTTTCATGGCTTTa                                                                                         <  1:3220600/150‑1 (MQ=255)
                                                            cTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTa                                                                                         >  1:2962522/1‑150 (MQ=255)
                                                              aaaGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATc                                                                                       >  1:2728189/1‑150 (MQ=255)
                                                              aaaGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATc                                                                                       <  1:4419178/150‑1 (MQ=255)
                                                                aGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGc                                                                                     <  1:797690/150‑1 (MQ=255)
                                                                 ggCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGcc                                                                                    >  1:5170344/1‑150 (MQ=255)
                                                                  gCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCt                                                                                   >  1:5308114/1‑150 (MQ=255)
                                                                      ggTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAgg                                                                               >  2:2139193/1‑150 (MQ=255)
                                                                      ggTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAgg                                                                               >  2:636031/1‑150 (MQ=255)
                                                                          aaaaTCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTAt                                                                           >  1:2789260/1‑150 (MQ=255)
                                                                            aaTCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCg                                                                         >  1:2061169/1‑150 (MQ=255)
                                                                             aTCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGt                                                                        >  1:3270666/1‑150 (MQ=255)
                                                                               cATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGttt                                                                      <  1:3607194/150‑1 (MQ=255)
                                                                                 tGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGc                                                                    <  1:1281706/150‑1 (MQ=255)
                                                                                 tGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGc                                                                    >  1:5220034/1‑150 (MQ=255)
                                                                                   gCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGccc                                                                  >  2:1284688/1‑150 (MQ=255)
                                                                                    ccAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGTGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCa                                                                 <  1:3363428/150‑1 (MQ=255)
                                                                                     cAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAg                                                                >  1:3097314/1‑150 (MQ=255)
                                                                                     cAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAg                                                                <  1:4806844/150‑1 (MQ=255)
                                                                                      aaTAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAgg                                                               >  1:3889458/1‑150 (MQ=255)
                                                                                       aTAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGa                                                              >  2:4443630/1‑150 (MQ=255)
                                                                                       aTAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGa                                                              <  1:231549/150‑1 (MQ=255)
                                                                                           aTGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAAcc                                                          >  2:4772036/1‑150 (MQ=255)
                                                                                            tGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCg                                                         <  1:1065042/150‑1 (MQ=255)
                                                                                                 cAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTtgt                                                    >  1:4620152/1‑150 (MQ=255)
                                                                                                  aGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTtgtg                                                   >  1:4990135/1‑150 (MQ=255)
                                                                                                     aTGGTCAGTTTGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTcc                                                <  2:3544700/150‑1 (MQ=255)
                                                                                                     aTGGTCAGTTAGCAATCCCGGGACTTCTTCCTCAACGTCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCTCAACAATCCAGCGTAAATAACGTTTGATGGATTTATCGCATGTGGTTATCGTTTGCCCAGGAAACCGCTTGTGTcc                                                <  2:4388148/150‑1 (MQ=255)
                                                                                                        gTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGt                                             >  1:5098315/1‑150 (MQ=255)
                                                                                                           aGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGtttt                                          >  2:2074458/1‑150 (MQ=255)
                                                                                                            gTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGttttt                                         >  1:1424620/1‑150 (MQ=255)
                                                                                                            gTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGttttt                                         >  2:5038732/1‑150 (MQ=255)
                                                                                                              tAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTa                                       >  2:4126324/1‑150 (MQ=255)
                                                                                                                gCAATCCCGGGACTTCTACCTCACCGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACg                                     <  2:2718535/150‑1 (MQ=255)
                                                                                                                gCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACg                                     <  2:28013/150‑1 (MQ=255)
                                                                                                                     cccGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTAt                                >  1:1150446/1‑150 (MQ=255)
                                                                                                                      ccGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTAtc                               <  2:3250003/150‑1 (MQ=255)
                                                                                                                       cGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTAtct                              <  1:4106159/150‑1 (MQ=255)
                                                                                                                       cGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTGTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTAtct                              <  2:4242846/150‑1 (MQ=255)
                                                                                                                        gggACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTAtctt                             <  1:980856/150‑1 (MQ=255)
                                                                                                                           aCTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTGTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATcttctt                          <  2:4891575/150‑1 (MQ=255)
                                                                                                                             ttCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGa                        >  2:3014448/1‑150 (MQ=255)
                                                                                                                                 aCCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAg                    >  1:165734/1‑150 (MQ=255)
                                                                                                                                 aCCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAg                    >  2:2008749/1‑150 (MQ=255)
                                                                                                                                     cAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATta                <  2:4970636/150‑1 (MQ=255)
                                                                                                                                      aaCGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATtat               <  2:2687636/150‑1 (MQ=255)
                                                                                                                                      aaCGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATtat               <  2:2085790/150‑1 (MQ=255)
                                                                                                                                      aaCGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATtat               <  2:3519931/150‑1 (MQ=255)
                                                                                                                                         ggCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATaaa            <  2:4998298/150‑1 (MQ=255)
                                                                                                                                             tGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATg        <  1:1972080/150‑1 (MQ=255)
                                                                                                                                             tGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATg        <  2:5320843/150‑1 (MQ=255)
                                                                                                                                             tGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATg        <  1:4941408/150‑1 (MQ=255)
                                                                                                                                               aaGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCg      >  1:2168070/1‑150 (MQ=255)
                                                                                                                                                aGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGt     >  1:3016329/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGTGCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGTTTT  >  CP009273/1097795‑1098091

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: