Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,147,478 G→T 5.8% A21A (GCC→GCA alkA ← DNA‑3‑methyladenine glycosylase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,147,4780GT 5.8% 304.6 / 3.6 91A21A (GCC→GCAalkADNA‑3‑methyladenine glycosylase 2
Reads supporting (aligned to +/- strand):  ref base G (38/47);  new base T (4/2);  total (42/49)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGC  >  NC_000913/2147329‑2147623
                                                                                                                                                     |                                                                                                                                                 
gCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACgg                                                                                                                                                   >  2:1340881/1‑150 (MQ=255)
    cAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGggcgg                                                                                                                                               <  2:1320583/150‑1 (MQ=255)
         tCTAAACCTGCGCTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACAACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCgaga                                                                                                                                          >  1:742925/1‑150 (MQ=255)
         tCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCGCCGCCCGTGCGGTgagc                                                                                                                                          <  1:209330/150‑2 (MQ=255)
         tCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCgagc                                                                                                                                          <  2:1168791/150‑2 (MQ=255)
          cTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGaa                                                                                                                                         >  2:351260/1‑150 (MQ=255)
            aaaCCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAAt                                                                                                                                       >  2:571654/1‑150 (MQ=255)
               ccTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATccc                                                                                                                                    >  2:824150/1‑150 (MQ=255)
               ccTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATccc                                                                                                                                    >  2:1077570/1‑150 (MQ=255)
                cTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCa                                                                                                                                   <  2:1474535/150‑1 (MQ=255)
                cTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCa                                                                                                                                   <  1:935264/150‑1 (MQ=255)
                   cACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACa                                                                                                                                <  2:1513156/150‑1 (MQ=255)
                    aCTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGCCCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACAt                                                                                                                               >  2:635098/1‑150 (MQ=255)
                       tAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCa                                                                                                                            >  2:163818/1‑150 (MQ=255)
                         aaTTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCCACATCCCCg                                                                                                                          >  2:228181/1‑150 (MQ=255)
                         aaTTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACg                                                                                                                          <  2:1081646/150‑1 (MQ=255)
                            ttATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGAcc                                                                                                                       <  1:163818/150‑1 (MQ=255)
                            ttATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGAcc                                                                                                                       <  1:735933/150‑1 (MQ=255)
                            ttATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGAcc                                                                                                                       <  2:586210/150‑1 (MQ=255)
                                 tgCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCa                                                                                                                  >  1:941881/1‑150 (MQ=255)
                                  gCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCAt                                                                                                                 >  1:234536/1‑150 (MQ=255)
                                     gagTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATAcg                                                                                                              <  2:489476/150‑1 (MQ=255)
                                       gTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATAcggc                                                                                                            >  1:92641/1‑150 (MQ=255)
                                       gTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATAcggc                                                                                                            <  2:1129631/150‑1 (MQ=255)
                                         aTGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATAcggcgg                                                                                                          >  1:748735/1‑150 (MQ=255)
                                            gCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTg                                                                                                       <  1:1160079/150‑1 (MQ=255)
                                              gggCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGcc                                                                                                     <  2:545283/150‑1 (MQ=255)
                                              gggCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGcc                                                                                                     <  1:824150/150‑1 (MQ=255)
                                                gCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCTGCTGCCAg                                                                                                   <  2:537064/150‑1 (MQ=255)
                                                gCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAg                                                                                                   <  2:531389/150‑1 (MQ=255)
                                                gCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCCCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAg                                                                                                   <  2:1192465/150‑1 (MQ=255)
                                                     atCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTg                                                                                                       <  2:1150348/141‑1 (MQ=255)
                                                     atCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTg                                                                                                       >  1:1150348/1‑141 (MQ=255)
                                                       ccGGAATAGCAGTCACCACGCCGCGATATTCGCCCCCCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCAACGCACGGGCGGCGAGAAATACCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAggg                                                                                            >  1:1251513/1‑150 (MQ=255)
                                                       ccGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGTCTGCCAGTTCAggg                                                                                            <  2:888201/150‑1 (MQ=255)
                                                       ccGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAggg                                                                                            <  2:1053037/150‑1 (MQ=255)
                                                       ccGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAggg                                                                                            <  1:1591628/150‑1 (MQ=255)
                                                        cGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGt                                                                                           >  1:1167793/1‑150 (MQ=255)
                                                          gAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGtat                                                                                         >  1:56127/1‑150 (MQ=255)
                                                           aaTAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGtata                                                                                        <  2:232589/150‑1 (MQ=255)
                                                           aaTAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGtata                                                                                        <  2:49011/150‑1 (MQ=255)
                                                               gCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATc                                                                                    <  2:989950/150‑1 (MQ=255)
                                                                cAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCg                                                                                   >  1:230672/1‑150 (MQ=255)
                                                                 aGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGc                                                                                  >  1:1421572/1‑150 (MQ=255)
                                                                  gTCACCACGCGTCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGTCCGTTTCCACACTGCTCACCGCACGGGCGGCTAGAAATCCCAACATCCACGACCAGTCATTCGGCGGCTGCCAGTTCAGTGTATACATCGct                                                                                 <  2:1426229/150‑2 (MQ=255)
                                                                     accacGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCAtct                                                                              >  2:1237337/1‑150 (MQ=255)
                                                                        acGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCt                                                                           <  2:986049/150‑1 (MQ=255)
                                                                         cgccgcGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCtt                                                                          <  2:1162030/150‑1 (MQ=255)
                                                                          gccgcGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTa                                                                         >  1:239915/1‑150 (MQ=255)
                                                                            cgcgATATTCGCCCACCGCCAGACTACGGGCATAATATCTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATGCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTAtt                                                                       <  2:742925/150‑1 (MQ=255)
                                                                              cgATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCa                                                                     <  2:1167793/150‑1 (MQ=255)
                                                                                 tatTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCTCACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATcc                                                                  <  1:1450921/150‑1 (MQ=255)
                                                                                 tatTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATcc                                                                  >  2:443819/1‑150 (MQ=255)
                                                                                    tCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCt                                                               <  1:516157/150‑1 (MQ=255)
                                                                                      gCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATc                                                                                    >  2:52617/1‑127 (MQ=255)
                                                                                      gCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATc                                                                                    <  1:52617/127‑1 (MQ=255)
                                                                                      gCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCttt                                                             >  2:42494/1‑150 (MQ=255)
                                                                                         cACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAg                                                          <  2:9033/150‑1 (MQ=255)
                                                                                         cACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAg                                                          <  1:503312/150‑1 (MQ=255)
                                                                                           ccGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCa                                                        <  1:1340881/150‑1 (MQ=255)
                                                                                                 gACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACg                                                  <  1:180031/150‑1 (MQ=255)
                                                                                                   cTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGTGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCTGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGtt                                                <  1:1077570/150‑1 (MQ=255)
                                                                                                   cTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGtt                                                >  1:164090/1‑150 (MQ=255)
                                                                                                     aCGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTAt                                              >  1:173566/1‑150 (MQ=255)
                                                                                                     aCGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTAt                                              >  2:865338/1‑150 (MQ=255)
                                                                                                      cGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCCCGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTAtt                                             >  2:440296/1‑150 (MQ=255)
                                                                                                       gggCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTc                                            <  2:922306/150‑1 (MQ=255)
                                                                                                           ataataACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGAc                                        >  2:1552949/1‑150 (MQ=255)
                                                                                                                    gTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCgctt                               >  2:449006/1‑150 (MQ=255)
                                                                                                                     tCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCgcttt                              >  2:857429/1‑150 (MQ=255)
                                                                                                                        gcgACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCgctttgct                           >  2:1171524/1‑150 (MQ=255)
                                                                                                                         cgACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCgctttgctt                          <  1:351260/150‑1 (MQ=255)
                                                                                                                            ccGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCGCTTTGCTTTCa                       <  2:371690/150‑1 (MQ=255)
                                                                                                                            ccGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCGCTTTGCTTTCa                       <  2:941881/150‑1 (MQ=255)
                                                                                                                            ccGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCGCTTTGCTTTCa                       <  1:571654/150‑1 (MQ=255)
                                                                                                                            ccGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCGCTTTGCTTTCa                       <  2:173566/150‑1 (MQ=255)
                                                                                                                              gTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCGCTTTGCTTTCata                     <  2:92641/150‑1 (MQ=255)
                                                                                                                                   cacaCTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCGCTTTGCTTTCATATTCCg                >  2:629615/1‑150 (MQ=255)
                                                                                                                                    ggatttcTCGCCGCCCGTGCGGTGAGAAATCCCAACATCCACGCCCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCGCTTTGCTTTCATATTCCgg               >  2:887617/7‑150 (MQ=255)
                                                                                                                                    ggatttcTCGCCGCCCGTGCGGTGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCGCTTTGCTTTCATATTCCgg               >  1:289031/7‑150 (MQ=255)
                                                                                                                                    acacTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCGCTTTGCTTTCATATTCCgg               <  1:635098/150‑1 (MQ=255)
                                                                                                                                     cacTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCGCTTTGCTTTCATATTCCGGt              <  1:1237337/150‑1 (MQ=255)
                                                                                                                                      atttcTCGCCGCCCGTGCGGTGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCGCTTTGCTTTCATATTCCGGtt             >  1:664519/5‑150 (MQ=255)
                                                                                                                                         gCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCgctttgcttt                         >  2:386925/1‑135 (MQ=255)
                                                                                                                                         gCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCgctttgcttt                         <  1:386925/135‑1 (MQ=255)
                                                                                                                                          cTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCGCTTTGCTTTCATATTCCGGTTGTcg         >  1:1244671/1‑150 (MQ=255)
                                                                                                                                          cTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCGCTTTGCTTTCATATTCCGGTTGTcg         >  1:796831/1‑150 (MQ=255)
                                                                                                                                            cgccGCCCGTGCGGTGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCGCTTTGCTTTCATATTCCGGTTGTcgcg       >  1:842396/3‑150 (MQ=255)
                                                                                                                                            cACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCGCTTTGCTTTCATATTCCGGTTGTcgcg       >  2:822052/1‑150 (MQ=255)
                                                                                                                                            cACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCGCTTTGCTTTCATATTCCGGTTGTcgcg       >  2:548075/1‑150 (MQ=255)
                                                                                                                                              ccGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGAc     <  2:748735/150‑1 (MQ=255)
                                                                                                                                                gCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACgg   >  2:352228/1‑150 (MQ=255)
                                                                                                                                                 cACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCGCTTTGCTTTCAAATTCCGGTTGTCGCGACGGc  <  1:42494/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTCAGCATAAACGTTATTCAGACGCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGC  >  NC_000913/2147329‑2147623

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