Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,386,303 A→T 14.3% I1105I (ATT→ATA yhdP ← DUF3971‑AsmA2 domains protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,386,3030AT14.3% 1031.8 / 145.9 350I1105I (ATT→ATAyhdPDUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand):  ref base A (169/131);  new base T (25/25);  total (194/156)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.31e-01

Only 100 of 366 total aligned reads displayed.

TGGAGTCAAAATAGAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAATTTTCTGCCCGCGCAGTTTTCC  >  CP009273/3386155‑3386452
                                                                                                                                                    |                                                                                                                                                     
  gAGTCAAAATAGAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTatgaatg                                                                                                                                                    <  2:1512895/150‑1 (MQ=255)
    gTCAAAATAGAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAatgat                                                                                                                                                  <  2:4346601/150‑1 (MQ=255)
     tCAAAATAGAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAatgatg                                                                                                                                                 >  2:2080423/1‑150 (MQ=255)
     tCAAAATAGAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAatgatg                                                                                                                                                 <  1:2602181/150‑1 (MQ=255)
      cAAAATAGAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGc                                                                                                                                                <  1:2267163/150‑1 (MQ=255)
       aaaaTAGAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGcc                                                                                                                                               >  2:91236/1‑150 (MQ=255)
          aTAGAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCAtt                                                                                                                                            <  1:443411/150‑1 (MQ=255)
            aGAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCa                                                                                                                                          >  1:1012266/1‑150 (MQ=255)
            aGAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCa                                                                                                                                          <  2:1075799/150‑1 (MQ=255)
             gAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTCTGAATGATGCCATACaa                                                                                                                                         >  2:2963351/1‑150 (MQ=255)
             gAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGAATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCaa                                                                                                                                         <  1:5133155/150‑1 (MQ=255)
              aaCCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCAAc                                                                                                                                        <  2:4898785/150‑1 (MQ=255)
               aCCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACg                                                                                                                                       <  1:1615173/150‑1 (MQ=255)
                  cTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATGCAACGTCg                                                                                                                                    >  2:3609413/1‑150 (MQ=255)
                   ttCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGc                                                                                                                                   >  2:5008123/1‑150 (MQ=255)
                     cGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCAACGTCGCtt                                                                                                                                 >  2:5021108/1‑150 (MQ=255)
                       ccAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCa                                                                                                                               >  1:4104769/1‑150 (MQ=255)
                        cAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCAACGTCGCTTCAt                                                                                                                              >  1:675/1‑150 (MQ=255)
                        cAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCAt                                                                                                                              >  1:3764667/1‑150 (MQ=255)
                        cAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCAt                                                                                                                              >  2:1774925/1‑150 (MQ=255)
                        cAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCAt                                                                                                                              >  2:4700754/1‑150 (MQ=255)
                          aaaGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCt                                                                                                                            <  1:4336099/150‑1 (MQ=255)
                          aaaGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCt                                                                                                                            <  1:1033611/150‑1 (MQ=255)
                                 ctctGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGcc                                                                                                                     <  1:3813968/150‑1 (MQ=255)
                                  tctGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCAACGTCGCTTCATCTGGCTGCCa                                                                                                                    <  2:1731175/150‑1 (MQ=255)
                                  tctGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCa                                                                                                                    >  1:3997078/1‑150 (MQ=255)
                                    tGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAgg                                                                                                                  <  2:858311/150‑1 (MQ=255)
                                     gAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAggg                                                                                                                 >  2:4842921/1‑150 (MQ=255)
                                      aaaTCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGt                                                                                                                >  2:1364804/1‑150 (MQ=255)
                                      aaaTCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGt                                                                                                                <  1:1107702/150‑1 (MQ=255)
                                       aaTCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTg                                                                                                               >  1:2134274/1‑150 (MQ=255)
                                          cAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCtt                                                                                                            <  2:2355911/150‑1 (MQ=255)
                                           aaaaCGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCttt                                                                                                           <  1:3030382/150‑1 (MQ=255)
                                           aaaaCGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGACAGGTTGc                                                                                                              >  2:4579595/1‑147 (MQ=255)
                                             aaCGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTgc                                                                                                         <  1:1501566/150‑1 (MQ=255)
                                               cGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTgcgc                                                                                                       >  2:2783180/1‑150 (MQ=255)
                                                gCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGcc                                                                                                      >  2:2614015/1‑150 (MQ=255)
                                                gCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGcc                                                                                                      <  2:3350470/150‑1 (MQ=255)
                                                gCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGcc                                                                                                      >  2:250648/1‑150 (MQ=255)
                                                  aGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAg                                                                                                    <  2:4961828/150‑1 (MQ=255)
                                                  aGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAg                                                                                                    <  2:2310913/150‑1 (MQ=255)
                                                  aGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAg                                                                                                    >  1:4584320/1‑150 (MQ=255)
                                                     ttACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAgtgt                                                                                                 <  2:1148220/150‑1 (MQ=255)
                                                        cGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTaaa                                                                                              <  2:748589/150‑1 (MQ=255)
                                                        cGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATGGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTaaa                                                                                              <  2:1473135/150‑1 (MQ=255)
                                                         gCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTGTGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAAt                                                                                             >  1:4939867/1‑150 (MQ=255)
                                                         gCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAAt                                                                                             <  2:3512674/150‑1 (MQ=255)
                                                         gCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAAt                                                                                             <  2:3816940/150‑1 (MQ=255)
                                                          cATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATc                                                                                            >  2:3944222/1‑150 (MQ=255)
                                                           aTCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCg                                                                                           >  1:1918469/1‑150 (MQ=255)
                                                              aGGGCATCTACGCTCAATAAGCGCAGCACTTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTaa                                                                                        >  2:2718171/1‑150 (MQ=255)
                                                                      tACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACAtt                                                                                <  2:3272800/150‑1 (MQ=255)
                                                                      tACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACAtt                                                                                <  2:2013973/150‑1 (MQ=255)
                                                                         gCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTaaa                                                                             >  2:3586465/1‑150 (MQ=255)
                                                                          cTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAAt                                                                            <  1:1172284/150‑1 (MQ=255)
                                                                           tCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATg                                                                           <  1:840544/150‑1 (MQ=255)
                                                                            cAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATga                                                                          >  1:508903/1‑150 (MQ=255)
                                                                             aaTAAGCGCAGCAATTGCCCGGCATGTCCGGGCTTGATGTCGGTAATTTCGCCTTTACCCAGTTGAGTATGACTGATGCCATTCAACGTCGCGTCATCTGGCTGCCCGGGTGCTTTGCGCCAGTGTAAATCGTCAGCCACAGTAAATgac                                                                         >  1:880356/1‑150 (MQ=255)
                                                                              aTAAGCGCAGCAAGTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATgacg                                                                        <  2:3285038/150‑1 (MQ=255)
                                                                               tAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATgacga                                                                       <  1:1278661/150‑1 (MQ=255)
                                                                                 aGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACt                                                                     <  1:3704163/150‑1 (MQ=255)
                                                                                  gcgcAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTg                                                                    <  2:1850353/150‑1 (MQ=255)
                                                                                   cgcAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTgg                                                                   <  2:4274951/150‑1 (MQ=255)
                                                                                    gcagcaATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGc                                                                  <  2:2881598/150‑1 (MQ=255)
                                                                                    gcagcaATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGc                                                                  <  2:290115/150‑1 (MQ=255)
                                                                                        caATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTAt                                                              <  1:1424531/150‑1 (MQ=255)
                                                                                         aaTTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATg                                                             >  1:1999425/1‑150 (MQ=255)
                                                                                         aaTTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATg                                                             <  2:1301393/150‑1 (MQ=255)
                                                                                         aaTTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATg                                                             >  2:1204344/1‑150 (MQ=255)
                                                                                         aaTTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATg                                                             >  1:2391173/1‑150 (MQ=255)
                                                                                            tGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGc                                                          <  1:1448988/150‑1 (MQ=255)
                                                                                              cccGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGcgt                                                        >  1:815795/1‑150 (MQ=255)
                                                                                               ccGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGcgtc                                                       <  1:3322096/150‑1 (MQ=255)
                                                                                               ccGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGcgtc                                                       <  2:3114020/150‑1 (MQ=255)
                                                                                                cGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGcgtcg                                                      >  1:2104052/1‑150 (MQ=255)
                                                                                                  gCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTg                                                    <  1:2228924/150‑1 (MQ=255)
                                                                                                   cATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGa                                                   >  1:567152/1‑150 (MQ=255)
                                                                                                   cATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGa                                                   >  1:3153043/1‑150 (MQ=255)
                                                                                                   cATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGa                                                   >  1:3269448/1‑150 (MQ=255)
                                                                                                    aTGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGTATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGac                                                  >  2:3733017/1‑150 (MQ=255)
                                                                                                     tGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGaca                                                 >  2:3075231/1‑150 (MQ=255)
                                                                                                      gTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGacac                                                <  2:3013845/150‑1 (MQ=255)
                                                                                                         cGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCaa                                             >  1:277435/1‑150 (MQ=255)
                                                                                                           gTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCaaaa                                           >  1:4139545/1‑150 (MQ=255)
                                                                                                            tATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCaaaaa                                          >  2:4090130/1‑150 (MQ=255)
                                                                                                             aTTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCaaaaaa                                         >  2:1658937/1‑150 (MQ=255)
                                                                                                              ttGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAAt                                        >  1:1933957/1‑150 (MQ=255)
                                                                                                                gATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTc                                      >  1:2908976/1‑150 (MQ=255)
                                                                                                                gATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTc                                      >  1:4535892/1‑150 (MQ=255)
                                                                                                                   ttCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTgc                                   <  2:1028273/150‑1 (MQ=255)
                                                                                                                             tCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCa                         <  1:3591364/150‑1 (MQ=255)
                                                                                                                                ccTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATTCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAAtt                      >  1:4852829/1‑150 (MQ=255)
                                                                                                                                ccTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCCAtt                      >  1:4480706/1‑150 (MQ=255)
                                                                                                                                ccTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAAtt                      <  2:4474209/150‑1 (MQ=255)
                                                                                                                                ccTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAAtt                      <  1:342735/150‑1 (MQ=255)
                                                                                                                                      cccAGTTGAGTATGTATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAATTTTCTGc                >  1:2839493/1‑150 (MQ=255)
                                                                                                                                         aGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAATTTTCTGCCcg             <  1:1760524/150‑1 (MQ=255)
                                                                                                                                             gAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAATTTTCTGCCCGCGCa         >  2:5047856/1‑150 (MQ=255)
                                                                                                                                                tATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAATTTTCTGCCCGCGCAGtt      >  2:2231425/1‑150 (MQ=255)
                                                                                                                                                tATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAATTTTCTGCCCGCGCAGtt      >  1:1640446/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGGAGTCAAAATAGAACCCTTCGCCAAAAGTGTCTCTGAAATCAAAACGCAGCTTACGCATCAGGGCATCTACGCTCAATAAGCGCAGCAATTGCCCGGCATGTCCGGTATTGATTTCGGTAATTTCGCCTTTACCCAGTTGAGTATGAATGATGCCATTCAACGTCGCTTCATCTGGCTGCCAGGGTGCTTTGCGCCAGTGTAAATCGTAATCCACATTAAATGACGACTGGCGTATGGGCGTCGTGACACCAAAAAATTCTGCGGCGGCATCAATTTTCTGCCCGCGCAGTTTTCC  >  CP009273/3386155‑3386452

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: