Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,441,822 T→A 5.0% V20D (GTT→GAT)  mnmC → fused 5‑methylaminomethyl‑2‑thiouridine‑forming enzyme methyltransferase and FAD‑dependent demodification enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,441,8220TA 5.0% 309.7 / 3.9 97V20D (GTT→GAT) mnmCfused 5‑methylaminomethyl‑2‑thiouridine‑forming enzyme methyltransferase and FAD‑dependent demodification enzyme
Reads supporting (aligned to +/- strand):  ref base T (45/47);  new base A (2/3);  total (47/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

AAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGTTTGTGGTAGCAGAGAGCGGCTTCGG  >  NC_000913/2441675‑2441966
                                                                                                                                                   |                                                                                                                                                
aaGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACAGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGtt                                                                                                                                                 >  6:211799/1‑149 (MQ=255)
   tCCGATCAGCCAATTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCAATACACCCTGCCAACCTCGAATTTAATGCTGAGGGTACACGTGTTTcc                                                                                                                                              >  4:63997/1‑149 (MQ=255)
      gaccacccAGTTAATAGCTAAATTTGCCCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCga                                                                                                                                           <  2:324252/143‑1 (MQ=255)
      gATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTccc                                                                                                                                             >  4:100793/1‑147 (MQ=255)
        tCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCgaga                                                                                                                                         <  6:68727/149‑1 (MQ=255)
          aGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGt                                                                                                                                                  >  1:316633/1‑138 (MQ=255)
          aGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGt                                                                                                                                                  <  2:316633/138‑1 (MQ=255)
           gCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGAttt                                                                                                                                      >  3:179097/1‑149 (MQ=255)
            ccATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAACCACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGAtttt                                                                                                                                     >  4:331172/1‑149 (MQ=255)
            ccATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGAtttt                                                                                                                                     <  1:9628/149‑1 (MQ=255)
             cATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTGAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCACCCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGAtttt                                                                                                                                     >  7:197519/1‑148 (MQ=255)
              aTTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGATACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGa                                                                                                                                   >  8:79278/1‑149 (MQ=255)
               tttAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGAc                                                                                                                                  <  3:198660/149‑1 (MQ=255)
               tttAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACCTGAAACACTACTCCATACAACCGGCCAACCTCGAATTTAATGCTGACGGTACACCTGTTTCCCGAGATTTTGAc                                                                                                                                  >  3:39737/1‑149 (MQ=255)
                   aTAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGAt                                                                                                                               <  8:155176/148‑1 (MQ=255)
                   aTAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGAt                                                                                                                               <  8:256819/148‑1 (MQ=255)
                   aTAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACAGCCTGTAAGACGAGTAACTTACGTGAAACACTAATCCAAACAACCTGCCAAACACGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATg                                                                                                                              >  2:364480/1‑149 (MQ=255)
                     agagAAATTTGGGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTc                                                                                                                            >  2:108484/1‑149 (MQ=255)
                      gagaAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTc                                                                                                                            >  5:80588/1‑148 (MQ=255)
                       agaAATTTGGGCAGCCTTACACACATCGCTAAGATCGAGCCTCCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTa                                                                                                                          <  1:108484/149‑1 (MQ=255)
                         aaaTTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTAtt                                                                                                                        >  1:284894/1‑149 (MQ=255)
                                cgcAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCaa                                                                                                                 <  2:263476/149‑1 (MQ=255)
                                 gcAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTACCAACCTCGAATTTAATGCTGAGGGTACACCTGATTCCCGAGATTTTGACGATGTCTATTTTTCCAAc                                                                                                                >  2:262515/1‑149 (MQ=255)
                                       ttACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAAc                                                                                                          >  5:99116/1‑149 (MQ=255)
                                           acacaTCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGc                                                                                                      <  4:336417/149‑1 (MQ=255)
                                            cacaTCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCt                                                                                                     <  2:284894/149‑1 (MQ=255)
                                             acaTCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTg                                                                                                    <  6:80588/149‑1 (MQ=255)
                                             acaTCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTg                                                                                                    >  3:108231/1‑149 (MQ=255)
                                             acaTCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTg                                                                                                    >  1:116007/1‑149 (MQ=255)
                                               aTCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTgg                                                                                                   >  7:108040/1‑148 (MQ=255)
                                                tCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGaa                                                                                                 <  7:155422/149‑1 (MQ=255)
                                                    tAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGt                                                                                                                             >  6:18647/1‑117 (MQ=255)
                                                    tAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAgaga                                                                                             >  8:100368/1‑149 (MQ=255)
                                                    tAAGATCGAGCCACCGCCTGCAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATACCGGGCTGGAAgaga                                                                                             >  1:334474/1‑149 (MQ=255)
                                                       gATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGAcgc                                                                                          <  4:66375/149‑1 (MQ=255)
                                                        aTCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACACCCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGaa                                                                                                 >  1:376925/1‑141 (MQ=255)
                                                        aTCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGaa                                                                                                 <  2:376925/141‑1 (MQ=255)
                                                         tCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGt                                                                                        >  1:299584/1‑149 (MQ=255)
                                                           gAGCCACCGCCTGTAAGACGTGTTACTTACGTGAAACACTTCTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGTTGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTa                                                                                      <  3:63997/149‑1 (MQ=255)
                                                            aGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTAt                                                                                     <  3:100793/149‑1 (MQ=255)
                                                             gCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGAcgc                                                                                          <  7:168064/143‑1 (MQ=255)
                                                             gCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGAcgc                                                                                          >  8:168064/1‑143 (MQ=255)
                                                             gCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGATTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTAt                                                                                     >  7:299115/1‑148 (MQ=255)
                                                             gCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGATTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTAt                                                                                     <  8:299115/148‑1 (MQ=255)
                                                                    catGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACATGTTTCCCGAGAATTTGACGATGTCTATTTATCCAACGATAACGGGCGGCAAGAGACGCGTTATGTATTTCt                                                                             >  2:296595/3‑149 (MQ=255)
                                                                     cTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCACCCTCGAAATTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTg                                                                            >  6:268679/1‑149 (MQ=255)
                                                                     cTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTg                                                                            >  8:181711/1‑149 (MQ=255)
                                                                        tAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGTTGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGa                                                                         <  5:154719/149‑1 (MQ=255)
                                                                           gACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGttttt                                                                               >  3:37433/1‑140 (MQ=255)
                                                                           gACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGttttt                                                                               <  4:37433/140‑1 (MQ=255)
                                                                             cGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCAGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCaa                                                                    <  2:116007/149‑1 (MQ=255)
                                                                              gAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAAc                                                                   >  8:129146/1‑149 (MQ=255)
                                                                                    cTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTa                                                             <  7:181711/149‑1 (MQ=255)
                                                                                       aCGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTaga                                                           <  4:179097/148‑1 (MQ=255)
                                                                                       aCGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTaga                                                           <  6:99116/148‑1 (MQ=255)
                                                                                        cGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTagag                                                          <  6:317932/148‑1 (MQ=255)
                                                                                        cGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAgg                                                         <  4:108231/149‑1 (MQ=255)
                                                                                                 ctactaCATACAACCTGCAAACCTCTAAATTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTc                                                <  1:364480/149‑1 (MQ=255)
                                                                                                 ctactCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTc                                                <  7:100368/149‑1 (MQ=255)
                                                                                                 ctactCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTc                                                <  3:331172/149‑1 (MQ=255)
                                                                                                    ctCCATACAACCTACCAACCTCGAATTTAATGCTGAGGGTACACCTGATTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCt                                              <  1:262515/148‑1 (MQ=255)
                                                                                                       cATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGa                                            <  8:108040/147‑1 (MQ=255)
                                                                                                       cATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGa                                            <  7:79278/147‑1 (MQ=255)
                                                                                                       cATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGa                                            <  1:289051/147‑1 (MQ=255)
                                                                                                       cATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGa                                            >  2:289051/1‑147 (MQ=255)
                                                                                                       cATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAg                                           >  5:349778/1‑148 (MQ=255)
                                                                                                        aTACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGtt                                                                                       <  2:322299/103‑1 (MQ=255)
                                                                                                        aTACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGtt                                                                                       >  1:322299/1‑103 (MQ=255)
                                                                                                        aTACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGa                                            >  1:85079/1‑146 (MQ=255)
                                                                                                        aTACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGa                                            <  2:85079/146‑1 (MQ=255)
                                                                                                          aCAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTAc                                                      <  1:248152/134‑1 (MQ=255)
                                                                                                          aCAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTAc                                                      >  2:248152/1‑134 (MQ=255)
                                                                                                           cAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATc                                       <  7:129146/148‑1 (MQ=255)
                                                                                                           cAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATc                                       <  7:339922/148‑1 (MQ=255)
                                                                                                              ccTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCcaca                                   <  5:211799/149‑1 (MQ=255)
                                                                                                               cTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACAt                                  >  6:27612/1‑149 (MQ=255)
                                                                                                                     aCCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGCCGGTGTCTATTTTTCCAACGATAACGGGCTGGAAGAGCCGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTGTCCTGAGCAGCCACCTCCTCTg                            >  6:211905/1‑149 (MQ=255)
                                                                                                                      ccTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATATTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGTGGCAACCAATTAGAGGTACGCTTTCATGAGCATCCACATCCTCTGt                           >  8:271957/1‑149 (MQ=255)
                                                                                                                       cTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGtt                          >  4:18096/1‑149 (MQ=255)
                                                                                                                        tctaaTTTAATGCTGAGGGTACACCTGTTTCCCGATATTTTTACTATGTCTATTTTTCAAAATATAAATTGATTTAAGATACGATTTATGTTTTTCTGGGAGGAAACCAATTAGAGGTAAGCTTTACTGACCATACACATActcttatt                         >  4:233768/4‑145 (MQ=255)
                                                                                                                          gAATTTAATGCTGAGCGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTATCCAACGATAACGGGCTGGACGAGACGCGTTATGTTTTTCTGGGCGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGTTtg                        <  6:47216/148‑1 (MQ=255)
                                                                                                                            aTTTATTGCTGAGTGTACACCTGTTTCCCGAGATGTTGCAGATGTCTATTTTTAGATCGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTTGAGGTACGCTTTCCTGAGCATCCTCATCCTCTGTTTGTgg                     <  4:39737/149‑1 (MQ=255)
                                                                                                                            aTTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGTTTGTgg                     <  6:349778/149‑1 (MQ=255)
                                                                                                                             tttAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATACCGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGTTTGTgg                     >  2:322867/1‑148 (MQ=255)
                                                                                                                             tttAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGTTTGTgg                     >  7:152686/1‑148 (MQ=255)
                                                                                                                                aaTGCTGATGGTACACCTGTTTCCAGAGATTTTGACGATGTCTATTTTTCGAACGATAACGGGCTGGAAGAGAGGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGTTTGTGGTAg                  <  1:296595/148‑1 (MQ=255)
                                                                                                                                 aTGCTGAGTGTACCCCTGATTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGTTTGTGGTAGc                 <  8:152686/148‑1 (MQ=255)
                                                                                                                                    cTGAGGGTACACCTGTTTCCCGAGATTTTGACGGTGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGTTTGTGGTAGCaga              <  5:211905/148‑1 (MQ=255)
                                                                                                                                    cTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGTTTGTGGTAGCa                >  8:244425/1‑146 (MQ=255)
                                                                                                                                              acCTGTTTCCCGAGATTTTGCCGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTATGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGTTTGTGGTAGCAGAGAGCGGCTTCg   >  7:355448/1‑149 (MQ=255)
                                                                                                                                              acCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGTTTGTGGTAGCagagag           <  5:64646/141‑1 (MQ=255)
                                                                                                                                              acCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGTTTGTGGTAGCAGAGAGCGGCTTCg   >  7:229524/1‑149 (MQ=255)
                                                                                                                                              acCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATAAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGTTTGTGGTAGCagagag           >  6:64646/1‑141 (MQ=255)
                                                                                                                                               ccTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGTTTGTGGTAGCAGAGAGCGGCTTCgg  <  5:268679/149‑1 (MQ=255)
                                                                                                                                                cTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTgg                                                                           <  6:315189/75‑1 (MQ=255)
                                                                                                                                                cTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTgg                                                                           >  5:315189/1‑75 (MQ=255)
                                                                                                                                                cTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCt                                                   <  1:348915/99‑1 (MQ=255)
                                                                                                                                                cTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCt                                                   >  2:348915/1‑99 (MQ=255)
                                                                                                                                                 tGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGtt                          >  4:227645/1‑123 (MQ=255)
                                                                                                                                                 tGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGtt                          <  3:227645/123‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
AAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGTTTGTGGTAGCAGAGAGCGGCTTCGG  >  NC_000913/2441675‑2441966

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