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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,870,8710TA 5.4% 291.8 / 5.1 92E145V (GAA→GTA) ‡truDtRNA pseudouridine(13) synthase
Reads supporting (aligned to +/- strand):  ref base T (38/49);  new base A (3/2);  total (41/51)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

CCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCT  >  NC_000913/2870723‑2871018
                                                                                                                                                     |                                                                                                                                                   
ccAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTc                                                                                                                                                    >  1:1266282/1‑150 (MQ=255)
 cAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAATTATCGATCAGACGTTGTTCa                                                                                                                                                   >  2:5110/1‑150 (MQ=255)
 cAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCa                                                                                                                                                   >  1:1516757/1‑150 (MQ=255)
 cAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCa                                                                                                                                                   >  1:1086664/1‑150 (MQ=255)
 cAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCa                                                                                                                                                   <  2:280957/150‑1 (MQ=255)
    aaaCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACg                                                                                                                                                <  1:487008/150‑1 (MQ=255)
     aaCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGt                                                                                                                                               <  2:568042/150‑1 (MQ=255)
      aCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTc                                                                                                                                              <  1:1630112/150‑1 (MQ=255)
        tCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCAt                                                                                                                                            >  1:110963/1‑150 (MQ=255)
        tCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGAACAACGTctg                                                                                                                                            <  2:959942/150‑3 (MQ=255)
          cGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATcg                                                                                                                                          >  2:529774/1‑150 (MQ=255)
           gTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATcgc                                                                                                                                         <  1:1147818/150‑1 (MQ=255)
           gTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATcgc                                                                                                                                         <  2:316768/150‑1 (MQ=255)
            tttATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATcgcg                                                                                                                                        <  1:1136678/150‑1 (MQ=255)
            tttATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGAACAACGTctgatcg                                                                                                                                        <  1:446871/150‑7 (MQ=255)
             ttATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGa                                                                                                                                       >  1:1201076/1‑150 (MQ=255)
                  gCGATCGCTCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACAGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCt                                                                                                                                  <  1:1021436/150‑1 (MQ=255)
                  gCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCt                                                                                                                                  >  1:223701/1‑150 (MQ=255)
                  gCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTCACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCt                                                                                                                                  >  1:745364/1‑150 (MQ=255)
                   cGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTc                                                                                                                                 <  2:1389243/150‑1 (MQ=255)
                      tCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACt                                                                                                                              <  1:1008676/150‑1 (MQ=255)
                       cgcgCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACtt                                                                                                                             <  2:975345/150‑1 (MQ=255)
                           cACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTcgcg                                                                                                                         <  1:341364/150‑1 (MQ=255)
                            aCCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCgcgc                                                                                                                        <  2:1121080/150‑1 (MQ=255)
                             ccGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCa                                                                                                                       <  1:707922/150‑1 (MQ=255)
                               ggAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGa                                                                                                                     <  1:1530708/150‑1 (MQ=255)
                                      tGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAggg                                                                                                              <  1:649005/150‑1 (MQ=255)
                                      gggTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACTGTGGGCACCGAAGGAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAggg                                                                                                              <  2:1631425/149‑1 (MQ=255)
                                       ggTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGt                                                                                                             >  1:1360104/1‑150 (MQ=255)
                                            gcgcCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAgg                                                                                                        <  2:1084820/150‑1 (MQ=255)
                                             cgcCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGc                                                                                                       <  1:260621/150‑1 (MQ=255)
                                             cgcCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGc                                                                                                       >  1:551992/1‑150 (MQ=255)
                                               cccAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGt                                                                                                     >  2:1331886/1‑150 (MQ=255)
                                                 cAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAGGCGTTAc                                                                                                  <  2:1193518/150‑1 (MQ=255)
                                                       tGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTc                                                                                             <  2:72460/150‑1 (MQ=255)
                                                        gcgcCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCa                                                                                            <  1:1166884/150‑1 (MQ=255)
                                                        gcgcCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCa                                                                                            <  1:1250002/150‑1 (MQ=255)
                                                         cgcCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAg                                                                                           <  2:55794/150‑1 (MQ=255)
                                                          gcCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAgc                                                                                          <  1:625722/150‑1 (MQ=255)
                                                            ccccTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCCCTTGGCGCAGAACCAGGGTAAGGGCGTTACCTTTCAgcgc                                                                                        >  1:623877/1‑150 (MQ=255)
                                                                 gCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTa                                                                                   <  1:523000/150‑1 (MQ=255)
                                                                  cAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTaa                                                                                  >  2:506296/1‑150 (MQ=255)
                                                                       tGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCa                                                                             >  2:1397959/1‑150 (MQ=255)
                                                                       tGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCa                                                                             >  2:1397964/1‑150 (MQ=255)
                                                                       tGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCa                                                                             <  1:735729/150‑1 (MQ=255)
                                                                           aCCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCtt                                                                         >  1:1440223/1‑150 (MQ=255)
                                                                               ccAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCATTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTgcgc                                                                     <  1:1342732/150‑1 (MQ=255)
                                                                                cAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCt                                                                    <  2:1086394/150‑1 (MQ=255)
                                                                                    cccAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTtgt                                                                 <  2:553289/150‑1 (MQ=255)
                                                                                    cccAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGACCCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTtgtg                                                                >  1:1542128/1‑150 (MQ=255)
                                                                                      cAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGcc                                                              >  2:96403/1‑150 (MQ=255)
                                                                                      cAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGcc                                                              >  2:783548/1‑150 (MQ=255)
                                                                                        aaaCGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTaa                                                                                  >  2:417650/1‑128 (MQ=255)
                                                                                        aaaCGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTaa                                                                                  <  1:417650/128‑1 (MQ=255)
                                                                                         aaCGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCcgcg                                                           <  2:1491495/150‑1 (MQ=255)
                                                                                            gTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCAt                                                        <  2:223701/150‑1 (MQ=255)
                                                                                            gTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCAt                                                        <  1:570885/150‑1 (MQ=255)
                                                                                             ttGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATa                                                       >  2:182008/1‑150 (MQ=255)
                                                                                              ggtgtaCCGAACTACTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATAc                                                      >  1:1159000/6‑150 (MQ=255)
                                                                                               gtgtaCCGAACTACTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACt                                                     >  1:1422357/5‑150 (MQ=255)
                                                                                                ggCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTc                                                    >  2:1278660/1‑150 (MQ=255)
                                                                                                 gCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAGCGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTcc                                                   <  2:795362/150‑1 (MQ=255)
                                                                                                 gCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTcc                                                   <  2:815697/150‑1 (MQ=255)
                                                                                                  cACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCa                                                  >  2:1460157/1‑150 (MQ=255)
                                                                                                       aaGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCAcc                                             <  2:1254260/150‑1 (MQ=255)
                                                                                                         gtagtTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTg                                           <  1:80481/150‑1 (MQ=255)
                                                                                                            gtTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCa                                        >  2:207985/1‑150 (MQ=255)
                                                                                                             ttCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAg                                       <  1:783548/150‑1 (MQ=255)
                                                                                                              gcGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGc                                      <  2:547165/149‑1 (MQ=255)
                                                                                                               cGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGcc                                     >  1:774516/1‑150 (MQ=255)
                                                                                                                   acacCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTcc                                 >  2:153356/1‑150 (MQ=255)
                                                                                                                    cacCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCa                                >  2:94281/1‑150 (MQ=255)
                                                                                                                     acCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCATAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAg                               >  1:603251/1‑150 (MQ=255)
                                                                                                                     acCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAg                               <  2:1440223/150‑1 (MQ=255)
                                                                                                                     acCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAg                               <  2:1086664/150‑1 (MQ=255)
                                                                                                                      ccTTTTACGCAATTATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGt                              <  1:5110/150‑1 (MQ=255)
                                                                                                                      ccTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGt                              <  2:1266282/150‑1 (MQ=255)
                                                                                                                         ttttACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTg                            >  1:399354/1‑150 (MQ=255)
                                                                                                                            tACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGaaa                         <  1:207985/150‑1 (MQ=255)
                                                                                                                             aCGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAg                        >  1:1394211/1‑150 (MQ=255)
                                                                                                                               gCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGgc                      <  2:110963/150‑1 (MQ=255)
                                                                                                                                  aaTATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCt                   >  2:443945/1‑150 (MQ=255)
                                                                                                                                    tatCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCa                 >  2:408001/1‑150 (MQ=255)
                                                                                                                                       cgcgatgACGTTGAACAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGAt              >  2:238627/7‑150 (MQ=255)
                                                                                                                                        gcgatgACGTTGAACAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATc             >  1:423481/6‑150 (MQ=255)
                                                                                                                                         cgatgACGTTGAACAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCg            >  2:1441590/5‑150 (MQ=255)
                                                                                                                                            aGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGc         <  1:1035512/150‑1 (MQ=255)
                                                                                                                                            aGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGc         <  1:506296/150‑1 (MQ=255)
                                                                                                                                            aGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGc         <  2:551992/150‑1 (MQ=255)
                                                                                                                                             gACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCa        >  1:432205/1‑150 (MQ=255)
                                                                                                                                             gACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCa        >  1:477974/1‑150 (MQ=255)
                                                                                                                                               cGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCAtt      >  2:1006379/1‑150 (MQ=255)
                                                                                                                                                   gttCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTCCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCt  >  1:1344988/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCT  >  NC_000913/2870723‑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: