Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,987,103 C→T 5.1% intergenic (+27/‑433) yqeG → / → yqeH putative transporter/putative LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,987,1030CT 5.1% 1068.0 / 28.7 316intergenic (+27/‑433)yqeG/yqeHputative transporter/putative LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (159/139);  new base T (12/4);  total (173/143)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 316 total aligned reads displayed.

TTTTCATCATGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCAT  >  NC_000913/2986955‑2987247
                                                                                                                                                    |                                                                                                                                                
      tcatGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTACTGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAAtt                                                                                                                                           >  1:1211932/1‑150 (MQ=255)
      tcatGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAAtt                                                                                                                                           >  1:4866556/1‑150 (MQ=255)
      tcatGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTGAGCTCATATCAGTTGTGGATGGTTCAAAATTGCTACAAtt                                                                                                                                           >  2:5685484/1‑150 (MQ=255)
       catGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAAttt                                                                                                                                          >  1:3696494/1‑150 (MQ=255)
        atGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAAtttt                                                                                                                                         <  1:4327930/150‑1 (MQ=255)
        atGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAAtttt                                                                                                                                         >  2:2432038/1‑150 (MQ=255)
        atGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAAtttt                                                                                                                                         >  1:2862430/1‑150 (MQ=255)
         tGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTg                                                                                                                                        <  2:2692599/150‑1 (MQ=255)
          gCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGTAGCAATTTTGa                                                                                                                                       <  1:2411430/150‑1 (MQ=255)
           ccTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGaa                                                                                                                                      <  1:5501018/150‑1 (MQ=255)
           ccTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGaa                                                                                                                                      <  1:5333812/150‑1 (MQ=255)
           ccTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCGGATCGTGGGTATCCTGTGCGTCTCGGGAATGTTCGTTAGCGAAGAGCCGGTGTGGATGGGTCAAAATTGCTACAATTTTGaa                                                                                                                                      >  1:2198270/1‑150 (MQ=255)
            cTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTTAAc                                                                                                                                     >  1:1154155/1‑150 (MQ=255)
            cTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAAc                                                                                                                                     >  2:2622244/1‑150 (MQ=255)
                   gTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAAccattcca                                                                                                                              <  2:2388736/150‑1 (MQ=255)
                      aaCGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCAttt                                                                                                                           >  2:5349747/1‑150 (MQ=255)
                       aCGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCAtttt                                                                                                                          >  1:591437/1‑150 (MQ=255)
                       aCGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCAGTGCGTATCGGGAATGTTATTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTTAAACATTCCAtttt                                                                                                                          >  1:1989882/1‑150 (MQ=255)
                       aCGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCAGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCAtttt                                                                                                                          >  2:5210677/1‑150 (MQ=255)
                              tcATCCCGGCGCTTAAACCCTGGCGTTCCAACGGAAATATGATTACGCTGATCGTGGGTATACTGTGCGTATCGGTAATGTTCTTTAGCTAAAATCAGTTGTGGATGGTTCAAAATTGCTACAAATTTTAACCATTACAATTTTCctaat                                                                                                                   >  1:3048540/1‑150 (MQ=255)
                               cATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCataata                                                                                                                  >  2:1244646/1‑150 (MQ=255)
                                     ggCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTAAGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCAtt                                                                                                            >  1:935645/1‑150 (MQ=255)
                                           tAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGttat                                                                                                      >  2:2737091/1‑150 (MQ=255)
                                             aaCCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGttatta                                                                                                    >  1:994813/1‑150 (MQ=255)
                                                ccTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTc                                                                                                 <  2:4979000/150‑1 (MQ=255)
                                                ccTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTc                                                                                                 >  1:3980147/1‑150 (MQ=255)
                                                 cTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTcc                                                                                                <  1:785866/150‑1 (MQ=255)
                                                 cTGGCGTTACATCGGAAATCTGAGTAGGCTGATCGTGGGTATCCTGTGAGTATCGGTTATGTTCTTTAGCTCATTTCTGTTGTGGTTGGTTCATTATTGCTTCAATTTTGTAGCATTCCTTTTTTCATTCTATTCATTTGTTattattgc                                                                                                <  2:4301945/150‑3 (MQ=255)
                                                    gCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCAcc                                                                                             <  1:1832041/150‑1 (MQ=255)
                                                    gCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCAcc                                                                                             >  2:4676260/1‑150 (MQ=255)
                                                    gCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCAcc                                                                                             <  1:4569051/150‑1 (MQ=255)
                                                     cGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCGGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCAccc                                                                                            <  2:3775763/150‑1 (MQ=255)
                                                          cATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTACACCCTATCa                                                                                       >  1:830430/1‑150 (MQ=255)
                                                           aTCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCAc                                                                                      >  1:3115772/1‑150 (MQ=255)
                                                            tCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACt                                                                                     >  1:3302741/1‑150 (MQ=255)
                                                            tCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACt                                                                                     <  1:229473/150‑1 (MQ=255)
                                                             cGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACtt                                                                                    <  1:51535/150‑1 (MQ=255)
                                                               gAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtct                                                                                  <  2:78839/150‑1 (MQ=255)
                                                                aaaTCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtctc                                                                                 <  1:208654/150‑1 (MQ=255)
                                                                 aaTCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtctct                                                                                <  1:3483525/150‑1 (MQ=255)
                                                                 aaTCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtctct                                                                                <  2:2914700/150‑1 (MQ=255)
                                                                 aaTCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtctct                                                                                <  2:605290/150‑1 (MQ=255)
                                                                  aTCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtctctc                                                                               <  1:1734463/150‑1 (MQ=255)
                                                                   tCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCg                                                                              <  2:1443698/150‑1 (MQ=255)
                                                                    cTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGt                                                                             <  1:1104393/150‑1 (MQ=255)
                                                                      gATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATTTTCTTTAGCTAATATCATTTGTGGATGGTTCATAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCa                                                                           <  1:731995/150‑1 (MQ=255)
                                                                      gATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCa                                                                           >  2:3598999/1‑150 (MQ=255)
                                                                        ttACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCTTTCCATTTTTCATAGTTTTGATTTGTTATTATTCCACCGTTTCAGTTCTCTCGtggtt                                                                         >  2:5430210/1‑146 (MQ=255)
                                                                         tACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCatat                                                                        >  2:3867920/1‑150 (MQ=255)
                                                                         tACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCatat                                                                        <  1:3230692/150‑1 (MQ=255)
                                                                         tACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCatat                                                                        <  2:5152980/150‑1 (MQ=255)
                                                                         tACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTATTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCatat                                                                        <  1:1044625/150‑1 (MQ=255)
                                                                          aCGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATc                                                                       <  1:1245175/150‑1 (MQ=255)
                                                                            gCTGATCGTGGGTATCCTGTGCGTTTAGGTAATGTTCTTTAGGTAATATCTGTTGTGGATGGTTCAATATTGCTACATTTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCt                                                                     <  1:2046205/150‑1 (MQ=255)
                                                                             cTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTa                                                                    <  2:3142540/150‑1 (MQ=255)
                                                                             cTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAACTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTa                                                                    >  2:1762012/1‑150 (MQ=255)
                                                                             atGATAGTGGGTATCCTGTGAGTATAGGTAATGTTCTTTAGCTAATATCAGTTTTGGATGGTTCAAAATTGATACAATTTTGAACCATTCCATTTATCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTa                                                                    >  2:165786/2‑150 (MQ=255)
                                                                                 tCGTGGGTATCCTGTGCGTGTCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGtat                                                                >  2:2281335/1‑150 (MQ=255)
                                                                                      ggTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCa                                                           >  1:3810751/1‑150 (MQ=255)
                                                                                      ggTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCa                                                           <  2:5890816/150‑1 (MQ=255)
                                                                                        tATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagt                                                         <  2:5561828/150‑1 (MQ=255)
                                                                                          tCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagtta                                                       <  1:4976723/150‑1 (MQ=255)
                                                                                            cTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagttagt                                                     >  1:2573291/1‑150 (MQ=255)
                                                                                            cTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagttagt                                                     >  2:2218006/1‑150 (MQ=255)
                                                                                            cTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCGAGTCTTATCagttagt                                                     >  2:2082720/1‑150 (MQ=255)
                                                                                              gtgCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATACTAGTATTATCAgttagtta                                                   >  1:628977/1‑150 (MQ=255)
                                                                                               tgCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTaa                                                  <  2:4170162/150‑1 (MQ=255)
                                                                                                  gTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAaca                                               <  1:4428223/150‑1 (MQ=255)
                                                                                                     tCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAacacac                                            <  1:2607407/150‑1 (MQ=255)
                                                                                                     tCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAacacac                                            >  2:3485498/1‑150 (MQ=255)
                                                                                                       ggTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACaa                                          >  1:3728249/1‑150 (MQ=255)
                                                                                                       ggTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACaa                                          >  1:4107458/1‑150 (MQ=255)
                                                                                                        gTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACaaa                                         >  2:5772349/1‑150 (MQ=255)
                                                                                                        gTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACaaa                                         <  2:104994/150‑1 (MQ=255)
                                                                                                        gTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAAAACACaaa                                         >  1:4145324/1‑150 (MQ=255)
                                                                                                          aaTGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAaga                                       <  2:4866556/150‑1 (MQ=255)
                                                                                                            ttttCTTTAGCTAATATCAGTTGTGTATGGTTCAAAATTGCTAAAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTTGTATTATCAGTTAGTTAAACACACAAAtaga                                     <  1:521870/148‑1 (MQ=255)
                                                                                                             ggtCTTTAGCTACTATCAGTTGTGGATGGTTCCAAAGTGCTACAATTTTGAACCAGTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTATCCACACAAAGAGAt                                    <  2:5078290/148‑1 (MQ=255)
                                                                                                             gTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGAt                                    >  2:4557946/1‑150 (MQ=255)
                                                                                                               tCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCt                                  <  2:3732944/150‑1 (MQ=255)
                                                                                                                      cTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATc                           >  2:5013685/1‑150 (MQ=255)
                                                                                                                       tAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCt                          <  1:5079514/150‑1 (MQ=255)
                                                                                                                         atatCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTa                        <  1:4628229/150‑1 (MQ=255)
                                                                                                                              aGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATaa                   <  1:1762410/150‑1 (MQ=255)
                                                                                                                               aaatggaaTGGTTCAAAATTGTAGCAATTTTGAACCATCCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAg                  >  2:91996/8‑150 (MQ=38)
                                                                                                                                     gATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTAtt            <  2:115119/150‑1 (MQ=255)
                                                                                                                                     gATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTAtt            >  1:2506029/1‑150 (MQ=255)
                                                                                                                                        ggTTCAAAATTGCTACACTTTTGAACCATTCCAATTTTCATAATATTCATTTGTTATTATTCCACCCTCTCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTaa         <  2:1957441/150‑1 (MQ=255)
                                                                                                                                        ggTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTaa         <  2:3014489/150‑1 (MQ=255)
                                                                                                                                          ttCAAAATTGTAGCAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagt                                                         >  1:3606150/1‑100 (MQ=255)
                                                                                                                                          ttCAAAATTGTAGCAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGa       >  1:3910362/1‑150 (MQ=255)
                                                                                                                                          ttCAAAATTGTAGCAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGa       >  1:3139214/1‑150 (MQ=255)
                                                                                                                                          ttCAAAATTGTAGCAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGGTAGTTAAACACACAAAGAGATCTGTACATCTTAAATGAGACTAttcgcgg       >  2:3664869/1‑145 (MQ=255)
                                                                                                                                          ttCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGa       >  2:179715/1‑150 (MQ=255)
                                                                                                                                          ttCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGa       >  2:1690558/1‑150 (MQ=255)
                                                                                                                                          ttCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGa       >  1:878017/1‑150 (MQ=255)
                                                                                                                                          ttCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGa       >  1:2748775/1‑150 (MQ=255)
                                                                                                                                           tCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGAt      >  2:1273943/1‑150 (MQ=255)
                                                                                                                                           tCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGAt      >  1:1658346/1‑150 (MQ=255)
                                                                                                                                              aaaTTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCa   <  2:1525423/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TTTTCATCATGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCAT  >  NC_000913/2986955‑2987247

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: