Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,931,898 G→T 6.1% M104I (ATG→ATT rbsR → transcriptional repressor of ribose metabolism

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,931,8980GT 6.1% 1187.6 / 28.9 318M104I (ATG→ATTrbsRtranscriptional repressor of ribose metabolism
Reads supporting (aligned to +/- strand):  ref base G (147/151);  new base T (12/7);  total (160/158)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 322 total aligned reads displayed.

AACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCGACAG  >  CP009273/3931750‑3932047
                                                                                                                                                    |                                                                                                                                                     
aaCTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGa                                                                                                                                                      >  1:3411172/1‑150 (MQ=255)
  cTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAAt                                                                                                                                                    <  2:325825/150‑1 (MQ=255)
    caatcaaACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCg                                                                                                                                                  >  2:3592445/1‑150 (MQ=255)
      atcaaACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCa                                                                                                                                                >  1:2447938/1‑150 (MQ=255)
       tcaaACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCaa                                                                                                                                               <  1:3098473/150‑1 (MQ=255)
        caaACACATACCAATGGCATGTTGATCACTGCCAGTACCAATCCTTTCTAGTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCACGTCCTATGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAAt                                                                                                                                              <  1:2374693/150‑1 (MQ=255)
         aaaCACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATTAATCGCAATc                                                                                                                                             <  1:5272446/150‑1 (MQ=255)
         aaaCACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATc                                                                                                                                             <  2:3065895/150‑1 (MQ=255)
          aaCACATACCATTGGCATGTTGATCTCTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCt                                                                                                                                            <  1:2966701/150‑1 (MQ=255)
          aaCACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTAGTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCt                                                                                                                                            <  2:1509455/150‑1 (MQ=255)
                 aCCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACg                                                                                                                                     <  2:4062984/150‑1 (MQ=255)
                    aTTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTg                                                                                                                                  <  1:3879478/150‑1 (MQ=255)
                    aTTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTg                                                                                                                                  >  1:4029550/1‑150 (MQ=255)
                       ggCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATg                                                                                                                               <  2:2781602/150‑1 (MQ=255)
                       ggCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGAGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCATTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATg                                                                                                                               >  2:2727930/1‑150 (MQ=255)
                         cATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAATGCTGATGCa                                                                                                                             >  1:3324226/1‑150 (MQ=255)
                         cATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCa                                                                                                                             >  1:817136/1‑150 (MQ=255)
                         cATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGGTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATTAATCGCAATCTGGAAACGCTGATGCa                                                                                                                             >  2:2688842/1‑150 (MQ=255)
                           tGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGATACGCTGATGCaaa                                                                                                                           >  2:3863586/1‑150 (MQ=255)
                            gTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCaaaa                                                                                                                          >  2:2062000/1‑150 (MQ=255)
                            gTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCaaaa                                                                                                                          <  2:2580192/150‑1 (MQ=255)
                              tGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATTAATCGCAATCTGGAAACGCTGATGCAAAAAc                                                                                                                        >  1:2915195/1‑150 (MQ=255)
                              tGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATTAATCGCAATCTGGAAACGCTGATGCAAAAAc                                                                                                                        >  2:4715319/1‑150 (MQ=255)
                              tGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAAc                                                                                                                        <  2:1716083/150‑1 (MQ=255)
                              tGATCACTGCCAGTACCAATCCATTCTATTCAGAACTGGAGCGTTGCGTAGAACGCAGCTGCTTCGAACGCGCTTATAGTATCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATAAATCGCAATCTGGAAACGCTCATGCTAAAAc                                                                                                                        >  2:3709957/1‑150 (MQ=255)
                                   aCTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGtt                                                                                                                   <  1:5133668/150‑1 (MQ=255)
                                   aCTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGtt                                                                                                                   <  2:1545811/150‑1 (MQ=255)
                                    cTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTg                                                                                                                  <  1:154949/150‑1 (MQ=255)
                                     tGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGa                                                                                                                 <  1:4723990/150‑1 (MQ=255)
                                      gCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTAGAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGAt                                                                                                                <  2:677266/150‑1 (MQ=255)
                                         aGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGGAACGCTGATGCAAAAACGCGTTGATGGc                                                                                                             >  2:440457/1‑150 (MQ=255)
                                          gTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCt                                                                                                            >  2:1864813/1‑150 (MQ=255)
                                            aCCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTg                                                                                                          >  1:2087267/1‑150 (MQ=255)
                                            aCCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTg                                                                                                          >  2:5314779/1‑150 (MQ=255)
                                             ccAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGCTGACCAGCGGATGAATCGCAATCTGGACACGCTGATGCAAAAACGCGTTCATGGCTTGc                                                                                                         >  2:1388562/1‑150 (MQ=255)
                                             ccAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGc                                                                                                         >  2:341035/1‑150 (MQ=255)
                                              cAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCt                                                                                                        <  1:366027/150‑1 (MQ=255)
                                               aaTCCTTTCTATTCAGTACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTg                                                                                                       <  1:2570032/150‑1 (MQ=255)
                                               aaTCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTg                                                                                                       >  1:1176496/1‑150 (MQ=255)
                                                    tttCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACt                                                                                                  >  2:229151/1‑150 (MQ=255)
                                                        tATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGc                                                                                              >  2:3805412/1‑150 (MQ=255)
                                                         aTTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCa                                                                                             >  2:4602173/1‑150 (MQ=255)
                                                         aTTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTAATGTGCa                                                                                             >  1:4519005/1‑150 (MQ=255)
                                                            cAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCg                                                                                          >  2:5102334/1‑150 (MQ=255)
                                                             aGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGa                                                                                         >  2:4274492/1‑150 (MQ=255)
                                                             aGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGa                                                                                         >  2:831485/1‑150 (MQ=255)
                                                             aGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGa                                                                                         >  1:1938356/1‑150 (MQ=255)
                                                              gAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGaa                                                                                        >  2:3630142/1‑150 (MQ=255)
                                                               aaCTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGaaa                                                                                       >  2:3593043/1‑150 (MQ=255)
                                                                aCTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAAc                                                                                      <  1:948270/150‑1 (MQ=255)
                                                                aCTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAAc                                                                                      <  1:2529423/150‑1 (MQ=255)
                                                                 cTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACg                                                                                     <  1:4420096/150‑1 (MQ=255)
                                                                 cTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAAGACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACg                                                                                     <  2:4758802/150‑1 (MQ=255)
                                                                  tGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGc                                                                                    <  2:3305843/150‑1 (MQ=255)
                                                                   ggTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCa                                                                                   >  2:2762297/1‑150 (MQ=255)
                                                                          ggCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCt                                                                            >  1:1134511/1‑150 (MQ=255)
                                                                           gCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCtt                                                                           <  1:4744016/150‑1 (MQ=255)
                                                                             gTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTcg                                                                         >  1:2785929/1‑150 (MQ=255)
                                                                             gTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTcg                                                                         >  1:3253237/1‑150 (MQ=255)
                                                                             gTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTAGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTcg                                                                         <  1:1877586/150‑1 (MQ=255)
                                                                              ttGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTcgc                                                                        <  2:4965181/150‑1 (MQ=255)
                                                                                 aaCGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTg                                                                     <  1:4597501/150‑1 (MQ=255)
                                                                                     cAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAAt                                                                 <  1:2208232/150‑1 (MQ=255)
                                                                                     cAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAAt                                                                 <  2:2365415/150‑1 (MQ=255)
                                                                                     cAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAAt                                                                 <  1:3769768/150‑1 (MQ=255)
                                                                                       gctgctTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCa                                                               <  2:1166058/150‑1 (MQ=255)
                                                                                         tgctTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATg                                                             >  2:1896891/1‑150 (MQ=255)
                                                                                          gctTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGc                                                            <  1:218076/150‑1 (MQ=255)
                                                                                          gctTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGc                                                            >  1:1071530/1‑150 (MQ=255)
                                                                                              cGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACg                                                        >  1:4269602/1‑150 (MQ=255)
                                                                                                aaCGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGtt                                                      >  2:4186721/1‑150 (MQ=255)
                                                                                                 aCGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTa                                                     >  1:2326917/1‑150 (MQ=255)
                                                                                                   gcgGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATc                                                   <  1:1193574/150‑1 (MQ=255)
                                                                                                    cgtttATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAATCGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATcc                                                  <  2:283474/147‑1 (MQ=255)
                                                                                                    cgGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATcc                                                  <  2:5271951/150‑1 (MQ=255)
                                                                                                                cGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATTAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTAc                                      >  2:4659142/1‑150 (MQ=255)
                                                                                                                cGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTAc                                      >  1:1665467/1‑150 (MQ=255)
                                                                                                                cGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTAc                                      >  2:2168890/1‑150 (MQ=255)
                                                                                                                 gTCCTTTGCAAGACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACt                                     <  1:3192783/150‑1 (MQ=255)
                                                                                                                      ttGCAATACCGAAGGCGATGAACAGCGGATGATTCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGtgat                                <  2:3735530/150‑1 (MQ=255)
                                                                                                                        gCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGtgatga                              <  2:2447938/150‑1 (MQ=255)
                                                                                                                          aaTACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGattatg                            >  1:3716818/1‑150 (MQ=255)
                                                                                                                            tACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGa                          <  2:493385/150‑1 (MQ=255)
                                                                                                                             aCCGAAGGCGATGAACAGCGGATTAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGAc                         <  2:636663/150‑1 (MQ=255)
                                                                                                                             aCCGAAGGCGATGAACAGCGGATTAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGAc                         <  1:4715359/150‑1 (MQ=255)
                                                                                                                                   ggCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCt                   <  1:3978385/150‑1 (MQ=255)
                                                                                                                                       aTGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGt               >  2:1081236/1‑150 (MQ=255)
                                                                                                                                        tGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGtt              <  2:5028816/150‑1 (MQ=255)
                                                                                                                                         gAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTc             <  1:30221/150‑1 (MQ=255)
                                                                                                                                         gAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTAATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTc             >  2:3346524/1‑150 (MQ=255)
                                                                                                                                             aGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGAATGGGCTCCGTTCGATg         <  2:3184927/150‑1 (MQ=255)
                                                                                                                                              gCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATgg        <  2:1629351/150‑1 (MQ=255)
                                                                                                                                               cGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGc       <  1:13150/150‑1 (MQ=255)
                                                                                                                                               cGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGc       >  2:3999581/1‑150 (MQ=255)
                                                                                                                                                  aTGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCGAc    <  1:3163962/150‑1 (MQ=255)
                                                                                                                                                  aTGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGAAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCGAc    <  1:568653/150‑1 (MQ=255)
                                                                                                                                                    gAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACTGTATCCGACAGTGCCTACTGTGAAGATGGACTGGGCTCCGTTCGATGTCGACAg  >  1:174085/1‑150 (MQ=255)
                                                                                                                                                    gAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCGACAg  >  1:4810038/1‑150 (MQ=255)
                                                                                                                                                    gAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCGACAg  <  1:4685652/150‑1 (MQ=255)
                                                                                                                                                    gAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCGACAg  >  2:1954446/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCGACAG  >  CP009273/3931750‑3932047

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: