Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,712,562 C→A 44.1% D31Y (GAT→TAT)  mliC ← inhibitor of c‑type lysozyme, membrane‑bound; putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,712,5620CA44.1% 189.5 / 248.5 215D31Y (GAT→TAT) mliCinhibitor of c‑type lysozyme, membrane‑bound; putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (53/67);  new base A (54/41);  total (107/108)
Fisher's exact test for biased strand distribution p-value = 7.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 216 total aligned reads displayed.

ATTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAGTGT  >  CP009273/1712417‑1712709
                                                                                                                                                 |                                                                                                                                                   
aTTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCgg                                                                                                                                                   >  2:2727418/1‑148 (MQ=255)
aTTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAgg                                                                                                                                                   >  2:838835/1‑148 (MQ=255)
aTTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAgg                                                                                                                                                   <  2:3124955/148‑1 (MQ=255)
 ttCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGt                                                                                                                                                  >  2:2021681/1‑148 (MQ=255)
  tCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGtt                                                                                                                                                 >  2:2884127/1‑148 (MQ=255)
  tCCGTCAGTGTAACGCGCGCCAGAGGCTGAAAGGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGtt                                                                                                                                                 <  2:1977169/148‑1 (MQ=255)
     gTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGc                                                                                                                                              <  1:1504150/148‑1 (MQ=255)
      tCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCa                                                                                                                                             <  2:213676/148‑1 (MQ=255)
       cAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCAt                                                                                                                                            >  2:1871668/1‑148 (MQ=255)
       cAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCAt                                                                                                                                            >  2:3582901/1‑148 (MQ=255)
               cgcgcgCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCaa                                                                                                                                    >  2:3159561/1‑148 (MQ=255)
               cgcgcgCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCaa                                                                                                                                    >  1:3325536/1‑148 (MQ=255)
               cgcgcgCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCaa                                                                                                                                    >  1:456368/1‑148 (MQ=255)
               cgcgcgCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCaa                                                                                                                                    <  1:734811/148‑1 (MQ=255)
                       agagGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGa                                                                                                                            <  1:835068/148‑1 (MQ=255)
                        gagGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGAt                                                                                                                           >  1:3523016/1‑148 (MQ=255)
                         agGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGAtt                                                                                                                          >  2:779514/1‑148 (MQ=255)
                          ggCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTa                                                                                                                         >  1:2808573/1‑148 (MQ=255)
                           gCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTaa                                                                                                                        >  2:2903918/1‑148 (MQ=255)
                            cTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTaaa                                                                                                                       <  2:2099336/148‑1 (MQ=255)
                            cTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTaaa                                                                                                                       <  2:157491/148‑1 (MQ=255)
                            cTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTaaa                                                                                                                       >  2:762390/1‑148 (MQ=255)
                              gAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAgg                                                                                                                     <  2:2046507/148‑1 (MQ=255)
                              gAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAgg                                                                                                                     <  1:50916/148‑1 (MQ=255)
                              gAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAgg                                                                                                                     <  1:113597/148‑1 (MQ=255)
                              gAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAgg                                                                                                                     <  2:544441/148‑1 (MQ=255)
                               aaaTGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGgc                                                                                                                    >  2:1091347/1‑148 (MQ=255)
                                 aTGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGgcgc                                                                                                                  >  2:159968/1‑148 (MQ=255)
                                  tGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGGCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCt                                                                                                                 >  1:973685/1‑148 (MQ=255)
                                   gCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTa                                                                                                                >  2:3179419/1‑148 (MQ=255)
                                        gTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGcc                                                                                                           >  1:103362/1‑148 (MQ=255)
                                          ttGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCgg                                                                                                         >  2:2267603/1‑148 (MQ=255)
                                               aTGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGgagaga                                                                                                    <  2:1444364/148‑1 (MQ=255)
                                               aTGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGgagaga                                                                                                    >  2:691230/1‑148 (MQ=255)
                                               aTGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGgagaga                                                                                                    <  1:1514578/148‑1 (MQ=255)
                                                tGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGaa                                                                                                   >  1:3761640/1‑148 (MQ=255)
                                                 gCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAg                                                                                                  <  2:801066/148‑1 (MQ=255)
                                                  cAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAgg                                                                                                 <  1:954445/148‑1 (MQ=255)
                                                   agtagtTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGa                                                                                                <  2:631788/148‑1 (MQ=255)
                                                   agtagtTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGa                                                                                                <  1:2515181/148‑1 (MQ=255)
                                                     tagtTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGAcc                                                                                              >  1:1917492/1‑148 (MQ=255)
                                                           aTTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCaaa                                                                                        <  2:1096812/148‑1 (MQ=255)
                                                             tATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAAt                                                                                      >  2:3645058/1‑148 (MQ=255)
                                                             tATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAAt                                                                                      >  2:1358203/1‑148 (MQ=255)
                                                             tATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTACAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAAt                                                                                      >  2:2685019/1‑148 (MQ=255)
                                                               tCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGa                                                                                    <  1:930722/148‑1 (MQ=255)
                                                               tCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGa                                                                                    <  1:1241601/148‑1 (MQ=255)
                                                               tCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTAATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGa                                                                                    >  1:3605265/1‑148 (MQ=255)
                                                                cGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGAt                                                                                   <  1:3074723/148‑1 (MQ=255)
                                                                cGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGAt                                                                                   >  1:2609950/1‑148 (MQ=255)
                                                                  tAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATaa                                                                                 <  1:3677817/148‑1 (MQ=255)
                                                                  tAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATaa                                                                                 <  2:37443/148‑1 (MQ=255)
                                                                  tAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATaa                                                                                 >  2:2433087/1‑148 (MQ=255)
                                                                  tAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGACGGAGAGAAGGACCGGCAAAATGATaa                                                                                 >  2:2090502/1‑148 (MQ=255)
                                                                     aCAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTa                                                                              >  2:1206526/1‑148 (MQ=255)
                                                                         aaCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAg                                                                          <  1:2445165/148‑1 (MQ=255)
                                                                          aCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGt                                                                         >  1:162696/1‑148 (MQ=255)
                                                                          aCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGt                                                                         <  1:3582901/148‑1 (MQ=255)
                                                                          aCTGACCTCCTGGCGCGGATCATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGt                                                                         <  2:3456445/148‑1 (MQ=255)
                                                                             gACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGtttt                                                                      <  2:2769524/148‑1 (MQ=255)
                                                                               cctcctGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGtttttt                                                                    <  1:3159561/148‑1 (MQ=255)
                                                                               cctcctGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGATTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGtttttt                                                                    >  2:1635307/1‑148 (MQ=255)
                                                                                ctcctGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTc                                                                   <  1:1353184/148‑1 (MQ=255)
                                                                                  cctGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCAt                                                                 <  2:2348937/148‑1 (MQ=255)
                                                                                   ctGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATa                                                                <  2:1804251/148‑1 (MQ=255)
                                                                                    tGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAg                                                               <  2:1275311/148‑1 (MQ=255)
                                                                                     ggCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGt                                                              <  2:2535377/148‑1 (MQ=255)
                                                                                       cgcgGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCa                                                            <  1:3626811/148‑1 (MQ=255)
                                                                                          ggATTATTCTGTTTGACCGTCAAAGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCTACCAGCTGATTAAAGGCGCTACAGCCGGAGAGTAGGACCGGCAAAATGATAAGTAACAGTTTTTTCATAGTCatat                                                         <  2:3599341/148‑1 (MQ=255)
                                                                                          ggATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCatat                                                         <  1:159968/148‑1 (MQ=255)
                                                                                           gATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATc                                                        <  2:2173678/148‑1 (MQ=255)
                                                                                                 tCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAg                                                  >  1:2126230/1‑148 (MQ=255)
                                                                                                    gTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACt                                               >  1:2708164/1‑148 (MQ=255)
                                                                                                       tGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTc                                            <  1:1651178/148‑1 (MQ=255)
                                                                                                        gACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTcc                                           <  1:655834/148‑1 (MQ=255)
                                                                                                         aCCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCt                                          <  1:779514/148‑1 (MQ=255)
                                                                                                         aCCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCt                                          <  1:2884127/148‑1 (MQ=255)
                                                                                                           cGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTgg                                        <  2:456906/148‑1 (MQ=255)
                                                                                                              cAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCt                                     >  2:3044390/1‑148 (MQ=255)
                                                                                                                aCGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTgg                                   >  1:3580810/1‑148 (MQ=255)
                                                                                                                      tttCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCa                             >  2:3009499/1‑148 (MQ=255)
                                                                                                                      tttCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCa                             >  2:1487502/1‑148 (MQ=255)
                                                                                                                      tttCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCa                             >  2:2832165/1‑148 (MQ=255)
                                                                                                                      tttCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCa                             >  2:1104635/1‑148 (MQ=255)
                                                                                                                       ttCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCaa                            <  1:1143555/148‑1 (MQ=255)
                                                                                                                           tcaCACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACg                        <  1:688626/148‑1 (MQ=255)
                                                                                                                               acTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGccct                    <  2:3384045/148‑1 (MQ=255)
                                                                                                                                cTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGccctc                   <  2:3306567/148‑1 (MQ=255)
                                                                                                                                 tGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGccctcc                  <  1:2949292/148‑1 (MQ=255)
                                                                                                                                 tGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGccctcc                  <  1:691230/148‑1 (MQ=255)
                                                                                                                                 tGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGccctcc                  <  1:762886/148‑1 (MQ=255)
                                                                                                                                    tATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAAAAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTa               <  1:2262599/148‑1 (MQ=255)
                                                                                                                                    tATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTa               >  2:850656/1‑148 (MQ=255)
                                                                                                                                     aTTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTaa              >  2:1305719/1‑148 (MQ=255)
                                                                                                                                     aTTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCCTAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTaa              <  1:1970214/148‑1 (MQ=255)
                                                                                                                                      ttCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAAc             >  1:2450841/1‑148 (MQ=255)
                                                                                                                                      ttCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAAc             >  1:530583/1‑148 (MQ=255)
                                                                                                                                           gCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTcc        >  2:1394819/1‑148 (MQ=255)
                                                                                                                                            cGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCa       <  1:2792401/148‑1 (MQ=255)
                                                                                                                                                tCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAgtg   <  1:2406865/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ATTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAGTGT  >  CP009273/1712417‑1712709

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: