Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,712,562 C→A 27.8% 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,5620CA27.8% 476.7 / 153.6 230D31Y (GAT→TAT) mliCinhibitor of c‑type lysozyme, membrane‑bound; putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (95/71);  new base A (32/32);  total (127/103)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

Only 100 of 233 total aligned reads displayed.

TAGATTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAGTGT  >  CP009273/1712414‑1712709
                                                                                                                                                    |                                                                                                                                                   
tAGATTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATc                                                                                                                                                     <  1:869116/149‑1 (MQ=255)
    ttCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCgggt                                                                                                                                                 >  2:2396791/1‑147 (MQ=255)
    ttCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGtt                                                                                                                                                 <  1:165879/149‑1 (MQ=255)
     tCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGttt                                                                                                                                                <  2:2906730/149‑1 (MQ=255)
        gTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCa                                                                                                                                             <  1:44177/149‑1 (MQ=255)
        gTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCa                                                                                                                                             <  1:2443826/149‑1 (MQ=255)
         tCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCAt                                                                                                                                            >  2:348040/1‑149 (MQ=255)
         tCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCAt                                                                                                                                            <  1:2806512/149‑1 (MQ=255)
            gtgtAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATgcg                                                                                                                                         >  2:3325992/1‑149 (MQ=255)
             tgtAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATgcgc                                                                                                                                        <  1:828763/149‑1 (MQ=255)
             tgtAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATgcgc                                                                                                                                        >  1:3134510/1‑149 (MQ=255)
             tgtAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGAAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATgcgc                                                                                                                                        <  1:2890122/149‑1 (MQ=255)
              gtAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCt                                                                                                                                       >  2:2776592/1‑149 (MQ=255)
              gtAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCt                                                                                                                                       >  1:2501851/1‑149 (MQ=255)
                       gcCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCt                                                                                                                              >  2:955516/1‑149 (MQ=255)
                        ccAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTg                                                                                                                             >  2:2127249/1‑149 (MQ=255)
                           gagGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGAtt                                                                                                                          >  2:2886053/1‑149 (MQ=255)
                            agGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTa                                                                                                                         >  2:1843354/1‑149 (MQ=255)
                            agGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTa                                                                                                                         >  1:1229905/1‑149 (MQ=255)
                              gCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTaaa                                                                                                                       >  2:2433591/1‑149 (MQ=255)
                                  aaaTGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGgcg                                                                                                                   >  2:2214242/1‑149 (MQ=255)
                                    aTGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCt                                                                                                                 <  1:1182570/149‑1 (MQ=255)
                                      gCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTAc                                                                                                               >  2:583942/1‑149 (MQ=255)
                                       cccTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACa                                                                                                              <  2:3365297/149‑1 (MQ=255)
                                             ttGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGga                                                                                                        <  1:1949568/149‑1 (MQ=255)
                                                 gaTGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGgagaga                                                                                                    <  2:113937/149‑1 (MQ=255)
                                                 gaTGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGgagaga                                                                                                    <  2:2943419/149‑1 (MQ=255)
                                                      agtagtTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGAc                                                                                               <  1:3107489/149‑1 (MQ=255)
                                                       gtagtTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGAcc                                                                                              <  2:2543464/149‑1 (MQ=255)
                                                       gtagtTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGAcc                                                                                              <  1:1637689/149‑1 (MQ=255)
                                                               ttATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAAt                                                                                      >  1:461179/1‑149 (MQ=255)
                                                               ttATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAAt                                                                                      <  2:3495231/149‑1 (MQ=255)
                                                               ttATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAAt                                                                                      >  1:2534684/1‑149 (MQ=255)
                                                                 aTCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGa                                                                                    <  1:1086960/149‑1 (MQ=255)
                                                                 aTCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGa                                                                                    <  1:3841485/149‑1 (MQ=255)
                                                                 aTCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGa                                                                                    <  1:3288550/149‑1 (MQ=255)
                                                                  tCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGAt                                                                                   >  2:2391028/1‑149 (MQ=255)
                                                                    gTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATaa                                                                                 <  2:3374693/149‑1 (MQ=255)
                                                                     tAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAAt                                                                                >  1:2123417/1‑149 (MQ=255)
                                                                        aCAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTaa                                                                             <  2:529210/149‑1 (MQ=255)
                                                                         cAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAAc                                                                            <  2:1000742/149‑1 (MQ=255)
                                                                          aaaaCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACa                                                                           <  1:1103330/149‑1 (MQ=255)
                                                                           aaaCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAg                                                                          >  1:80809/1‑149 (MQ=255)
                                                                            aaCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGt                                                                         <  1:2213943/149‑1 (MQ=255)
                                                                             aCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGtt                                                                        <  2:2840780/149‑1 (MQ=255)
                                                                               tGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGtttt                                                                      >  1:3225987/1‑149 (MQ=255)
                                                                               tGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGtttt                                                                      >  1:2406438/1‑149 (MQ=255)
                                                                                gACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGttttt                                                                     <  1:3950624/149‑1 (MQ=255)
                                                                                 aCCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGtttttt                                                                    >  2:1869296/1‑149 (MQ=255)
                                                                                 aCCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGtttttt                                                                    >  2:3953509/1‑149 (MQ=255)
                                                                                 aCCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGtttttt                                                                    >  1:3159063/1‑149 (MQ=255)
                                                                                    tcctGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCAt                                                                 <  1:1508441/149‑1 (MQ=255)
                                                                                    tcctGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCAt                                                                 <  2:4031496/149‑1 (MQ=255)
                                                                                     cctGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATa                                                                <  1:1843354/149‑1 (MQ=255)
                                                                                      ctGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAg                                                               <  1:2309773/149‑1 (MQ=255)
                                                                                        ggCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTc                                                             >  2:3616752/1‑149 (MQ=255)
                                                                                        ggCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTc                                                             >  1:2679999/1‑149 (MQ=255)
                                                                                          cgcgGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCat                                                           <  1:4067712/149‑1 (MQ=255)
                                                                                          cgcgGATTATTCAGTTTGACCGTCAACGGTTTTCCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCat                                                           <  1:683630/149‑1 (MQ=255)
                                                                                           gcgGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCata                                                          >  1:1429618/1‑149 (MQ=255)
                                                                                            cgGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCatat                                                         <  1:502439/149‑1 (MQ=255)
                                                                                             ggATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATc                                                        >  1:3443031/1‑149 (MQ=255)
                                                                                             ggATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATc                                                        >  1:2102587/1‑149 (MQ=255)
                                                                                             ggATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATc                                                        >  1:1346589/1‑149 (MQ=255)
                                                                                              gATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATcc                                                       <  2:1886519/149‑1 (MQ=255)
                                                                                                ttattCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCg                                                     >  2:861227/1‑149 (MQ=255)
                                                                                                ttattCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCg                                                     <  1:291363/149‑1 (MQ=255)
                                                                                                ttattCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCg                                                     >  2:1097570/1‑149 (MQ=255)
                                                                                                 tattCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGa                                                    >  2:212497/1‑149 (MQ=255)
                                                                                                 tattCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGa                                                    >  2:3904109/1‑149 (MQ=255)
                                                                                                  attCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGaa                                                   >  1:3181631/1‑149 (MQ=255)
                                                                                                  attCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTCCAGCCGGAGAGAAGGACCGGCACAATGATCATTCACAGTTTTGTCGTAGTCATATCCCGac                                                   >  1:2917687/1‑148 (MQ=255)
                                                                                                  attCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGaa                                                   <  1:3821857/149‑1 (MQ=255)
                                                                                                       gTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACtt                                              >  1:1463626/1‑149 (MQ=255)
                                                                                                       gTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACtt                                              >  2:1262606/1‑149 (MQ=255)
                                                                                                        tttGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACttt                                             >  2:3882422/1‑149 (MQ=255)
                                                                                                          tGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTcc                                           <  2:1443889/149‑1 (MQ=255)
                                                                                                          tGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTcc                                           <  2:2294944/149‑1 (MQ=255)
                                                                                                           gACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCt                                          <  2:714064/149‑1 (MQ=255)
                                                                                                                tCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCt                                     <  1:3627847/149‑1 (MQ=255)
                                                                                                                tCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCt                                     <  1:1342258/149‑1 (MQ=255)
                                                                                                                tCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCt                                     <  2:2102587/149‑1 (MQ=255)
                                                                                                                   aCGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCCTATCCCGAAGACTTTCCTGGTCTGGa                                  >  1:1324434/1‑149 (MQ=255)
                                                                                                                       tttttCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGc                              >  2:2025669/1‑149 (MQ=255)
                                                                                                                         tttCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCaa                            <  2:807941/149‑1 (MQ=255)
                                                                                                                           tcatcaCACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATa                          >  1:1702609/1‑149 (MQ=255)
                                                                                                                            catcaCACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATAc                         >  2:2547701/1‑149 (MQ=255)
                                                                                                                             atcaCACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATGGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACg                        >  2:3785935/1‑149 (MQ=255)
                                                                                                                             atcaCACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACg                        >  1:609346/1‑149 (MQ=255)
                                                                                                                             atcaCACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACg                        >  2:1700107/1‑149 (MQ=255)
                                                                                                                                   cTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGccctcc                  >  1:3696086/1‑149 (MQ=255)
                                                                                                                                    tGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGccctccc                 <  2:1771579/149‑1 (MQ=255)
                                                                                                                                       tATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTaa              >  2:2245293/1‑149 (MQ=255)
                                                                                                                                         ttCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACg            >  2:1235535/1‑149 (MQ=255)
                                                                                                                                                 tATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAgt    <  2:2501851/149‑1 (MQ=255)
                                                                                                                                                 tATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAgt    >  1:1668871/1‑149 (MQ=255)
                                                                                                                                                 tATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAgt    >  2:3767609/1‑149 (MQ=255)
                                                                                                                                                  aTCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAgtg   <  1:2776592/149‑1 (MQ=255)
                                                                                                                                                   tCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAgtgt  <  1:227801/149‑1 (MQ=255)
                                                                                                                                                   tCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAgtgt  >  2:310091/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TAGATTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAGTGT  >  CP009273/1712414‑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: