Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,423,620 A→G 6.0% noncoding (36/76 nt) thrV ← tRNA‑Thr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,6200AG 6.0% 300.0 / 7.2 99noncoding (36/76 nt)thrVtRNA‑Thr
Reads supporting (aligned to +/- strand):  ref base A (43/50);  new base G (4/2);  total (47/52)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

CGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACC  >  NC_000913/3423474‑3423768
                                                                                                                                                  |                                                                                                                                                    
cgcTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCa                                                                                                                                                   <  2:1417588/150‑1 (MQ=255)
    aCGGCGTTTCACTTCTTAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTACTCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAggg                                                                                                                                               <  1:66411/150‑1 (MQ=255)
    aCGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAggg                                                                                                                                               >  1:1119143/1‑150 (MQ=255)
    aCGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAggg                                                                                                                                               <  1:2215727/150‑1 (MQ=255)
    aCGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAggg                                                                                                                                               <  2:1261596/150‑1 (MQ=255)
     cGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGt                                                                                                                                              <  1:1953659/150‑1 (MQ=255)
       gCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTgc                                                                                                                                            <  1:1841067/150‑1 (MQ=255)
       gCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTgc                                                                                                                                            <  2:1202347/150‑1 (MQ=255)
        cGTTTCACTACTGAGATCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTgcg                                                                                                                                           <  1:1374783/150‑1 (MQ=255)
         gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTGGTAAGGGTgagg                                                                                                                                          <  2:2005683/150‑4 (MQ=255)
         gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTgcgc                                                                                                                                          <  2:1203698/150‑1 (MQ=255)
         gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTgcgc                                                                                                                                          <  2:261741/150‑1 (MQ=255)
         gTTTCACTTCTGAGTTCGGCATGGGGTCAGGGGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTgcgc                                                                                                                                          <  2:1718402/150‑1 (MQ=255)
         gTTTCACTTCTGAGTGCGGCATGGGGGCAGGGGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGGGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTgcgc                                                                                                                                          <  2:360392/150‑1 (MQ=255)
            tCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAATAGTGGTGCTGATACCCAGAGTCGAAGTGGGGACCTCACCCTTACCAAGGGTGCGctct                                                                                                                                       <  1:807579/150‑1 (MQ=255)
            tCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTGGTAAGGGTgaggtcc                                                                                                                                       <  2:2152542/150‑7 (MQ=37)
            tCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACAGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTgaggtcc                                                                                                                                       <  2:1650965/150‑7 (MQ=255)
             cACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTa                                                                                                                                      <  1:1367997/150‑1 (MQ=255)
              aCTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTAc                                                                                                                                     <  2:582653/150‑1 (MQ=255)
               cTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTAcc                                                                                                                                    <  1:2025033/150‑1 (MQ=255)
               cTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTAcc                                                                                                                                    <  1:253750/150‑1 (MQ=255)
               cTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTAcc                                                                                                                                    <  1:607294/150‑1 (MQ=255)
               cTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTAcc                                                                                                                                    <  1:732706/150‑1 (MQ=255)
                    gAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTg                                                                                                                               <  2:2190449/150‑1 (MQ=255)
                    cagatCGGCAGGGGGTCAGGCGGGACCACCGCGCTACTGCCTCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGGGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTg                                                                                                                               <  1:1352963/146‑1 (MQ=38)
                     aGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGa                                                                                                                              <  2:1875265/150‑1 (MQ=255)
                      gTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAg                                                                                                                             <  1:1364652/150‑1 (MQ=255)
                       ttCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGc                                                                                                                            <  2:1245801/150‑1 (MQ=255)
                         cGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCa                                                                                                                          <  1:2128358/150‑1 (MQ=255)
                         cGGCATGGGGGCAGGTGGGACCACCGCGCTACTGCCGCCATACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCa                                                                                                                          <  1:1725468/150‑1 (MQ=255)
                           gCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCata                                                                                                                        <  2:197330/150‑1 (MQ=255)
                             aTGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCAGCCGAACTTTAACCTAAAAAGTGGTGCGGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATc                                                                                                                      <  1:25674/150‑1 (MQ=255)
                             aTGGGGTCAGGTGGGACAACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATc                                                                                                                      <  1:1123428/150‑1 (MQ=255)
                              tGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCa                                                                                                                     >  1:1482884/1‑150 (MQ=255)
                               ggggTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAg                                                                                                                    <  1:2062125/150‑1 (MQ=255)
                                        tGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTaaa                                                                                                           >  2:851835/1‑150 (MQ=255)
                                           gACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAAttt                                                                                                        <  1:468179/150‑1 (MQ=255)
                                             ccaccGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGa                                                                                                      <  1:832388/150‑1 (MQ=255)
                                               accGCGCTACTTCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACATAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACATCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATg                                                                                                    <  2:2181438/150‑1 (MQ=255)
                                               accGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATg                                                                                                    <  1:1675597/150‑1 (MQ=255)
                                               accGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATg                                                                                                    >  2:1481488/1‑150 (MQ=255)
                                                   ggcTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTg                                                                                                <  2:1885284/149‑1 (MQ=255)
                                                   cgcTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTg                                                                                                <  1:1883263/150‑1 (MQ=255)
                                                   cgcTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTg                                                                                                <  2:1407987/150‑1 (MQ=255)
                                                     cTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGc                                                                                              <  1:599130/150‑1 (MQ=255)
                                                       aCTGCCGCCAGACAAATTCTTTTCTAATCTGCCTAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAg                                                                                            <  1:1897392/150‑1 (MQ=255)
                                                       aCTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGTCAg                                                                                            >  1:667030/1‑150 (MQ=255)
                                                          gccgccAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCGACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTc                                                                                         >  2:1614933/1‑150 (MQ=255)
                                                            cgccACACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTccc                                                                                       <  1:1391470/150‑1 (MQ=255)
                                                             gccAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCt                                                                                      <  2:1119143/150‑1 (MQ=255)
                                                                   cAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACtctc                                                                                >  2:629070/1‑150 (MQ=255)
                                                                        tCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCCCCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCCGCACGCGAAATTTGATGCCTGGCAGTTACCTACTCTCGCATg                                                                           >  2:166732/1‑150 (MQ=255)
                                                                         cTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATgg                                                                          >  1:1954893/1‑150 (MQ=255)
                                                                         cTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATgg                                                                          >  1:564340/1‑150 (MQ=255)
                                                                           tttCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATgggg                                                                        >  1:1049750/1‑150 (MQ=255)
                                                                            ttCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGga                                                                       <  1:851835/150‑1 (MQ=255)
                                                                                 aTCTGCCGAACTTTAACCTAAACAGTGGGGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCGCTACCAACTGAGCCATATCAGCACGCTAAATTGGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGAccc                                                                  <  1:1614933/150‑1 (MQ=255)
                                                                                 aTCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGAccc                                                                  >  1:665501/1‑150 (MQ=255)
                                                                                 aTCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGAccc                                                                  >  1:194579/1‑150 (MQ=255)
                                                                                 aTCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGAccc                                                                  >  1:216545/1‑150 (MQ=255)
                                                                                  tCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATgggg                                                                        >  1:765164/1‑143 (MQ=255)
                                                                                  tCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATgggg                                                                        <  2:765164/143‑1 (MQ=255)
                                                                                   cTGCCGAACTTTAACCTAAAAAGTGGAGCTGATACCCAGAGTCGAACTGGGGACCTCACACTTACCAAGGGTGCCCTCTACAAACTGAGCCATATCAGCACCCTAAATTTTATGCCTGGCAGTTCCCTACTCTCGCATGTGGATACCCca                                                                >  2:1216490/1‑150 (MQ=38)
                                                                                      ccGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCacac                                                             >  2:1464482/1‑150 (MQ=255)
                                                                                       cGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACt                                                            <  2:1482884/150‑1 (MQ=255)
                                                                                        gAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCATAACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTa                                                           >  1:409329/1‑150 (MQ=255)
                                                                                           cTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCCCACTACCa                                                        >  2:1892856/1‑150 (MQ=255)
                                                                                           cTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCa                                                        >  1:377053/1‑150 (MQ=255)
                                                                                             ttAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCGGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATc                                                      >  1:222920/1‑150 (MQ=255)
                                                                                                  cTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCCCGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGgcgc                                                 >  1:369448/1‑150 (MQ=255)
                                                                                                   tAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCt                                                >  1:1552668/1‑150 (MQ=255)
                                                                                                    aaaaaGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTa                                               >  2:1585024/1‑150 (MQ=255)
                                                                                                             tgCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTc                                      >  2:1274410/1‑150 (MQ=255)
                                                                                                              gCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCa                                     >  1:256638/1‑150 (MQ=255)
                                                                                                                tGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAAGTGCGCCATATCAGCCCGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCCTGGGGAGACCCCACACTACCATCGGCGCTACGGCGTGTCCCt                                   >  1:324858/1‑150 (MQ=25)
                                                                                                                tGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACt                                   >  1:2003106/1‑150 (MQ=255)
                                                                                                                tGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACt                                   >  1:588818/1‑150 (MQ=255)
                                                                                                                 gATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACtt                                  >  2:980915/1‑150 (MQ=255)
                                                                                                                 gATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCGAAATTTGATGGCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACAATCGGACCTACGGCGTTTCACtt                                  >  2:1809230/1‑150 (MQ=37)
                                                                                                                    actcAGAGTCGATCTGTCTACCTCACCTATACCAATCCTGCGCTCTACCTCCTGCGCCATATCCGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTg                               <  1:1980689/147‑1 (MQ=17)
                                                                                                                     cccAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGa                              >  1:1625350/1‑150 (MQ=255)
                                                                                                                      ccAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAg                             >  2:930799/1‑150 (MQ=255)
                                                                                                                      ccAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAg                             >  2:2155584/1‑150 (MQ=255)
                                                                                                                           gTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCgg                        <  1:1481488/150‑1 (MQ=255)
                                                                                                                              gAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCAt                     >  2:1174170/1‑150 (MQ=255)
                                                                                                                              gAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCGACGGCGTTTCACTTCTGAGTTCGGCAt                     >  1:1614379/1‑150 (MQ=255)
                                                                                                                                  tGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGCCCCCACACTACCATCGGCGCTACGGCGTTTCCCTTCTGAGTTCGGCATgggg                 >  1:2046050/1‑150 (MQ=255)
                                                                                                                                   ggggACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGt                >  2:1525630/1‑150 (MQ=255)
                                                                                                                                   ggggACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGt                >  1:695475/1‑150 (MQ=255)
                                                                                                                                      gagcgcACCCTTGGTAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAg             >  2:688654/6‑150 (MQ=255)
                                                                                                                                      gACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGGACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACACCTCACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAg             <  2:409329/150‑1 (MQ=255)
                                                                                                                                       agcgcACCCTTGGTAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAgg            >  1:1083975/5‑150 (MQ=255)
                                                                                                                                        ccTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGt           >  2:1900461/1‑150 (MQ=255)
                                                                                                                                         cgcACCCTTGGTAAGGGTGCGCTCTCCCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCGACTCTCGCATGGGGAGCCCCCACACTACCCTCGGCGCTACGGCGTTTCCCTTCTGAGGTCGGCATGGGGTCAGGTg          >  1:1478944/3‑150 (MQ=18)
                                                                                                                                         cgcACCCTTGGTAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACCTCTGAGTTCGGCATGGGGTCTGGTg          >  1:319036/3‑150 (MQ=18)
                                                                                                                                         cTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGTCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTg          <  2:667030/150‑1 (MQ=255)
                                                                                                                                           cACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTggg        >  1:2108645/1‑150 (MQ=255)
                                                                                                                                             cccTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGac      >  1:596055/1‑150 (MQ=255)
                                                                                                                                                 tACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGaccacc  >  1:1286962/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACC  >  NC_000913/3423474‑3423768

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