Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,712,562 C→A 63.2% 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,5620CA63.2% 120.4 / 297.1 231D31Y (GAT→TAT) mliCinhibitor of c‑type lysozyme, membrane‑bound; putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (45/40);  new base A (78/68);  total (123/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 233 total aligned reads displayed.

ATTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAGT  >  CP009273/1712417‑1712707
                                                                                                                                                 |                                                                                                                                                 
agtCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCgg                                                                                                                                                 <  2:1404015/146‑1 (MQ=255)
aTTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAgg                                                                                                                                                 >  2:2413828/1‑148 (MQ=255)
 ttCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGt                                                                                                                                                <  1:3359206/148‑1 (MQ=255)
 ttCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGt                                                                                                                                                >  2:1377302/1‑148 (MQ=255)
    cGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTg                                                                                                                                             >  1:3461349/1‑148 (MQ=255)
          tgtAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATgcg                                                                                                                                       >  1:2667076/1‑148 (MQ=255)
               cgcgcgCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCaa                                                                                                                                  >  1:2568336/1‑148 (MQ=255)
                       agagGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGa                                                                                                                          >  1:102055/1‑148 (MQ=255)
                       agagGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGa                                                                                                                          <  1:482794/148‑1 (MQ=255)
                       agagGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGa                                                                                                                          <  2:2964808/148‑1 (MQ=255)
                        gagGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGAt                                                                                                                         <  2:1963091/148‑1 (MQ=255)
                            cTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTaaa                                                                                                                     >  1:1606011/1‑148 (MQ=255)
                             tGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAg                                                                                                                    <  2:769918/148‑1 (MQ=255)
                              gAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAgg                                                                                                                   >  1:2602454/1‑148 (MQ=255)
                                aaTGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGgcg                                                                                                                 <  1:1338683/148‑1 (MQ=255)
                                   gcacTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTa                                                                                                              >  2:3159374/4‑148 (MQ=255)
                                   gCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTa                                                                                                              <  2:2171847/148‑1 (MQ=255)
                                              gaTGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGgagag                                                                                                   >  2:1193080/1‑148 (MQ=255)
                                                 gCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAg                                                                                                >  2:2434521/1‑148 (MQ=255)
                                                 gCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAg                                                                                                <  1:647115/148‑1 (MQ=255)
                                                   agtagtTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGa                                                                                              <  1:1001447/148‑1 (MQ=255)
                                                     tagtTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAATGAcc                                                                                            >  2:3240383/1‑148 (MQ=255)
                                                     tagtTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGAcc                                                                                            >  1:524283/1‑148 (MQ=255)
                                                     tagtTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGAcc                                                                                            >  1:3115476/1‑148 (MQ=255)
                                                      agtTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCg                                                                                           <  2:102055/148‑1 (MQ=255)
                                                      agtTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCg                                                                                           >  1:2283037/1‑148 (MQ=255)
                                                       gtTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCgg                                                                                          >  2:3520180/1‑148 (MQ=255)
                                                       gtTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCgg                                                                                          <  1:2582336/148‑1 (MQ=255)
                                                         tGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCa                                                                                        >  2:796017/1‑148 (MQ=255)
                                                           aTTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCaaa                                                                                      <  1:1850510/148‑1 (MQ=255)
                                                           aTTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCaaa                                                                                      <  2:3461349/148‑1 (MQ=255)
                                                           aTTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCaaa                                                                                      >  1:1292044/1‑148 (MQ=255)
                                                             tATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAAt                                                                                    >  1:2022624/1‑148 (MQ=255)
                                                                cGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCACAATGAt                                                                                 >  1:1201537/1‑148 (MQ=255)
                                                                cGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGAt                                                                                 <  1:1788505/148‑1 (MQ=255)
                                                                 gTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATa                                                                                <  1:560988/148‑1 (MQ=255)
                                                                 gTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATa                                                                                <  1:276486/148‑1 (MQ=255)
                                                                  tAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATaa                                                                               >  2:2571812/1‑148 (MQ=255)
                                                                  tAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATaa                                                                               <  2:877454/148‑1 (MQ=255)
                                                                  tAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATaa                                                                               >  1:1056584/1‑148 (MQ=255)
                                                                  tAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATaa                                                                               >  1:139882/1‑148 (MQ=255)
                                                                   aaaCAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAAt                                                                              >  1:3367218/1‑148 (MQ=255)
                                                                    aaCAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAAtt                                                                             <  2:3379432/148‑1 (MQ=255)
                                                                    aaCAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAAtt                                                                             >  1:3053110/1‑148 (MQ=255)
                                                                      cAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTaa                                                                           >  1:555200/1‑148 (MQ=255)
                                                                      cAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTaa                                                                           <  2:2495474/148‑1 (MQ=255)
                                                                      cAAAACTGACCTCCTGGCGCGGATTATTCAGTGTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTaa                                                                           >  2:2628720/1‑148 (MQ=255)
                                                                       aaaaCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCCACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGATTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATGAAc                                                                          >  1:3540121/1‑148 (MQ=255)
                                                                         aaCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAg                                                                        >  1:1079622/1‑148 (MQ=255)
                                                                         aaCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAg                                                                        >  1:2140581/1‑148 (MQ=255)
                                                                          aCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGt                                                                       >  2:3289908/1‑148 (MQ=255)
                                                                          aCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGt                                                                       <  1:1539203/148‑1 (MQ=255)
                                                                           cTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGtt                                                                      >  1:1245142/1‑148 (MQ=255)
                                                                             gACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGtttt                                                                    >  2:1268319/1‑148 (MQ=255)
                                                                              aCCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGttttt                                                                   >  1:1564721/1‑148 (MQ=255)
                                                                              aCCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGttttt                                                                   <  2:1898375/148‑1 (MQ=255)
                                                                                ctcctGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTc                                                                 <  2:1228866/148‑1 (MQ=255)
                                                                                  cctGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCAt                                                               >  2:595493/1‑148 (MQ=255)
                                                                                  cctGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCAt                                                               <  1:2023287/148‑1 (MQ=255)
                                                                                  cctGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCAt                                                               >  2:2262192/1‑148 (MQ=255)
                                                                                   ctGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATa                                                              <  1:386378/148‑1 (MQ=255)
                                                                                   ctGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATa                                                              <  2:404821/148‑1 (MQ=255)
                                                                                    tGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAg                                                             >  2:1815171/1‑148 (MQ=255)
                                                                                      gcgcGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTc                                                           >  2:728189/1‑148 (MQ=255)
                                                                                      gcgcGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTc                                                           >  2:1309221/1‑148 (MQ=255)
                                                                                           gATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATc                                                      <  1:85503/148‑1 (MQ=255)
                                                                                            attattCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATcc                                                     <  2:1725442/148‑1 (MQ=255)
                                                                                               attCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGa                                                  <  2:3177276/148‑1 (MQ=255)
                                                                                               attCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGa                                                  <  2:1056637/148‑1 (MQ=255)
                                                                                               attCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGa                                                  <  2:1568234/148‑1 (MQ=255)
                                                                                                      ttGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACttt                                           >  2:1071688/1‑148 (MQ=255)
                                                                                                        gACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTcc                                         <  2:2039707/148‑1 (MQ=255)
                                                                                                        gACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTcc                                         <  1:216282/148‑1 (MQ=255)
                                                                                                        gACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTcc                                         >  2:629311/1‑148 (MQ=255)
                                                                                                          ccGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTg                                       <  1:2622107/148‑1 (MQ=255)
                                                                                                             tCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTc                                    <  2:524283/148‑1 (MQ=255)
                                                                                                             tCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTc                                    <  2:3115476/148‑1 (MQ=255)
                                                                                                                aCGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTgg                                 >  1:561788/1‑148 (MQ=255)
                                                                                                                   gTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAgg                              >  2:3133888/1‑148 (MQ=255)
                                                                                                                    tttttCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAggg                             >  1:2351096/1‑148 (MQ=255)
                                                                                                                     ttttCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGc                            >  2:2608156/1‑148 (MQ=255)
                                                                                                                     ttttCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGc                            >  1:2770756/1‑148 (MQ=255)
                                                                                                                     ttttCATCACACTGGTATCCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGc                            >  1:552416/1‑148 (MQ=255)
                                                                                                                      tttCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCa                           <  2:2946738/148‑1 (MQ=255)
                                                                                                                       ttCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCaa                          <  2:1739274/148‑1 (MQ=255)
                                                                                                                        tcatcaCACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAAt                         <  2:1623810/148‑1 (MQ=255)
                                                                                                                         catcaCACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATa                        >  1:1059515/1‑148 (MQ=255)
                                                                                                                          atcaCACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAAACAGCTGATTAAATGCGCTACAGCAGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTATCAGAGTAATATCCCGAAGAATTACCAGGTATGGAGGGCAATAc                       >  1:716471/1‑148 (MQ=255)
                                                                                                                           tcaCACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACg                      >  2:188738/1‑148 (MQ=255)
                                                                                                                           tcaCACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACg                      >  2:1827843/1‑148 (MQ=255)
                                                                                                                            cacaCTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGc                     >  1:3211322/1‑148 (MQ=255)
                                                                                                                                 tGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGccctcc                <  1:2425354/148‑1 (MQ=255)
                                                                                                                                 tGGTAGTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGccctcc                <  2:1759660/148‑1 (MQ=255)
                                                                                                                                  ggTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGccctccc               <  2:2786395/148‑1 (MQ=255)
                                                                                                                                     aTTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTaa            <  1:2128610/148‑1 (MQ=255)
                                                                                                                                       tCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACg          <  1:3159374/148‑1 (MQ=255)
                                                                                                                                           gCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTcc      >  2:1489364/1‑148 (MQ=255)
                                                                                                                                            cGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCa     <  2:1844901/148‑1 (MQ=255)
                                                                                                                                              tataGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAg   >  2:657752/1‑148 (MQ=255)
                                                                                                                                               aTCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAgt  <  1:1676326/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
ATTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAGT  >  CP009273/1712417‑1712707

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: