Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,870,871 T→A 6.3% E145V (GAA→GTA) ‡ truD ← tRNA(Glu) pseudouridine(13) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,870,8710TA 6.3% 344.1 / 11.2 95E145V (GAA→GTA) ‡truDtRNA(Glu) pseudouridine(13) synthase
Reads supporting (aligned to +/- strand):  ref base T (46/43);  new base A (3/3);  total (49/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-01

AAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATT  >  NC_000913/2870725‑2871014
                                                                                                                                                  |                                                                                                                                               
aaaaaCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTc                                                                                                                                                >  4:938951/1‑148 (MQ=255)
 aaaaCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCa                                                                                                                                               >  2:979605/1‑148 (MQ=255)
  aaaCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCaa                                                                                                                                              <  3:938951/148‑1 (MQ=255)
    aCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACg                                                                                                                                            >  4:555639/1‑148 (MQ=255)
       ccGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCa                                                                                                                                         >  3:539080/1‑148 (MQ=255)
       ccGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCa                                                                                                                                         >  3:886395/1‑148 (MQ=255)
            tATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATcgcg                                                                                                                                    <  1:979605/148‑1 (MQ=255)
                   aTCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTc                                                                                                                             >  1:272678/1‑148 (MQ=255)
                           ccGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTcgcg                                                                                                                     >  1:210626/1‑148 (MQ=255)
                           ccGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTcgcg                                                                                                                     >  2:1057998/1‑148 (MQ=255)
                                  aTTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTc                                                                                                                                                >  2:95478/1‑114 (MQ=255)
                                  aTTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTc                                                                                                                                                <  1:95478/114‑1 (MQ=255)
                                  aTTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAAcc                                                                                                              <  2:210626/148‑1 (MQ=255)
                                    tGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTc                                                                                                                                                >  1:143959/1‑112 (MQ=255)
                                    tGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTc                                                                                                                                                <  2:143959/112‑1 (MQ=255)
                                    tGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGa                                                                                                                                   <  1:279653/125‑1 (MQ=255)
                                    tGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGa                                                                                                                                   >  2:279653/1‑125 (MQ=255)
                                    tGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAg                                                                                                            <  2:272678/148‑1 (MQ=255)
                                    tGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAg                                                                                                            <  4:886395/148‑1 (MQ=255)
                                    tGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAg                                                                                                            <  4:539080/148‑1 (MQ=255)
                                             cccAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGAACAACGTctg                                                                                                                                        >  4:1056205/1‑109 (MQ=255)
                                             cccAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGAACAACGTctg                                                                                                                                        <  3:1056205/111‑3 (MQ=255)
                                              ccAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGc                                                                                                                               >  4:1027550/1‑119 (MQ=255)
                                              ccAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGc                                                                                                                               <  3:1027550/119‑1 (MQ=255)
                                                       cgcCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTc                                                                                         <  3:555639/148‑1 (MQ=255)
                                                          ccccTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTaaa                                                                                                      <  1:145734/132‑1 (MQ=255)
                                                          ccccTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTaaa                                                                                                      >  2:145734/1‑132 (MQ=255)
                                                           cccTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTc                                                                                         <  2:1066664/144‑1 (MQ=255)
                                                           cccTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTc                                                                                         >  1:1066664/1‑144 (MQ=255)
                                                               gCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGcc                                                                                 >  1:698688/1‑148 (MQ=255)
                                                                  aGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACgttgtt                                                                                                                                                 <  3:603210/81‑1 (MQ=255)
                                                                  aGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACgttgtt                                                                                                                                                 >  4:603210/1‑81 (MQ=255)
                                                                   gTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTaaa                                                                             <  1:1057998/148‑1 (MQ=255)
                                                                    ttGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCg                                                                                                  >  3:538293/1‑126 (MQ=255)
                                                                    ttGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCg                                                                                                  <  4:538293/126‑1 (MQ=255)
                                                                         aCCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCt                                                                                            <  1:239027/127‑1 (MQ=255)
                                                                         aCCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCt                                                                                            >  2:239027/1‑127 (MQ=255)
                                                                             ccAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTgc                                                                   >  3:849114/1‑148 (MQ=255)
                                                                                   ccAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTtgt                                                             <  2:698688/148‑1 (MQ=255)
                                                                                    cAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATcgcg                                                                                                                                    >  3:339712/1‑76 (MQ=255)
                                                                                    cAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATcgcg                                                                                                                                    <  4:339712/76‑1 (MQ=255)
                                                                                    cAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTcg                                                                                                                       <  4:1036830/89‑1 (MQ=255)
                                                                                    cAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTcg                                                                                                                       <  3:33059/89‑1 (MQ=255)
                                                                                    cAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTcg                                                                                                                       >  4:33059/1‑89 (MQ=255)
                                                                                    cAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTcg                                                                                                                       >  3:1036830/1‑89 (MQ=255)
                                                                                     aaaaCGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCa                                                                                                             <  1:1104196/98‑1 (MQ=255)
                                                                                     aaaaCGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCa                                                                                                             >  2:1104196/1‑98 (MQ=255)
                                                                                      aaaCGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCAt                                                                                                                                        <  1:696880/70‑1 (MQ=255)
                                                                                      aaaCGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCAt                                                                                                                                        >  2:696880/1‑70 (MQ=255)
                                                                                      aaaCGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCa                                                                                                             <  1:224799/97‑1 (MQ=255)
                                                                                      aaaCGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTaaa                                                                                                      <  2:130995/104‑1 (MQ=255)
                                                                                      aaaCGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTaaa                                                                                                      >  1:130995/1‑104 (MQ=255)
                                                                                      aaaCGTTGGGCACCGAAGGAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCa                                                                                                             >  2:224799/1‑97 (MQ=255)
                                                                                                  ccGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCttt                                                                                          <  2:845367/104‑1 (MQ=255)
                                                                                                  ccGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCttt                                                                                          >  1:845367/1‑104 (MQ=255)
                                                                                                   cGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTc                                                                                         <  2:1030777/104‑1 (MQ=255)
                                                                                                   cGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTc                                                                                         >  1:1030777/1‑104 (MQ=255)
                                                                                                     aaGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCa                                           <  4:849114/148‑1 (MQ=255)
                                                                                                           ttCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCg                                                                                                  >  1:996936/1‑87 (MQ=255)
                                                                                                           ttCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCg                                                                                                  <  2:996936/87‑1 (MQ=255)
                                                                                                           ttCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTc                                                                                         <  3:15888/96‑1 (MQ=255)
                                                                                                           ttCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTc                                                                                         >  4:15888/1‑96 (MQ=255)
                                                                                                            tCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCa                                    >  3:661567/1‑148 (MQ=255)
                                                                                                            tCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCa                                    >  4:175257/1‑148 (MQ=255)
                                                                                                                tACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGcc                                 <  3:712110/147‑1 (MQ=255)
                                                                                                                tACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGcc                                 >  4:712110/1‑147 (MQ=255)
                                                                                                                 acacCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTtgtg                                                            <  1:732695/119‑1 (MQ=255)
                                                                                                                 acacCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTtgtg                                                            >  2:732695/1‑119 (MQ=255)
                                                                                                                     cTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTATCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAg                           <  4:661567/148‑1 (MQ=255)
                                                                                                                       tttACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTa                                                                                               <  3:869762/78‑1 (MQ=255)
                                                                                                                       tttACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTa                                                                                               >  4:869762/1‑78 (MQ=255)
                                                                                                                       tttACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTAc                                                                                              <  3:860280/79‑1 (MQ=255)
                                                                                                                       tttACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTAc                                                                                              >  4:860280/1‑79 (MQ=255)
                                                                                                                       tttACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAgcgcg                                                                                   >  3:830468/1‑90 (MQ=255)
                                                                                                                       tttACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAgcgcg                                                                                   <  4:830468/90‑1 (MQ=255)
                                                                                                                       tttACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCt                                                                      >  4:1143512/1‑103 (MQ=255)
                                                                                                                       tttACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCt                                                                      <  3:1143512/103‑1 (MQ=255)
                                                                                                                       tttACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCAc                                          <  3:658478/131‑1 (MQ=255)
                                                                                                                       tttACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCAc                                          >  4:658478/1‑131 (MQ=255)
                                                                                                                       tttACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGtc                         >  3:97210/1‑147 (MQ=255)
                                                                                                                             cAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTgc                                                                   >  4:1061155/1‑100 (MQ=255)
                                                                                                                             cAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTgc                                                                   <  3:1061155/100‑1 (MQ=255)
                                                                                                                               aaTATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGcc                                                                                 <  4:105028/84‑1 (MQ=255)
                                                                                                                               aaTATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGcc                                                                                 >  3:105028/1‑84 (MQ=255)
                                                                                                                                 tatCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCTCCTAAATGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCAcc                                         >  4:689882/1‑122 (MQ=255)
                                                                                                                                 tatCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCAcc                                         <  3:689882/122‑1 (MQ=255)
                                                                                                                                       tCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCa                                                                                                             <  2:601882/48‑1 (MQ=255)
                                                                                                                                       tCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCa                                                                                                             >  1:601882/1‑48 (MQ=255)
                                                                                                                                         tgACGTTGAACAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGcc                                                          >  1:652690/2‑97 (MQ=255)
                                                                                                                                         tgACGTTGAACAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGcc                                                          <  2:652690/96‑1 (MQ=255)
                                                                                                                                         tgACGTTGAACAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCAcc                                         <  3:523894/113‑1 (MQ=255)
                                                                                                                                         tgACGTTGAACAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCAcc                                         >  4:523894/2‑114 (MQ=255)
                                                                                                                                          gACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGc                                                           >  1:828362/1‑95 (MQ=255)
                                                                                                                                          gACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGc                                                           <  2:828362/95‑1 (MQ=255)
                                                                                                                                          gACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGaaa                     <  3:1098011/133‑1 (MQ=255)
                                                                                                                                          gACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGaaa                     >  4:1098011/1‑133 (MQ=255)
                                                                                                                                            cGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCa    >  3:489402/1‑148 (MQ=255)
                                                                                                                                              ttgttCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCAtt  <  4:97210/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
AAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATT  >  NC_000913/2870725‑2871014

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: