Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,407,898 A→C 5.5% Q113P (CAG→CCG) ‡ alaA → glutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,407,8980AC 5.5% 460.2 / 9.6 149Q113P (CAG→CCG) ‡alaAglutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase
Reads supporting (aligned to +/- strand):  ref base A (71/69);  new base C (7/2);  total (78/71)
Fisher's exact test for biased strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 153 total aligned reads displayed.

TACCGCTCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCG  >  NC_000913/2407749‑2408047
                                                                                                                                                     |                                                                                                                                                      
tACCGCTCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGTTGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGACCTTATCGTTCAGGCAATGCa                                                                                                                                                        >  2:1319596/1‑150 (MQ=255)
tACCGCTCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGCTATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCa                                                                                                                                                        >  1:1049030/1‑150 (MQ=255)
  ccGCTCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAgg                                                                                                                                                      <  2:569783/150‑1 (MQ=255)
  ccGCTCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGCCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGAAgg                                                                                                                                                      >  2:1303521/1‑150 (MQ=255)
  ccGCTCAAGGGTATTGAGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAgg                                                                                                                                                      >  1:1843601/1‑150 (MQ=255)
     cTCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCAt                                                                                                                                                   >  2:1102476/1‑150 (MQ=255)
     cTCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGACAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCAt                                                                                                                                                   >  1:1376488/1‑150 (MQ=255)
        aaGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGc                                                                                                                                               <  1:1773720/150‑1 (MQ=255)
        aaGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGc                                                                                                                                               >  1:736521/1‑150 (MQ=255)
          gggTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTg                                                                                                                                             <  1:523894/150‑1 (MQ=255)
          cggaaTTGCGAATCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACAAGTAAATTGTCATGCGAGATGTTACCGTTTAAGATATTTACATCGTAAATGGTGTATCGGAGCTTATGTTTCAGGCAATGCAGGCATTGCTg                                                                                                                                             <  2:248047/146‑1 (MQ=255)
           ggTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCCTGCATTGcctg                                                                                                                                            <  1:1675415/150‑4 (MQ=255)
             tATTGCGATTCCAACGGTCTTTACTCAGCGCGTAAAGCCATCATGAAGCACTACCTGGCTCGTGTCATGCGTGATGTTCCCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAAc                                                                                                                                          <  2:1296905/150‑1 (MQ=255)
             tATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGCAATGCCTGCATTGCCTGAAc                                                                                                                                          <  1:1690092/150‑1 (MQ=255)
               ttGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAg                                                                                                                                        >  2:1215593/1‑150 (MQ=255)
                tGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCTTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTAAGTCAATGCAGTCATTTCTGAAAAGc                                                                                                                                       >  2:1249270/1‑150 (MQ=255)
                       ccAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACg                                                                                                                                >  1:1707956/1‑150 (MQ=255)
                             gTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGt                                                                                                                          <  1:702701/150‑1 (MQ=255)
                              tCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGtt                                                                                                                         >  1:2009417/1‑150 (MQ=255)
                                tttACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCATGCATTGCTGAACAGCGGGGACGAAATGTTgg                                                                                                                       >  1:137450/1‑150 (MQ=255)
                                tttACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTgg                                                                                                                       <  1:419478/150‑1 (MQ=255)
                                  tACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCGCGGGGCATGCGGGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGtt                                                                                                                     <  2:1568206/150‑1 (MQ=255)
                                    cTCCGCGCGTAACGCCATCATGCTGCACTACCAGGCTCGTGACATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGTAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTcc                                                                                                                   <  2:951746/150‑1 (MQ=255)
                                    cTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATATACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCAATGCTGAACAGCGGGGACGAAATGTTGGTTcc                                                                                                                   <  1:497526/150‑1 (MQ=255)
                                          gcgTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCAcc                                                                                                             <  1:597414/150‑1 (MQ=255)
                                          gcgTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCGTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCAcc                                                                                                             <  1:332323/150‑1 (MQ=255)
                                             tAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGa                                                                                                          >  2:634766/1‑150 (MQ=255)
                                               aaGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGAtt                                                                                                        >  2:763457/1‑150 (MQ=255)
                                                aGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTa                                                                                                       >  1:1523569/1‑150 (MQ=255)
                                                aGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTa                                                                                                       >  2:699365/1‑150 (MQ=255)
                                                 gCCATCATGCAGCACTACCAGGCTCGTGGCATGCGGGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTAc                                                                                                      <  1:256400/150‑1 (MQ=255)
                                                        tGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCActct                                                                                               <  2:1579412/150‑1 (MQ=255)
                                                        tGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCActct                                                                                               >  2:116658/1‑150 (MQ=255)
                                                          cagcaCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCTGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTgg                                                                                             >  2:37546/1‑150 (MQ=255)
                                                          cagcaCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTgg                                                                                             >  2:1952131/1‑150 (MQ=255)
                                                              aCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACcg                                                                                         <  1:1215593/150‑1 (MQ=255)
                                                                 aCCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGTTTCCTGCAACAGATTACACACTCTGGACCgcgg                                                                                      >  1:1952273/1‑150 (MQ=255)
                                                                   cAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCgcggcg                                                                                    >  1:562858/1‑150 (MQ=255)
                                                                   cAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATCGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCgcggcg                                                                                    <  2:1045293/150‑1 (MQ=255)
                                                                    aGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGcggcgg                                                                                   <  2:1843601/150‑1 (MQ=255)
                                                                     ggCTCGTGGCATGCGTGATGTTACCGGGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAATTGTTGGTTCCTGCACCAGATTACCCACTCTGGTCCGCGGCGGt                                                                                  <  1:173779/150‑1 (MQ=255)
                                                                      gCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGtt                                                                                 <  2:689617/150‑1 (MQ=255)
                                                                            ggCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCtt                                                                           >  2:1691488/1‑150 (MQ=255)
                                                                            ggCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTATCGCtt                                                                           >  2:889347/1‑150 (MQ=255)
                                                                             gCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCttt                                                                          <  1:2013527/150‑1 (MQ=255)
                                                                               ctGCGTGATGATACCGTGGAAGATAATTACCTCGGCAATGGTGTATCGGAGCTTATCGTTCCGGCCATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTcc                                                                        <  1:365121/149‑1 (MQ=255)
                                                                               aTGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTcc                                                                        <  1:1409545/150‑1 (MQ=255)
                                                                               aTGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTcc                                                                        >  1:1647767/1‑150 (MQ=255)
                                                                                   gTGATTTTACCGTGGAATATATTTACATCGGCAATGGTGTATCGTAGCTTATCGTTCAGGCAATGCAGGCATTGCGGAACAGCGGGGACGAAATTTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCg                                                                    <  1:1870240/150‑1 (MQ=255)
                                                                                   gTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCg                                                                    <  2:1738409/150‑1 (MQ=255)
                                                                                    tGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCgg                                                                   <  2:173558/150‑1 (MQ=255)
                                                                                         ttACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCAACAGATTACCCACTCTGGACCGCGGCGGGTTCGCTTTCCAGCGGTAAAg                                                              >  2:363574/1‑150 (MQ=255)
                                                                                          tACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGCACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTGTCGAGCGGTAAAGc                                                             >  1:297567/1‑150 (MQ=255)
                                                                                              gTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCGGGTTACCCAATCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTg                                                         >  2:296994/1‑150 (MQ=255)
                                                                                              gTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTg                                                         >  1:1973331/1‑150 (MQ=255)
                                                                                               tGGAAGATCTTTACATCGGCAATGGGGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGCTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGGTTCGCTTTCCAGCGGTAAAGCGGTGc                                                        <  1:1741369/150‑1 (MQ=255)
                                                                                               tGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCCATGCAGGCAGTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGc                                                        <  1:1930730/150‑1 (MQ=255)
                                                                                                      tatTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCt                                                 >  1:112018/1‑150 (MQ=255)
                                                                                                      tatTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCt                                                 <  1:1938809/150‑1 (MQ=255)
                                                                                                         ttACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTg                                              >  2:527528/1‑150 (MQ=255)
                                                                                                          tACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGc                                             <  2:1121290/150‑1 (MQ=255)
                                                                                                          tACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGc                                             <  1:1944106/150‑1 (MQ=255)
                                                                                                          tACATCGGCAATGGTGCATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGc                                             <  1:881030/150‑1 (MQ=255)
                                                                                                           aCATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCg                                            <  2:1977850/150‑1 (MQ=255)
                                                                                                           aCATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCg                                            <  2:1271222/150‑1 (MQ=255)
                                                                                                           aCATCGACAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCg                                            <  2:1376488/150‑1 (MQ=255)
                                                                                                            cATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGa                                           <  1:1149914/150‑1 (MQ=255)
                                                                                                             aTCGGCAATGGTGTATCGGGGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTGCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGAt                                          <  2:1835522/150‑1 (MQ=255)
                                                                                                               cGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACACACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCAGTATCTTTGCGATGa                                        >  2:335895/1‑150 (MQ=255)
                                                                                                                ggCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGaa                                       >  1:1449644/1‑150 (MQ=255)
                                                                                                                  cAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATc                                     >  2:1534295/1‑150 (MQ=255)
                                                                                                                   aaTGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTTGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATcc                                    >  2:1311691/1‑150 (MQ=255)
                                                                                                                       gtgtATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTg                                <  2:1990507/150‑1 (MQ=255)
                                                                                                                       gtgtATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTg                                <  1:116658/150‑1 (MQ=255)
                                                                                                                         gtATCGGAGCTTATCGTTCAGGCAATGCAGTCATTGCTGAACAGCGTGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGGTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGAc                              <  1:329822/150‑1 (MQ=255)
                                                                                                                         gtATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGAc                              <  2:1950057/150‑1 (MQ=255)
                                                                                                                           aTCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTg                            >  1:1771386/1‑150 (MQ=255)
                                                                                                                              ggAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGtt                         >  2:1917617/1‑150 (MQ=255)
                                                                                                                                 gCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTccc                      <  2:1336523/150‑1 (MQ=255)
                                                                                                                                  cccgctgTTCAGCAATGCCTGCATTGCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTcccc                     >  2:1616001/7‑149 (MQ=255)
                                                                                                                                   ttATCGTTCAGGCAATGCAGTCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCAGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCgg                    <  2:1024418/150‑1 (MQ=255)
                                                                                                                                   ttATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCgg                    <  2:1519945/150‑1 (MQ=255)
                                                                                                                                   ccgctgTTCAGCAATGCCTGCATTGCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCgg                    >  2:1689553/6‑150 (MQ=255)
                                                                                                                                     gctgTTCAGCAATGCCTGCATTGCCGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGAcc                 >  1:1119125/4‑150 (MQ=255)
                                                                                                                                     aTCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAATTGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCggcc                  >  1:411475/1‑148 (MQ=255)
                                                                                                                                       cGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCt                >  1:1491662/1‑150 (MQ=255)
                                                                                                                                        gTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTc               <  2:600387/150‑1 (MQ=255)
                                                                                                                                        gTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTc               >  1:45141/1‑150 (MQ=255)
                                                                                                                                           cAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCgat            <  1:1456262/150‑1 (MQ=255)
                                                                                                                                           cAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCgat            <  1:190267/150‑1 (MQ=255)
                                                                                                                                             ggCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCgatga          <  2:1588604/150‑1 (MQ=255)
                                                                                                                                             ggCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCgatga          <  2:830931/150‑1 (MQ=255)
                                                                                                                                               cAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGata        >  1:1214965/1‑150 (MQ=255)
                                                                                                                                               cAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGata        >  1:296788/1‑150 (MQ=255)
                                                                                                                                               cAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTACTGCACCAGATTACACACTCTGGACCGCGGCGGTTTCGCTTTCCACCGGTAAAGCGGTGCATTATCAATGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGata        >  2:160852/1‑150 (MQ=255)
                                                                                                                                                   gCAGGCATTGCTGCACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCAGTAGCTTTGCGAGGAATCCTCTGACTGGTTCCCGCCCCGCGAGGAGATTcg    >  1:1008827/1‑150 (MQ=255)
                                                                                                                                                   gCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTcg    >  1:838639/1‑150 (MQ=255)
                                                                                                                                                   gCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGAGCATTATCTTTGCGATGAATCCTCTGACTGGGTCCCGGACCACGATGATATTcg    >  2:747581/1‑150 (MQ=255)
                                                                                                                                                     aGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTcgcg  >  1:859269/1‑150 (MQ=255)
                                                                                                                                                     aGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTcgcg  >  2:1926520/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TACCGCTCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCG  >  NC_000913/2407749‑2408047

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: