Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,504,880 A→C 7.5% intergenic (+97/+55) frlR → / ← yhfS putative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,504,8800AC 7.5% 505.6 / 7.1 183intergenic (+97/+55)frlR/yhfSputative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base A (67/98);  new base C (12/6);  total (79/104)
Fisher's exact test for biased strand distribution p-value = 4.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 186 total aligned reads displayed.

TTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTG  >  NC_000913/3504733‑3505029
                                                                                                                                                   |                                                                                                                                                     
ttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCaaa                                                                                                                                                     <  1:2423428/150‑1 (MQ=255)
 tctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATATGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAg                                                                                                                                                    <  2:2543755/150‑1 (MQ=255)
 tctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAg                                                                                                                                                    >  2:4001007/1‑150 (MQ=255)
    gCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAgggg                                                                                                                                                 <  2:3402973/150‑1 (MQ=255)
     ccGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCACACCCACACCCACACCCCCAAGGGGc                                                                                                                                                >  1:3611294/1‑150 (MQ=38)
      cgcgCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCg                                                                                                                                               <  2:3176817/150‑1 (MQ=255)
       gcgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACAACCAAACCCACACCCCAAAAGGGCGa                                                                                                                                              >  2:810622/1‑150 (MQ=39)
          ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAggg                                                                                                                                           <  1:24081/150‑1 (MQ=255)
            actCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGa                                                                                                                                         <  1:1988930/148‑1 (MQ=255)
              tCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGAcc                                                                                                                                       <  1:1871147/150‑1 (MQ=255)
                  aTCACCTTAACTATTGATAATAAAAGACATTAACCGTCGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCGAACCCCCACCCCCAAGGGGCGGGGGGGCCgggc                                                                                                                                   >  1:2032811/1‑147 (MQ=38)
                   tCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCa                                                                                                                                  <  1:3835760/150‑1 (MQ=255)
                   tCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCa                                                                                                                                  <  1:146010/150‑1 (MQ=255)
                      ccTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAAACCACTCCCCAAAGGGGCGAGGGGGCCGGGCAGTa                                                                                                                               >  1:1749167/1‑150 (MQ=38)
                               ttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCGCACCCAAACCCACAACCCAAAAGGGCGAGGGGGCCCGTCAGTACGGAAATTa                                                                                                                      >  1:2889799/1‑150 (MQ=37)
                                cgATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAg                                                                                                                     <  1:772105/149‑1 (MQ=255)
                                 gATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCACACCCAAAACAACACCCCCAACGGGCGCGGGGACCGGTCAGTAAGGAATTTatg                                                                                                                    >  2:1569729/1‑148 (MQ=37)
                                 gATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGc                                                                                                                    <  2:189555/150‑1 (MQ=255)
                                 gATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCTAACCCACACCCCAAAGGGGCGTGGGGACCCGTCAGTACCGAAATTTGc                                                                                                                    >  2:1692926/1‑150 (MQ=38)
                                  ataataAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCTCACCCAAAACCACAACCCAAAGGGGGGAGGGGGCCGGTAAGTTCGGAAATTAgga                                                                                                                   >  1:4021761/1‑148 (MQ=37)
                                  ataataAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCt                                                                                                                   <  2:2204797/150‑1 (MQ=255)
                                   taataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCACACCCGAACCCACACCCCAAAGGGGCGCGGGGACCCGTCAGTTCGGAAGTTAGcat                                                                                                                  >  2:3599689/1‑148 (MQ=37)
                                      taaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCTAACCCTCTACCCAAAGGGGCGAGGGGACCGGTCAGTACGGAATTTAGGATGAc                                                                                                               >  2:4088829/1‑150 (MQ=39)
                                       aaaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCACACCCGCTCCCCCAAGGGGCGAGGGGCCCCTTCACTACGGAAGTTAGCTTGACg                                                                                                              >  2:2876163/1‑150 (MQ=38)
                                          aGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTg                                                                                                           <  1:1076270/150‑1 (MQ=255)
                                          aGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTg                                                                                                           <  1:1829268/150‑1 (MQ=255)
                                           gACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGc                                                                                                          <  1:3727100/150‑1 (MQ=255)
                                               ttAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATc                                                                                                      <  1:3675790/150‑1 (MQ=255)
                                               ttAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATc                                                                                                      <  1:674848/150‑1 (MQ=255)
                                                tAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGGTCAGTACGGAATTTAGCTTGACGGTGCTATcc                                                                                                     <  2:1972212/150‑1 (MQ=255)
                                                 taCCGTAGGCAGGATAAGATGCGCCAGCATCGCATCAGGCTATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCTCACACTAACCCTCTCCACAAAGGGGCGAGGGGTCCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCa                                                                                                    <  2:758704/149‑1 (MQ=255)
                                                 aaCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCa                                                                                                    <  1:3965835/150‑1 (MQ=255)
                                                 aaCCGTAGGCCGGATAAGATGCGCCAGCAACGCATCCGGCGATGCTGGCGCGATGAATATTACATCCCGTACCTTCCCCTCACCCTAACCCTCTCCCCAAACGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCa                                                                                                    <  2:2510098/150‑1 (MQ=39)
                                                    cGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCa                                                                                                 <  2:243626/150‑1 (MQ=255)
                                                    cGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCa                                                                                                 <  1:2690369/150‑1 (MQ=255)
                                                    cGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACACCCAAAACCACTACCCAAAGGGGCGCGGGGGCCGTTCAGTACGGAATTTAGGTTGACCGTGCTATACAGCa                                                                                                 >  2:124407/1‑150 (MQ=37)
                                                      tAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCAtt                                                                                               <  2:2925995/150‑1 (MQ=255)
                                                        ggCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTaa                                                                                             <  1:3367441/150‑1 (MQ=255)
                                                              aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                                                       <  2:2773028/150‑1 (MQ=255)
                                                              aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                                                       <  2:2772608/150‑1 (MQ=255)
                                                               tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                                      <  1:486625/150‑1 (MQ=255)
                                                                    tGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATa                                                                                 <  2:2256245/150‑1 (MQ=255)
                                                                      cgcCAGCATCGCCGGATGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAACCCACACCCCAAAGGGGCGAGGGGACCCGTCAGTACGGAAGTTAGCTTGACGGTGCGATCCAGCAGTAAATAATGGCaaaaat                                                                               >  2:424608/1‑146 (MQ=25)
                                                                       gcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCGCCCCCGAACCCGCTCCCCCCCGGGGGGGGGGGGCCGGGCCGTACGGAAATTAGCTTTGCGGTGGTTGCCCGCCTTAAAAAGTGGGCATActc                                                                              >  2:1205037/1‑150 (MQ=255)
                                                                         cAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACtctc                                                                            <  2:2321790/150‑1 (MQ=255)
                                                                           gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                          <  1:3511144/150‑1 (MQ=255)
                                                                           gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                          <  2:3694230/150‑1 (MQ=255)
                                                                           gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                          <  2:3694228/150‑1 (MQ=255)
                                                                           gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                          <  2:1151572/150‑1 (MQ=255)
                                                                              tcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTATCTTGAAGGTGCTATCCAGCATTAAATAATGGCAATACACACAAGCa                                                                       >  2:341682/1‑150 (MQ=38)
                                                                                     cGGCGATTCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCg                                                                <  2:1291862/150‑1 (MQ=255)
                                                                                       gCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgca                                                              <  1:301341/150‑1 (MQ=255)
                                                                                            gCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCAcc                                                         <  1:921416/150‑1 (MQ=255)
                                                                                              tGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGt                                                       <  1:1589252/150‑1 (MQ=255)
                                                                                                gcgcGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCACTACCCACAGGGGCGAGGGGACCGGTCAGTACGGAATTTATCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTAACCCAAAATACGCAGCACCGtct                                                     >  2:506003/1‑150 (MQ=255)
                                                                                                gcgcGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAAACCACTCCCCAAAGGGGCCACGGGACCGGTCAGTACCGAATTTAGCATGGCGGTGCTATCCACCAATAAAAACGGGCAAAACACACACCCAAAAAACGCAGCACCctct                                                     >  2:2155942/1‑150 (MQ=255)
                                                                                                   cgTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCtt                                                  <  2:2133828/150‑1 (MQ=255)
                                                                                                         aTTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGc                                            <  1:4072241/150‑1 (MQ=255)
                                                                                                          ttttACATCCCGTACGTTCCCCACACCCAAACCCACTCCCCCAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCATCAATAAAAAAAGGCAAAACTCTCACGCACAAACCCCACCACCCACTCATCACCGCt                                           >  2:2209029/1‑150 (MQ=255)
                                                                                                            ttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAACATTAAAAACAGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCAGc                                         >  2:3466791/1‑150 (MQ=255)
                                                                                                              aCATCCCGTACGTTCCCCTCACCCTAACCCTCGCCCCCCAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGa                                       <  2:2717730/150‑1 (MQ=255)
                                                                                                                   ccGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCgg                                  <  2:3449294/150‑1 (MQ=255)
                                                                                                                    cGTACGTTCCCCTCACCCTAACCCTCTCCCCAAATGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCggg                                 >  2:238314/1‑150 (MQ=255)
                                                                                                                    cGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCggg                                 >  2:3079054/1‑150 (MQ=255)
                                                                                                                    cGTACGTTCCCCTCACCCTAAACCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAAATTAGCTTGAAGGTGCTATCCACCATTAAATAATGGCAAAACTCTCACGCAAAATCCGCACCACCCTCTCTTCACCGCTGCGATTCggg                                 >  2:399064/1‑150 (MQ=38)
                                                                                                                    cGTACGTTCCCCTCACCCTAAACCTCTCCACAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATACAGCAATAAAAACTGGCAATAATATCACGCAAAATACGCAGCACCGTCTCTTCACCGCTGCGATTCggg                                 >  2:138232/1‑150 (MQ=38)
                                                                                                                    cGTACGTTCCCCTCACCCTAAACCACTACCCCAAGGGGCGAGGGGACCGTTCCGTACGGAAATTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATAATCTCACGCAAAATACGCCCCACCCTCTCTTCACCGCTGCGATTCtgg                                 >  2:2498038/1‑150 (MQ=255)
                                                                                                                       tgggtcaccTCACGCTATAGCTCTCCCCCCAGGGGCGAGGGGACCGTTCAGTACGGATTTTAGCTTGACGGTGCTATCCAGCATTTATTTCTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCAGCGCTTCGATTCGGGTTa                              <  2:3731542/143‑1 (MQ=255)
                                                                                                                        cGTTCCCCTCACCCTAACCCTCTCCCCCCAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTaa                             <  2:2623169/150‑1 (MQ=255)
                                                                                                                         gTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCAc                                               <  2:3310309/131‑1 (MQ=255)
                                                                                                                         gTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCAc                                               >  1:3310309/1‑131 (MQ=255)
                                                                                                                          ttCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCCCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATg                           <  1:1207562/150‑1 (MQ=255)
                                                                                                                             cccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGa                        <  2:3280079/150‑1 (MQ=255)
                                                                                                                              ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                       >  2:834014/1‑150 (MQ=255)
                                                                                                                              ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                       >  2:878109/1‑150 (MQ=255)
                                                                                                                              ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                       >  2:2604593/1‑150 (MQ=255)
                                                                                                                              ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                       >  2:2048284/1‑150 (MQ=255)
                                                                                                                              ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                       >  1:3035549/1‑150 (MQ=255)
                                                                                                                              ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                       >  1:526818/1‑150 (MQ=255)
                                                                                                                               cTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAAt                      <  1:2227238/150‑1 (MQ=255)
                                                                                                                                tcccccTACCCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATc                     <  2:451683/147‑1 (MQ=255)
                                                                                                                                tCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATc                     >  2:1455258/1‑150 (MQ=255)
                                                                                                                                 cACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCg                    <  1:1455258/150‑1 (MQ=255)
                                                                                                                                   cccTAACCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGCATTAACCCGAATCGca                  <  2:1208388/150‑1 (MQ=255)
                                                                                                                                   cccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGca                  <  2:2546226/150‑1 (MQ=255)
                                                                                                                                   cccTAACCCACTCCCCAAAGGGGCGAGGGGACCGTTCAGAACGGAAATTACCTTGACGGAGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGca                  <  1:3973502/150‑1 (MQ=255)
                                                                                                                                    ccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcac                 >  1:2153185/1‑150 (MQ=255)
                                                                                                                                      taccccTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACaa               <  2:656064/147‑1 (MQ=255)
                                                                                                                                      tAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACaa               >  2:1879266/1‑150 (MQ=255)
                                                                                                                                          ccTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGtt           >  2:3770564/1‑150 (MQ=255)
                                                                                                                                          ccTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGtt           >  2:2826530/1‑150 (MQ=255)
                                                                                                                                           ctctCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGTTGCTATCCAGCATTAAGTACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc          >  1:2542801/1‑150 (MQ=255)
                                                                                                                                            tctcCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCt         >  2:265214/1‑150 (MQ=255)
                                                                                                                                            tctcCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCt         <  1:2848108/150‑1 (MQ=255)
                                                                                                                                            tctcCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCt         >  2:2593688/1‑150 (MQ=255)
                                                                                                                                                 ccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtg    >  1:2487271/1‑150 (MQ=255)
                                                                                                                                                 ccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtg    >  2:1847939/1‑150 (MQ=255)
                                                                                                                                                  cAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtgt   <  1:2208756/150‑1 (MQ=255)
                                                                                                                                                  cAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtgt   <  1:1015697/150‑1 (MQ=255)
                                                                                                                                                   aaaGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtgtg  >  1:979742/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTG  >  NC_000913/3504733‑3505029

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: