Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,385,315 C→A 5.3% A266S (GCA→TCA)  ycjY ← S9 homolog non‑peptidase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,385,3150CA 5.3% 1251.3 / 35.8 319A266S (GCA→TCA) ycjYS9 homolog non‑peptidase family protein
Reads supporting (aligned to +/- strand):  ref base C (167/135);  new base A (7/10);  total (174/145)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 319 total aligned reads displayed.

TTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCTGCGCAAAGTAGCGTAACCCGG  >  CP009273/1385168‑1385462
                                                                                                                                                   |                                                                                                                                                   
ttGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCt                                                                                                                                                    <  2:657021/149‑1 (MQ=255)
ttGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCt                                                                                                                                                    <  2:1122410/149‑1 (MQ=255)
 tGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGATc                                                                                                                                                   <  2:124947/149‑1 (MQ=255)
   cGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGa                                                                                                                                                 <  2:5004688/149‑1 (MQ=255)
     gATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATc                                                                                                                                               >  2:2229484/1‑149 (MQ=255)
     gATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATc                                                                                                                                               >  1:1072154/1‑149 (MQ=255)
     gATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATc                                                                                                                                               <  2:2241702/149‑1 (MQ=255)
        ttGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTa                                                                                                                                            <  2:2832722/149‑1 (MQ=255)
          gAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACa                                                                                                                                          >  1:705292/1‑149 (MQ=255)
            atatGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAgg                                                                                                                                        <  1:4219730/149‑1 (MQ=255)
               tGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGgtcg                                                                                                                                     >  1:3786410/1‑149 (MQ=255)
                gTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGgtcgt                                                                                                                                    >  2:3508703/1‑149 (MQ=255)
                   tCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCAc                                                                                                                                 <  1:190802/149‑1 (MQ=255)
                     aTTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTc                                                                                                                               >  2:1300049/1‑149 (MQ=255)
                     aTTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTc                                                                                                                               >  1:3808142/1‑149 (MQ=255)
                     aTTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTc                                                                                                                               >  1:4285010/1‑149 (MQ=255)
                     aTTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTACAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTc                                                                                                                               >  2:4583329/1‑149 (MQ=255)
                          agagTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCa                                                                                                                          <  1:1007544/149‑1 (MQ=255)
                             gTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCAttt                                                                                                                       >  1:2410626/1‑149 (MQ=255)
                             gTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCAttt                                                                                                                       >  1:369237/1‑149 (MQ=255)
                              tttCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTg                                                                                                                      >  2:3854254/1‑149 (MQ=255)
                               ttCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGc                                                                                                                     <  1:3910278/149‑1 (MQ=255)
                               ttCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGc                                                                                                                     <  1:899098/149‑1 (MQ=255)
                                tCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCt                                                                                                                    >  1:4181235/1‑149 (MQ=255)
                                 ccTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTc                                                                                                                   <  2:3736851/149‑1 (MQ=255)
                                  cTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTcc                                                                                                                  >  2:3889481/1‑149 (MQ=255)
                                   tCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTccc                                                                                                                 <  2:4758604/149‑1 (MQ=255)
                                      aaGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCcctgc                                                                                                              >  1:1538214/1‑149 (MQ=255)
                                         aaCGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGATCGATCGTACAGGTCGTCACTCATCCATTTGCTCcctgcctg                                                                                                           <  2:1370339/149‑1 (MQ=255)
                                            ggAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCt                                                                                                        >  2:3410042/1‑149 (MQ=255)
                                            ggAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCt                                                                                                        <  1:2696324/149‑1 (MQ=255)
                                                   aTACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGcc                                                                                                 >  1:742783/1‑149 (MQ=255)
                                                   aTACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGcc                                                                                                 <  2:1257113/149‑1 (MQ=255)
                                                    tACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCca                                                                                                >  1:3468897/1‑149 (MQ=255)
                                                       tGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACaa                                                                                             <  1:1019588/149‑1 (MQ=255)
                                                         aaaTTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATc                                                                                           >  2:4037045/1‑149 (MQ=255)
                                                              gCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGt                                                                                      >  1:597101/1‑149 (MQ=255)
                                                                   aGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGTTGCCCGCCACAATCTGCGTCGGCt                                                                                 >  2:4499245/1‑149 (MQ=255)
                                                                   aGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCt                                                                                 >  2:4662459/1‑149 (MQ=255)
                                                                    gCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTg                                                                                >  2:72313/1‑149 (MQ=255)
                                                                    gCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTg                                                                                <  1:3362929/149‑1 (MQ=255)
                                                                         aTAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCa                                                                           <  1:2224973/149‑1 (MQ=255)
                                                                         aTAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCa                                                                           <  2:842631/149‑1 (MQ=255)
                                                                         aTAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCa                                                                           <  1:4818802/149‑1 (MQ=255)
                                                                            ggCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGt                                                                        >  1:961251/1‑149 (MQ=255)
                                                                             gCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTa                                                                       >  2:5115626/1‑149 (MQ=255)
                                                                                 tGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACt                                                                   >  2:3239916/1‑149 (MQ=255)
                                                                                    cGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTcc                                                                <  2:327188/149‑1 (MQ=255)
                                                                                     gTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCg                                                               >  1:2827280/1‑149 (MQ=255)
                                                                                          aaaGATCCATATGGTTTGCGCCTTCAATAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCata                                                          <  2:2410176/149‑1 (MQ=255)
                                                                                          aaaGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCata                                                          <  2:5011648/149‑1 (MQ=255)
                                                                                             gATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATgg                                                       >  1:1162950/1‑149 (MQ=255)
                                                                                             gATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATgg                                                       >  1:251336/1‑149 (MQ=255)
                                                                                               tCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTa                                                     >  1:4099413/1‑149 (MQ=255)
                                                                                                ccATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTTGGCTGAGTCAGGTACACTTCCGCCATATGGTaa                                                    <  1:242119/149‑1 (MQ=255)
                                                                                                ccATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTaa                                                    <  2:2594360/149‑1 (MQ=255)
                                                                                                  atatGGTGTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGc                                                  <  1:5164896/149‑1 (MQ=255)
                                                                                                   tatGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGca                                                 >  2:3067085/1‑149 (MQ=255)
                                                                                                       gTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGcatcat                                             <  2:1075822/149‑1 (MQ=255)
                                                                                                         ttGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAg                                           <  1:4522291/149‑1 (MQ=255)
                                                                                                         ttGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAg                                           <  2:692095/149‑1 (MQ=255)
                                                                                                             gcCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGtaa                                       >  2:251055/1‑149 (MQ=255)
                                                                                                              ccTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGtaat                                      >  2:4325909/1‑149 (MQ=255)
                                                                                                              ccTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGtaat                                      >  2:1456412/1‑149 (MQ=255)
                                                                                                              ccTTCAACAATGTGAGAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGtaat                                      <  1:3696042/149‑1 (MQ=255)
                                                                                                               cTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGtaata                                     <  2:2331648/149‑1 (MQ=255)
                                                                                                               cTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGtaata                                     >  2:4331092/1‑149 (MQ=255)
                                                                                                               cTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGATCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGtaata                                     >  1:1829352/1‑149 (MQ=255)
                                                                                                                 tCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTaataat                                   <  2:2796340/149‑1 (MQ=255)
                                                                                                                   aacaaTGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCt                                 <  2:2098071/149‑1 (MQ=255)
                                                                                                                     caaTGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTgg                               <  1:1326947/149‑1 (MQ=255)
                                                                                                                     caaTGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTgg                               <  2:365654/149‑1 (MQ=255)
                                                                                                                      aaTGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGt                              <  1:2465815/149‑1 (MQ=255)
                                                                                                                       aTGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGtt                             >  2:868167/1‑149 (MQ=255)
                                                                                                                       aTGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGtt                             <  2:365020/149‑1 (MQ=255)
                                                                                                                         gtgATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTaa                           >  2:3614881/1‑149 (MQ=255)
                                                                                                                           gATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAgg                         <  2:938863/149‑1 (MQ=255)
                                                                                                                            aTAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGc                        >  1:3257391/1‑149 (MQ=255)
                                                                                                                             tAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCt                       >  1:4553348/1‑149 (MQ=255)
                                                                                                                             tAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCt                       <  2:533242/149‑1 (MQ=255)
                                                                                                                               gCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCTgc                     >  1:4734611/1‑149 (MQ=255)
                                                                                                                                 gTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCTgcgc                   <  1:1767161/149‑1 (MQ=255)
                                                                                                                                  tttATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCTGCGCa                  >  1:1226197/1‑149 (MQ=255)
                                                                                                                                    tATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCTGCGCaaa                <  2:4181235/149‑1 (MQ=255)
                                                                                                                                     aTCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCTGCGCAAAg               >  2:1540260/1‑149 (MQ=255)
                                                                                                                                      tCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCTGCGCAAAGt              >  1:1712026/1‑149 (MQ=255)
                                                                                                                                        cTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCTGCGCAAAGTAg            <  2:568625/149‑1 (MQ=255)
                                                                                                                                        cTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCTGCGCAAAGTAg            <  2:3313021/149‑1 (MQ=255)
                                                                                                                                        cTGGCTTGATGCTAGATCGTACAGGTCGTCACTCATCCATTTGATCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGTCTGCGCAAAGTAg            >  1:2300107/1‑149 (MQ=255)
                                                                                                                                         tGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCTGCGCAAAGTAGc           >  1:4015926/1‑149 (MQ=255)
                                                                                                                                          ggCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCTGCGCAAAGTAGCg          >  2:1846204/1‑149 (MQ=255)
                                                                                                                                          ggCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCTGCGCAAAGTAGCg          >  2:1729444/1‑149 (MQ=255)
                                                                                                                                          ggCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCAGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCTGCGCAAAGTAGCg          >  1:3112365/1‑149 (MQ=255)
                                                                                                                                             ttGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCTGCGCAAAGTAGCGTaa       >  1:74407/1‑149 (MQ=255)
                                                                                                                                               gATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCTGCGCAAAGTAGCGTAAcc     >  2:3950362/1‑149 (MQ=255)
                                                                                                                                               gATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCGGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCTGCGCAAAGTAGCGTAAcc     <  1:4799637/149‑1 (MQ=255)
                                                                                                                                                aTGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCTGCGCAAAGTAGCGTAAccc    <  1:5035619/149‑1 (MQ=255)
                                                                                                                                                  gCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCTGCGCAAAGTAGCGTAACCCgg  >  1:3638279/1‑149 (MQ=255)
                                                                                                                                                  gCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCTGCGCAAAGTAGCGTAACCCgg  <  1:2938182/149‑1 (MQ=255)
                                                                                                                                                  gCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCTGCGCAAAGTAGCGTAACCCgg  <  2:4242245/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGCCTTCAACAATGTGATAGCGTTTATCCTGGCTTGATGCTCGATCGTACAGGTCGTCACTCATCCATTTGCTCCCTGCCTGGCTGCCCGCCACAATCTGCGTCGGCTGAGTCAGGTACACTTCCGCCATATGGTAAGCATCATAGGTAATAATCTGGTTAAGGCTGCGCAAAGTAGCGTAACCCGG  >  CP009273/1385168‑1385462

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