Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,712,562 C→A 46.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,5620CA46.1% 134.2 / 239.6 195D31Y (GAT→TAT) mliCinhibitor of c‑type lysozyme, membrane‑bound; putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (51/54);  new base A (46/44);  total (97/98)
Fisher's exact test for biased strand distribution p-value = 7.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 200 total aligned reads displayed.

AGATTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAGTGT  >  CP009273/1712415‑1712709
                                                                                                                                                   |                                                                                                                                                   
aGATTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATc                                                                                                                                                     <  2:1369256/148‑1 (MQ=255)
  aTTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCgg                                                                                                                                                   <  2:938364/148‑1 (MQ=255)
  aTTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAgg                                                                                                                                                   >  2:3222062/1‑148 (MQ=255)
   ttCCGTCAGTGTAACGCGCGCCAGAGGCTGAGATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGt                                                                                                                                                  >  1:1644758/1‑148 (MQ=255)
   ttCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGt                                                                                                                                                  <  2:511030/148‑1 (MQ=255)
      cGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTAAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTg                                                                                                                                               >  2:3040070/1‑148 (MQ=255)
        tCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCATCGTATAGGTTTGCa                                                                                                                                             >  1:887018/1‑148 (MQ=255)
         cAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCAt                                                                                                                                            <  1:357621/148‑1 (MQ=255)
         cAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCAt                                                                                                                                            <  1:595856/148‑1 (MQ=255)
            tgtAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATgcg                                                                                                                                         <  1:2322309/148‑1 (MQ=255)
              tAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCATCGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCt                                                                                                                                       <  2:3016289/148‑1 (MQ=255)
                aCGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCa                                                                                                                                     >  1:1124521/1‑148 (MQ=255)
                 cgcgcgCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCaa                                                                                                                                    >  1:1120414/1‑148 (MQ=255)
                 cgcgcgCCAGAGGCTGAAATGCCCTGTTTGAGATGAAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCaa                                                                                                                                    <  2:2169120/148‑1 (MQ=255)
                          gagGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGAt                                                                                                                           >  1:2097012/1‑148 (MQ=255)
                           agGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGAtt                                                                                                                          >  2:3043266/1‑148 (MQ=255)
                             gCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTaa                                                                                                                        >  2:3255822/1‑148 (MQ=255)
                                gAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAgg                                                                                                                     <  1:1090654/148‑1 (MQ=255)
                                   aTGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGgcgc                                                                                                                  >  1:2395654/1‑148 (MQ=255)
                                     gCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTa                                                                                                                >  2:811258/1‑148 (MQ=255)
                                     gCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGAGTAAAGGCGCTa                                                                                                                >  1:2662668/1‑148 (MQ=255)
                                        cTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAg                                                                                                             <  1:941316/148‑1 (MQ=255)
                                         tGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGc                                                                                                            >  2:174191/1‑148 (MQ=255)
                                         tGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGc                                                                                                            >  1:2595196/1‑148 (MQ=255)
                                         tGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGc                                                                                                            >  1:2849808/1‑148 (MQ=255)
                                             tGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGga                                                                                                        >  2:1873712/1‑148 (MQ=255)
                                                 aTGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGgagaga                                                                                                    >  2:1103536/1‑148 (MQ=255)
                                                  tGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGaa                                                                                                   <  2:1391078/148‑1 (MQ=255)
                                                  tGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGaa                                                                                                   <  2:2213049/148‑1 (MQ=255)
                                                    cAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAgg                                                                                                 <  1:1320506/148‑1 (MQ=255)
                                                    cAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAgg                                                                                                 <  1:31589/148‑1 (MQ=255)
                                                    cAGTAGTGGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAgg                                                                                                 <  2:1726140/148‑1 (MQ=255)
                                                     agtagtTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGa                                                                                                <  1:903207/148‑1 (MQ=255)
                                                       tagtTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGAcc                                                                                              >  1:509671/1‑148 (MQ=255)
                                                        agtTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCg                                                                                             >  1:1588006/1‑148 (MQ=255)
                                                        agatGATTATCGTAAACAAAACTGACCTCATGTCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACNCTGGTATGCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCg                                                                                             >  1:1683010/4‑148 (MQ=255)
                                                         gtTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCgg                                                                                            >  2:838217/1‑148 (MQ=255)
                                                          ttGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGc                                                                                           >  1:45922/1‑148 (MQ=255)
                                                           tGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCa                                                                                          >  2:546907/1‑148 (MQ=255)
                                                             aTTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCaaa                                                                                        <  2:2596037/148‑1 (MQ=255)
                                                             aTTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTGGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCaaa                                                                                        <  2:2280608/148‑1 (MQ=255)
                                                                aTCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATg                                                                                     <  2:2850027/148‑1 (MQ=255)
                                                                 tCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGa                                                                                    >  2:41211/1‑148 (MQ=255)
                                                                 tCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGa                                                                                    <  1:2254357/148‑1 (MQ=255)
                                                                  cGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGAt                                                                                   >  1:1815927/1‑148 (MQ=255)
                                                                  cGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGAt                                                                                   >  2:1333288/1‑148 (MQ=255)
                                                                  cGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGAt                                                                                   >  1:3117820/1‑148 (MQ=255)
                                                                       aCAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTa                                                                              <  2:1929896/148‑1 (MQ=255)
                                                                          aaaCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACa                                                                           >  2:2629411/1‑148 (MQ=255)
                                                                           aaCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAg                                                                          >  1:2780139/1‑148 (MQ=255)
                                                                           aaCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAg                                                                          >  1:3077427/1‑148 (MQ=255)
                                                                           aaCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAg                                                                          >  1:1434864/1‑148 (MQ=255)
                                                                            aCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGt                                                                         <  2:246988/148‑1 (MQ=255)
                                                                            aCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGt                                                                         >  2:2347016/1‑148 (MQ=255)
                                                                            aCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGt                                                                         >  2:1471583/1‑148 (MQ=255)
                                                                             cTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGtt                                                                        <  1:197949/148‑1 (MQ=255)
                                                                             cTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGtt                                                                        >  1:635631/1‑148 (MQ=255)
                                                                              tGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGttt                                                                       <  1:1884208/148‑1 (MQ=255)
                                                                               gACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGtttt                                                                      <  2:1788130/148‑1 (MQ=255)
                                                                                  ctcctGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTc                                                                   <  2:1746712/148‑1 (MQ=255)
                                                                                  ctcctGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTc                                                                   >  2:2884205/1‑148 (MQ=255)
                                                                                    cctGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCAt                                                                 <  2:395989/148‑1 (MQ=255)
                                                                                      tGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAg                                                               <  1:799298/148‑1 (MQ=255)
                                                                                          gcgGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCat                                                           <  1:1004807/148‑1 (MQ=255)
                                                                                          gcgGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCat                                                           <  1:2932363/148‑1 (MQ=255)
                                                                                            ggATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCatat                                                         <  1:3043214/148‑1 (MQ=255)
                                                                                             gATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATc                                                        >  1:1430753/1‑148 (MQ=255)
                                                                                                tattCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCg                                                     >  1:572167/1‑148 (MQ=255)
                                                                                                 attCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGa                                                    >  1:253357/1‑148 (MQ=255)
                                                                                                 attCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGa                                                    <  2:2101118/148‑1 (MQ=255)
                                                                                                    cAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGa                                                 >  1:404416/1‑148 (MQ=255)
                                                                                                        ttGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACttt                                             >  2:3106221/1‑148 (MQ=255)
                                                                                                        ttGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACttt                                             >  1:323495/1‑148 (MQ=255)
                                                                                                        ttGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACTAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACttt                                             >  1:1287409/1‑148 (MQ=255)
                                                                                                         tGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTc                                            <  1:2191976/148‑1 (MQ=255)
                                                                                                             cGTCAACGGTTTTTCATCACACTGGTATTCCTGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTgg                                        >  2:2695650/1‑148 (MQ=255)
                                                                                                              gTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGt                                       >  1:2404937/1‑148 (MQ=255)
                                                                                                                  aCGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTgg                                   <  1:964932/148‑1 (MQ=255)
                                                                                                                   cGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGa                                  <  2:655852/148‑1 (MQ=255)
                                                                                                                   cGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGa                                  >  2:3082345/1‑148 (MQ=255)
                                                                                                                   cGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGa                                  <  1:137318/148‑1 (MQ=255)
                                                                                                                        tttCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATCATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTACTGGTCTGGAGGGCa                             >  1:1755331/1‑148 (MQ=255)
                                                                                                                         ttCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCaa                            <  2:1174686/148‑1 (MQ=255)
                                                                                                                          tcatcaCACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAAt                           <  2:1275557/148‑1 (MQ=255)
                                                                                                                          tcatcaCACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAAt                           <  1:3222062/148‑1 (MQ=255)
                                                                                                                          tcatcaCACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAAt                           <  1:153277/148‑1 (MQ=255)
                                                                                                                              cacaCTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGc                       <  1:2480587/148‑1 (MQ=255)
                                                                                                                                cacTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGccc                     >  2:2938647/1‑148 (MQ=255)
                                                                                                                                  cTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGccctc                   >  1:2996352/1‑148 (MQ=255)
                                                                                                                                   tGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGccctcc                  <  1:2751226/148‑1 (MQ=255)
                                                                                                                                   tGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGccctcc                  <  2:1124521/148‑1 (MQ=255)
                                                                                                                                   tGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGccctcc                  <  2:1305472/148‑1 (MQ=255)
                                                                                                                                      tATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTa               <  2:1649585/148‑1 (MQ=255)
                                                                                                                                      tATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTa               <  2:1588006/148‑1 (MQ=255)
                                                                                                                                         tCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACg            >  2:2035656/1‑148 (MQ=255)
                                                                                                                                         tCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACg            >  1:2856137/1‑148 (MQ=255)
                                                                                                                                          ccAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGt           <  2:1120414/148‑1 (MQ=255)
                                                                                                                                                 aTCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAgt    >  2:1560208/1‑148 (MQ=255)
                                                                                                                                                   cGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAgtgt  >  2:2667745/1‑148 (MQ=255)
                                                                                                                                                   cGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAgtgt  >  2:1264813/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGATTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAGTGT  >  CP009273/1712415‑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: