Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,870,870 T→A 5.5% E145D (GAA→GAT) ‡ truD ← tRNA pseudouridine(13) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,870,8700TA 5.5% 292.9 / 5.0 91E145D (GAA→GAT) ‡truDtRNA pseudouridine(13) synthase
Reads supporting (aligned to +/- strand):  ref base T (48/38);  new base A (2/3);  total (50/41)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

ACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCT  >  NC_000913/2870722‑2871018
                                                                                                                                                     |                                                                                                                                                    
aCCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACgttgtt                                                                                                                                                     >  2:247666/1‑150 (MQ=255)
aCCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACgttgtt                                                                                                                                                     >  1:303158/1‑150 (MQ=255)
 ccAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTc                                                                                                                                                    >  2:525440/1‑150 (MQ=255)
 ccAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTCCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGCCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTc                                                                                                                                                    >  2:492065/1‑150 (MQ=255)
  cAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCa                                                                                                                                                   <  1:759616/150‑1 (MQ=255)
    aaaaCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAAc                                                                                                                                                 <  1:1109869/150‑1 (MQ=255)
     aaaCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACg                                                                                                                                                <  2:378541/150‑1 (MQ=255)
          ccGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGAACAACGTctga                                                                                                                                           <  1:858639/150‑4 (MQ=255)
           cGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATcg                                                                                                                                          <  1:614479/150‑1 (MQ=255)
             tttATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATcgcg                                                                                                                                        >  1:308659/1‑150 (MQ=255)
             tttATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGAACAACGTctgatcg                                                                                                                                        <  1:955015/150‑7 (MQ=255)
             tttATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGAACAACGTctgatcg                                                                                                                                        <  1:960572/150‑7 (MQ=255)
               tATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGAt                                                                                                                                      <  1:912984/150‑1 (MQ=255)
               tATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGAt                                                                                                                                      >  1:451997/1‑150 (MQ=255)
                 ttGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTg                                                                                                                                    <  1:891629/150‑1 (MQ=255)
                   gCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCt                                                                                                                                  >  1:606663/1‑150 (MQ=255)
                                ggAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGa                                                                                                                     >  2:912694/1‑150 (MQ=255)
                                ggAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGa                                                                                                                     <  1:247666/150‑1 (MQ=255)
                                ggAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGa                                                                                                                     >  2:1265633/1‑150 (MQ=255)
                                ggAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGa                                                                                                                     >  2:1270949/1‑150 (MQ=255)
                                 gAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGaa                                                                                                                    <  2:606663/150‑1 (MQ=255)
                                   gTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAAcc                                                                                                                  >  1:191233/1‑150 (MQ=255)
                                        ggTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAAcc                                                                                                                  >  1:756547/1‑145 (MQ=255)
                                        ggTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAAcc                                                                                                                  <  2:756547/145‑1 (MQ=255)
                                          tttGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTaa                                                                                                           >  2:1345042/1‑150 (MQ=255)
                                             gcgcCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAgg                                                                                                        <  2:575738/150‑1 (MQ=255)
                                             gcgcCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAgg                                                                                                        >  2:594263/1‑150 (MQ=255)
                                              cgcCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTATTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGc                                                                                                       >  1:128700/1‑150 (MQ=255)
                                               gcCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCg                                                                                                      >  2:313693/1‑150 (MQ=255)
                                                 ccAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGtt                                                                                                    <  2:38050/150‑1 (MQ=255)
                                                  cAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTa                                                                                                   >  2:195324/1‑150 (MQ=255)
                                                          cgcCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGGTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAg                                                                                           <  1:1173434/150‑1 (MQ=255)
                                                            cccccTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAgcg                                                                                         <  2:1196596/150‑1 (MQ=255)
                                                               ccTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAgcgcgc                                                                                      >  1:73874/1‑150 (MQ=255)
                                                                  gCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTa                                                                                   >  2:745888/1‑150 (MQ=255)
                                                                  gCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTa                                                                                   >  1:849527/1‑150 (MQ=255)
                                                                  gCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTa                                                                                   <  2:1093311/150‑1 (MQ=255)
                                                                   cAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTaa                                                                                  >  1:513309/1‑150 (MQ=255)
                                                                     aGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAAc                                                                                >  1:1451777/1‑150 (MQ=255)
                                                                      gTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACg                                                                               <  2:337487/150‑1 (MQ=255)
                                                                        tGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCa                                                                             <  2:334054/150‑1 (MQ=255)
                                                                        tGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCa                                                                             <  2:339530/150‑1 (MQ=255)
                                                                          cTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGc                                                                           <  1:247282/150‑1 (MQ=255)
                                                                           tACCGCCAATCCCAAAACGTTTGTCAACTAAGTAGTTCGGTACACCTTTTACGCAAATATCCATCAGACTTTGTTCAACGTCATCGCTATTGCTCACTTCGCTCAGAACCATGGTAAAGTCGTTACCTTTCAGCGCGCCTAAACGCAGCt                                                                          <  2:265684/150‑1 (MQ=255)
                                                                              cGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTgc                                                                       >  2:633210/1‑150 (MQ=255)
                                                                                ccAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTgcgc                                                                     <  1:1345042/150‑1 (MQ=255)
                                                                                ccAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTgcgc                                                                     <  2:1285757/150‑1 (MQ=255)
                                                                                 cAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCt                                                                    >  2:329058/1‑150 (MQ=255)
                                                                                     cccAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTtgtg                                                                <  1:487146/150‑1 (MQ=255)
                                                                                     cccAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTtgtg                                                                >  2:1100765/1‑150 (MQ=255)
                                                                                      ccAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAgcgcg                                                                                       <  1:1219526/126‑1 (MQ=255)
                                                                                      ccAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAgcgcg                                                                                       >  2:1219526/1‑126 (MQ=255)
                                                                                      ccAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGc                                                               >  1:1463900/1‑150 (MQ=255)
                                                                                          aaCGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCTCTTTGTGTTTATGATGGGTATAGGCGTTAGCTTTCAGCGCGCCTAAACGCATCTTGTGCTTGTGGcgtt                                                           <  2:303158/150‑3 (MQ=255)
                                                                                           aCGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCgcgc                                                          <  2:619229/150‑1 (MQ=255)
                                                                                              ttGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATa                                                       <  1:1462494/150‑1 (MQ=255)
                                                                                                gggCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACt                                                     >  1:649843/1‑150 (MQ=255)
                                                                                                 ggCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTc                                                    >  1:59906/1‑150 (MQ=255)
                                                                                                   cACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCa                                                  >  1:645228/1‑150 (MQ=255)
                                                                                                   cACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCa                                                  >  2:1059849/1‑150 (MQ=255)
                                                                                                      cGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCa                                               >  1:898646/1‑150 (MQ=255)
                                                                                                      cGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCa                                               <  1:745888/150‑1 (MQ=255)
                                                                                                         agtagtTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCt                                            >  1:1320348/1‑150 (MQ=255)
                                                                                                             gtTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCa                                        >  2:1474447/1‑150 (MQ=255)
                                                                                                              ttCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAg                                       <  2:980942/150‑1 (MQ=255)
                                                                                                                cGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGcc                                     <  2:1451777/150‑1 (MQ=255)
                                                                                                                cGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGcc                                     >  1:389393/1‑150 (MQ=255)
                                                                                                                cGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGcc                                     >  1:395354/1‑150 (MQ=255)
                                                                                                                 ggTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCt                                    >  2:485191/1‑150 (MQ=255)
                                                                                                                  gTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCtt                                   >  1:535190/1‑150 (MQ=255)
                                                                                                                     cacCTTTTACGCAAAATATCGATGAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTCGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTcc                                 <  2:789762/150‑1 (MQ=255)
                                                                                                                      acCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAg                               <  1:735741/150‑1 (MQ=255)
                                                                                                                        cTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCTCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCt                                    <  2:1486834/143‑1 (MQ=255)
                                                                                                                        cTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCt                                    >  1:1486834/1‑143 (MQ=255)
                                                                                                                             aCGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAg                        <  1:485191/150‑1 (MQ=255)
                                                                                                                             aCGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAg                        >  1:746935/1‑150 (MQ=255)
                                                                                                                                cAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGgcg                     >  1:801693/1‑150 (MQ=255)
                                                                                                                                cAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGgcg                     >  1:524731/1‑150 (MQ=255)
                                                                                                                                cAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCCTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGgcc                     >  1:1370636/1‑149 (MQ=255)
                                                                                                                                    atatCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTc                  <  1:746766/150‑1 (MQ=255)
                                                                                                                                    atatCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTc                  <  2:1068720/150‑1 (MQ=255)
                                                                                                                                     tatCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCAATTCGCGCAGAACCAGGGTAAATGCGTTACCTTTCAGAGCGCCTAAACGCAGCTTGCGCTTGTGACGCGCATACTCCAGCACCTGGCAGACTTCCAGTTGAAAGGCGCTCa                 >  2:915163/1‑150 (MQ=255)
                                                                                                                                      atCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTTCAACTGGAAGgctgccag                <  1:156849/150‑7 (MQ=255)
                                                                                                                                        cGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGAt              <  2:898646/150‑1 (MQ=255)
                                                                                                                                         gcgatgACGTTGAACAACGTCATCGCGATTTCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATc             >  2:644625/6‑150 (MQ=255)
                                                                                                                                         gATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATc             >  2:822700/1‑150 (MQ=255)
                                                                                                                                         gATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATc             >  1:1051118/1‑150 (MQ=255)
                                                                                                                                          aTCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCACAGATCg            >  2:251524/1‑150 (MQ=255)
                                                                                                                                            cAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCggg          >  2:578172/1‑150 (MQ=255)
                                                                                                                                              gACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCa        >  1:800350/1‑150 (MQ=255)
                                                                                                                                              gACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCa        >  2:876635/1‑150 (MQ=255)
                                                                                                                                              gACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCa        <  1:736835/150‑1 (MQ=255)
                                                                                                                                                cGTTGAACAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCAtt      >  2:648564/1‑150 (MQ=255)
                                                                                                                                                   tgttCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTcc   >  1:807929/1‑150 (MQ=255)
                                                                                                                                                    gttCAACTTCATCGCGATTGCTCAGTTCGCTCAGAACCAGGTTAAAGGCGTTAGCTTTCAGCGCTCCTAAACGCATCTTGCGCTTGTGCCGCGCATACTACAGCACGTGGCAGCCTTCCATTTGAAAGGCGCTCATATCGGGCATTTCCt  <  2:557637/150‑1 (MQ=255)
                                                                                                                                                    gttCAACGTCATCGCGATTGCTCAGTTCGCGCAGAACCAGGTTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGAGCTTGTGGCGCGCATACTCCAGCACCTGGCAGCCTTCGAGTTGAAAGGCGCTCAGATCGGGCATTTCCt  <  2:550225/150‑1 (MQ=255)
                                                                                                                                                    gttCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCt  <  1:492065/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCT  >  NC_000913/2870722‑2871018

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: