Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,863,258 A→T 7.4% N163I (AAC→ATC) ‡ gapA → glyceraldehyde‑3‑phosphate dehydrogenase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,863,2580AT 7.4% 197.6 / 5.1 77N163I (AAC→ATC) ‡gapAglyceraldehyde‑3‑phosphate dehydrogenase A
Reads supporting (aligned to +/- strand):  ref base A (38/32);  new base T (2/4);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

TGCGAAGAAAGTGGTTATGACTGGTCCGTCTAAAGACAACACTCCGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCT  >  NC_000913/1863109‑1863406
                                                                                                                                                     |                                                                                                                                                    
tGCGAAGAAAGTGGTTATGACTGGTCCGTCTAAAGACAACACTCCGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCaa                                                                                                                                                      >  2:972085/1‑150 (MQ=255)
tGCGAAGAAAGTGGATATGACTGGTCCGGCTAAAGACAACACTCCGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCaa                                                                                                                                                      <  1:103213/150‑1 (MQ=255)
 gCGAAGAAAGTGGTTATGACTGGACCGTCTAAAGACAACACTACGATGTTCGTTAAAGGCGCTAACATCGACAAAAATGATGGCCAGGACATCGTAACCAACGCTTACTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTAATCAAc                                                                                                                                                     >  1:1458784/1‑150 (MQ=255)
  cGAAGAAAGTGGTTATGACTGGTCTGTCTAAAGACAACACTCCGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTCCACCCCCAACTGCCTGGCTCCGCTGGCTAAAGTTTTCAACg                                                                                                                                                    >  2:328431/1‑150 (MQ=255)
   gaagaaAGTGGTTATGACTGGTCCGTCTAAAGACAACACTCCGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGa                                                                                                                                                   >  1:347721/1‑150 (MQ=255)
          gTGGTTATGACTGGTCCGTCTAAAGACAACACTCCGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCGTTGATAACttt                                                                                                                                            <  1:241448/150‑2 (MQ=255)
           tGGTTATGACTGGTCCGTCTAAAGACAACACTCCGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCGTTGATAACttta                                                                                                                                           <  2:767104/150‑3 (MQ=255)
           tGGTTATGACTGGTCCGTCTAAAGACAACACTCCGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCGTTGATAACttta                                                                                                                                           <  2:1288572/150‑3 (MQ=255)
           tGGTTATGACTGGTCCGTCTAAAGACAACACTCCGATGTTCGTTAAAGGCGCTAACTTCGACAAATATACTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCg                                                                                                                                           >  1:1666775/1‑150 (MQ=255)
            ggatATGACTGGTACGTCTAAAGACAACACTCCGATGTTCGTTAAATGCGCTAACTTCGACAAATAAGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTCGCTCCGCTGGCTAAAGTTATCGTTGATAACtttag                                                                                                                                          <  1:250520/147‑4 (MQ=255)
            ggTTATGACTGGTCCGTCTAAAGACAACACTCCGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCgg                                                                                                                                          <  1:183043/150‑1 (MQ=255)
              ttATGACTGGTCCGTCTAAAGACAACACTCCGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGca                                                                                                                                        <  1:1511390/150‑1 (MQ=255)
                aTGACTGGTCCGTCTAAAGACAACACTCCGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGcatc                                                                                                                                      >  2:878681/1‑150 (MQ=255)
                          cGTCTAAAGACAACACTCCGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTc                                                                                                                            <  1:651160/150‑1 (MQ=255)
                            tCTAAAGACAACACTCCGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCtg                                                                                                                          <  2:361287/150‑1 (MQ=255)
                                 aGACAACACTCCGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGac                                                                                                                     >  2:1522581/1‑150 (MQ=255)
                                   acaacaCTCCGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGacca                                                                                                                   >  2:1687380/1‑150 (MQ=255)
                                    caacaCTCCGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGTACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGaccac                                                                                                                  <  1:750169/150‑1 (MQ=255)
                                    caacaCTCCGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGaccac                                                                                                                  <  1:83551/150‑1 (MQ=255)
                                         cTCCGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTc                                                                                                             >  1:1162952/1‑150 (MQ=255)
                                         cTCCGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTc                                                                                                             >  2:563627/1‑150 (MQ=255)
                                           ccGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCTACGAAGGTCTGCTGACCACCGTTCAc                                                                                                           >  2:266696/1‑150 (MQ=255)
                                           ccGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCAc                                                                                                           >  1:622296/1‑150 (MQ=255)
                                            cGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACg                                                                                                          <  2:1488954/150‑1 (MQ=255)
                                             gATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGc                                                                                                         <  1:412205/150‑1 (MQ=255)
                                                 ttCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGctact                                                                                                     >  2:960418/1‑150 (MQ=255)
                                                        aaGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCtt                                                                                              >  2:318421/1‑149 (MQ=255)
                                                           gcgcTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTc                                                                                           <  2:1116974/150‑1 (MQ=255)
                                                            cgcTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCa                                                                                          <  2:1808/150‑1 (MQ=255)
                                                            cgcTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCa                                                                                          <  1:780053/150‑1 (MQ=255)
                                                             gcTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAg                                                                                         <  2:1233083/150‑1 (MQ=255)
                                                             gcTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAg                                                                                         >  1:1520075/1‑150 (MQ=255)
                                                                aaCTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGaaa                                                                                      <  1:403731/150‑1 (MQ=255)
                                                                   ttCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAAcc                                                                                   <  2:801710/150‑1 (MQ=255)
                                                                   ttCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAAcc                                                                                   <  1:1422300/150‑1 (MQ=255)
                                                                     cGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCCCCGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTCCTCAGAAAACCGt                                                                                 >  2:64222/1‑150 (MQ=255)
                                                                       aCAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTg                                                                               <  1:1687380/150‑1 (MQ=255)
                                                                           atatGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCtgat                                                                                                                        <  2:1649188/105‑1 (MQ=255)
                                                                           atatGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCtgat                                                                                                                        >  1:1649188/1‑105 (MQ=255)
                                                                           atatGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATgg                                                                           <  2:347721/150‑1 (MQ=255)
                                                                            tatGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGc                                                                          <  2:622296/150‑1 (MQ=255)
                                                                             atGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCTACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACGCAGCAAACCGTTGATGGcc                                                                         >  1:509950/1‑150 (MQ=255)
                                                                                          cATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAACCCGTTGATGGCCCGTCTCACAAAGa                                                            >  2:239263/1‑150 (MQ=255)
                                                                                            tCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACt                                                          <  1:328431/150‑1 (MQ=255)
                                                                                            tCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACt                                                          <  1:563627/150‑1 (MQ=255)
                                                                                                ttCCAACGCTTCCTGTCCCACCAACTGCCTGGCTCCGCTGGCTAAATTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGgcg                                                      <  1:1701706/150‑1 (MQ=255)
                                                                                                 tCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGgcgc                                                     >  2:426312/1‑150 (MQ=255)
                                                                                                     aCGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCgcggcg                                                 <  1:878681/150‑1 (MQ=255)
                                                                                                              gCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCgcggcg                                        <  2:1520075/150‑1 (MQ=255)
                                                                                                               caccacCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGgcgc                                       >  2:659624/1‑150 (MQ=255)
                                                                                                               caccacCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGgcgc                                       >  2:707558/1‑150 (MQ=255)
                                                                                                                  caccaACTGCCTGGCTCCGCTGGCTAAATATATCAACGATAACTTCGGCATCATCGAAGGGCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTc                                    <  1:1623341/150‑1 (MQ=255)
                                                                                                                   accaACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAACACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTcc                                   >  1:738722/1‑150 (MQ=255)
                                                                                                                   accaACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTcc                                   <  2:1472039/150‑1 (MQ=255)
                                                                                                                     caACTGGCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACAGATACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGAGCTTCCCa                                 >  2:39920/1‑150 (MQ=255)
                                                                                                                     caACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCa                                 <  2:1666775/150‑1 (MQ=255)
                                                                                                                     caACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCa                                 <  1:426312/150‑1 (MQ=255)
                                                                                                                       aCTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGa                               >  1:1421281/1‑150 (MQ=255)
                                                                                                                           ccTGGCTCCGCTGGCTAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAACCCGTTGCTGGCCCGTCGCACAAAGACTGGCGCCGCGGCCGCGGCGCTTCCCAGACcaac                          >  2:927293/1‑148 (MQ=255)
                                                                                                                             tGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGTCCGCGGCGCTTCCCAGAAcctcc                         >  2:243233/1‑149 (MQ=255)
                                                                                                                             tGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAAcatca                         >  1:204659/1‑150 (MQ=255)
                                                                                                                             tGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAACGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGACGCGGCGCTTCCCAGAAcatca                         >  2:1027569/1‑150 (MQ=255)
                                                                                                                              ggCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATAGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTATCACAAAGACTGGCGCGGCGGCCGCGGCGCGTCCCAGAAcatcat                        >  1:692862/1‑150 (MQ=255)
                                                                                                                              ggCTCCGCTGGCTAAAGTTATCAACGATAACATCGGCATCATCGAAGGTATGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGTCCCGGCGCACAAAGACAGGCGCGGCGGCCGCGGCGCTTCCCAGACcagcat                        >  1:686777/1‑150 (MQ=255)
                                                                                                                               gCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGACcataata                       >  2:1261916/1‑149 (MQ=255)
                                                                                                                                  ccGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTAGCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCATCCCAGAACATCATcccc                    >  2:1699010/1‑149 (MQ=255)
                                                                                                                                    gCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTc                  >  1:1032441/1‑150 (MQ=255)
                                                                                                                                    gCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTc                  <  1:960418/150‑1 (MQ=255)
                                                                                                                                     cTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTcc                 >  2:983889/1‑150 (MQ=255)
                                                                                                                                        gCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTCCCGCTACTCAGATAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCAAGAACATCATCCCGTCctca              >  1:317629/1‑149 (MQ=255)
                                                                                                                                        gCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCctct              >  2:547183/1‑150 (MQ=255)
                                                                                                                                         ccgaaGTTATCGTTGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTa             >  2:1274941/4‑150 (MQ=255)
                                                                                                                                         ccgaaGTTATCGTTGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTa             >  1:634892/4‑150 (MQ=255)
                                                                                                                                          tAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGTGCTTCCCAGAACATCATCCCGTCCTCTAc            <  1:625437/150‑1 (MQ=255)
                                                                                                                                           aaaGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTAcc           <  2:316180/150‑1 (MQ=255)
                                                                                                                                           aaaGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTAcc           >  2:1255907/1‑150 (MQ=255)
                                                                                                                                             aGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTTCCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCgg         <  1:239263/150‑1 (MQ=255)
                                                                                                                                             aGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCAGTCCTCTACCgg         >  1:70133/1‑150 (MQ=255)
                                                                                                                                                   cAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGtgctgc   >  1:1685729/1‑150 (MQ=255)
                                                                                                                                                    aaCGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTgctgct  >  2:548991/1‑150 (MQ=255)
                                                                                                                                                    aaCGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTgctgct  <  2:1421281/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGCGAAGAAAGTGGTTATGACTGGTCCGTCTAAAGACAACACTCCGATGTTCGTTAAAGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCT  >  NC_000913/1863109‑1863406

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: