Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,125,409 T→A 5.0% E17D (GAA→GAT)  glcF ← glycolate oxidase 4Fe‑4S iron‑sulfur cluster subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,125,4090TA 5.0% 398.4 / 5.0 141E17D (GAA→GAT) glcFglycolate oxidase 4Fe‑4S iron‑sulfur cluster subunit
Reads supporting (aligned to +/- strand):  ref base T (60/71);  new base A (2/5);  total (64/77)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 143 total aligned reads displayed.

TTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAG  >  NC_000913/3125260‑3125558
                                                                                                                                                     |                                                                                                                                                     
   aGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTcc                                                                                                                                                    >  1:1533892/1‑150 (MQ=255)
   aGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGGCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTcc                                                                                                                                                    >  2:256566/1‑150 (MQ=255)
    gCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCa                                                                                                                                                   >  2:1328687/1‑150 (MQ=255)
          cTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCATTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGCATGCTGTCGGCTTCCAgcgcgc                                                                                                                                             >  2:109940/1‑150 (MQ=255)
           ttCGTTGCCTTCCAGCACCTGGTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGGTCAGCGCCCAGAAGCGGCCAGGTTGGGCAGGTTGCGGTTCAAAAGCCGCTGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAgcgcgcg                                                                                                                                            <  1:1916324/150‑1 (MQ=255)
           ttCGTTGCCGTCCAGCACCTGTTTAATCAGCTTCATGCGCCCGCCCGGGCCGTCCAGGTCTTCGCCCAGAAGCTGATTCGTTGGGCAGGTTGCGGTACAAAAACCGGATTGAACACAGGCGCGCTGGAAGCCGACAGCATCCTgcgcgcc                                                                                                                                            <  1:1272780/150‑2 (MQ=255)
             cGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCTGGAAGCCGACAGCATCCTgcgcgcctg                                                                                                                                          <  1:885508/150‑4 (MQ=255)
             cGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCTGGAAGCCGACAGCATCCTgcgcgcctg                                                                                                                                          <  1:1181194/150‑4 (MQ=255)
              gtggCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAgcgcgcgcgc                                                                                                                                         <  1:1522764/147‑1 (MQ=255)
              ggtGCCTTCCAGCACCTGTTTAAACACCTAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAgcgcgcgcgt                                                                                                                                         <  1:431151/148‑2 (MQ=255)
              gTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAgcgcgcgcgt                                                                                                                                         >  1:1037252/1‑149 (MQ=255)
              gTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAgcgcgcgcgc                                                                                                                                         >  2:732923/1‑150 (MQ=255)
              gTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCACGCGCGGCCCGTCCAGTTCATCGCCCAGCAGCAGATATGTTGGGCAGGGAGCGGTACAAAATACGCAGTGATCACAGGCGCGCAGGATGGTGTCGGCCTCCCgcgcgcgcgc                                                                                                                                         >  2:304385/1‑150 (MQ=255)
                  ccTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGCTGCTGTCGGCTTCCCGCGCGCGCGCGTTc                                                                                                                                     >  1:257398/1‑150 (MQ=255)
                  ccTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTAAAAAATCCGCAGTGCACACCGGCGCGCAGGATGCTGTCGCCTTCCATCGCGCGCGCGTTc                                                                                                                                     >  2:315399/1‑150 (MQ=255)
                   cTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCt                                                                                                                                    >  2:1880563/1‑150 (MQ=255)
                     tCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGc                                                                                                                                  >  2:1004897/1‑150 (MQ=255)
                      ccAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGcc                                                                                                                                 >  2:1150710/1‑150 (MQ=255)
                         gCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCa                                                                                                                              <  2:1583557/150‑1 (MQ=255)
                         gCACCTGGTTAATCTGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGGTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCa                                                                                                                              <  2:1860028/150‑1 (MQ=255)
                           aCCTGTTTAAACAGATAGATGCGCCCGCGCGGCCCGCCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCAtc                                                                                                                            <  1:433430/150‑1 (MQ=255)
                             cTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCAtctc                                                                                                                          >  2:41590/1‑150 (MQ=255)
                              tGTTTAATCAGATAGATGCGCCCGCGCGGCCCATCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATctct                                                                                                                         >  2:21925/1‑150 (MQ=255)
                                 ttAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCa                                                                                                                      >  2:571221/1‑150 (MQ=255)
                                 ttAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCa                                                                                                                      >  1:242662/1‑150 (MQ=255)
                                 ttAATCAGATAGAGGAGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCa                                                                                                                      <  2:1520420/150‑1 (MQ=255)
                                  tAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAg                                                                                                                     >  1:929796/1‑150 (MQ=255)
                                  tAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAg                                                                                                                     <  1:631804/150‑1 (MQ=255)
                                     tCAGCTAGATGCGCCCGCGCGGCCCGTCCAGTTCAGCGCCCCGACGCTGATAGGTAGGGCAGGTTGCGGGACAAACTCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGGTCTGCCGCATCTCTTCAGTTa                                                                                                                  <  1:297280/150‑1 (MQ=255)
                                        gatagatGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCTGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAAtt                                                                                                               <  1:230922/150‑1 (MQ=255)
                                        gatagatGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAAtt                                                                                                               >  2:1116260/1‑150 (MQ=255)
                                        gatagatGCGCCCGCGCGGCCCGGCCAGATCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAAtt                                                                                                               <  1:510802/150‑1 (MQ=255)
                                           agatGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAACTCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCGTCCCGCGCGCGCGCGGTCTGCCGCAACTCTTCAGTTATTTggg                                                                                                            >  1:1630419/1‑150 (MQ=255)
                                           agatGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTggg                                                                                                            <  1:443675/150‑1 (MQ=255)
                                              tGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAACTCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGttt                                                                                                         >  1:1287354/1‑150 (MQ=255)
                                               gcgcCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAATTCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTg                                                                                                        <  1:1311833/150‑1 (MQ=255)
                                               gcgcCCGCGCGGCCCGGCCAGTTCATCGCCCAGAAGCTGATAGGTGGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTg                                                                                                        <  1:1806459/150‑1 (MQ=255)
                                                cgcCCGCGCGGCCCGTCCAGTTCCTCGCCCCGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGc                                                                                                       <  1:1972286/150‑1 (MQ=255)
                                                cgcCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGc                                                                                                       <  1:391433/150‑1 (MQ=255)
                                                cgcCCGCGCGGCCCGGCCAGTTCATCGCCCAGAAGCTGATAGGTGGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGGTCTGCCGCATCTCTTCAGTTAATTGGGTTTGc                                                                                                       <  2:654250/150‑1 (MQ=255)
                                                 gcCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCa                                                                                                      <  1:753905/150‑1 (MQ=255)
                                                 gcCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTGGCa                                                                                                      >  1:866557/1‑150 (MQ=255)
                                                   ccGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGATGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGAGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATa                                                                                                    <  1:964334/150‑1 (MQ=255)
                                                      cgcgGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGcc                                                                                                 >  2:1919491/1‑150 (MQ=255)
                                                         ggaccGTCGAGTTCATCGCCCAGAAGGTGATAGGTTGGGCAGGTTGCGGTACAAATTCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTGGAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGc                                                                                              <  2:1951222/147‑1 (MQ=255)
                                                         ggCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGTGTTTGCATAGCCTGc                                                                                              <  2:193736/150‑1 (MQ=255)
                                                         ggCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGc                                                                                              <  1:922548/150‑1 (MQ=255)
                                                                ccAGTTCATCGCCCAGAAGCTGATAGGGTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCaa                                                                                       <  1:1306192/150‑1 (MQ=255)
                                                                     tCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCCGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCAGAGCCTGCTCCGCCAAGTTc                                                                                  >  2:818414/1‑150 (MQ=255)
                                                                     tCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTc                                                                                  >  1:860045/1‑150 (MQ=255)
                                                                      cATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTcc                                                                                 <  2:1749867/150‑1 (MQ=255)
                                                                          gCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGt                                                                             >  1:1843786/1‑150 (MQ=255)
                                                                             cAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACa                                                                          <  2:770387/150‑1 (MQ=255)
                                                                              agcaGCTGTTAGTTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACAt                                                                         <  1:793116/147‑1 (MQ=255)
                                                                               gAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATg                                                                        <  1:1328687/150‑1 (MQ=255)
                                                                               gAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATg                                                                        >  2:738314/1‑150 (MQ=255)
                                                                                    tGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGAcc                                                                   >  2:1884061/1‑150 (MQ=255)
                                                                                        agtttGGGCAGGTTGCGGTACAACATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCTTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCgggg                                                               <  2:152744/147‑1 (MQ=255)
                                                                                        aGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCgggg                                                               <  1:1004897/150‑1 (MQ=255)
                                                                                        aGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCgggg                                                               <  2:1606996/150‑1 (MQ=255)
                                                                                        aGGTTGGGCAGGTTGCGGGACAAAATCCGCAGTGAACACAGGCGCGCCGGATGCTGGCGGCTTCCAGCGCGCGCGCGTTCTGCCGCAGCTCTTCAGTTAATTGGGTTTGCATAGCCTGCGCCTCAAAGTTCCGCGTACATGCGACCgggg                                                               <  2:366195/150‑1 (MQ=255)
                                                                                         ggTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGt                                                              >  2:1857682/1‑150 (MQ=255)
                                                                                          gTTGGGCAGGGTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGGCGGCGTCCAGCGCGCGCGCGTGCTGCCGCATCGCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGtt                                                             <  1:1348147/150‑1 (MQ=255)
                                                                                           tttggCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTa                                                            <  1:317184/147‑1 (MQ=255)
                                                                                           ttGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCCAAGTTCCGCGTACATGCGACCGGGGGTa                                                            >  1:829614/1‑150 (MQ=255)
                                                                                              ggCATGTTTCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAac                                                         <  2:173731/150‑1 (MQ=255)
                                                                                                  ggTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCTTCTCTTCAGTTAATTGGGTTTGCCTAGCCTGCTCCTCAAAGTTCCGGGTACAAGCGACCGGGGTTAAACAcgc                                                     <  2:1595556/150‑1 (MQ=255)
                                                                                                         gTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAgg                                              <  2:906661/150‑1 (MQ=255)
                                                                                                           aCAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGt                                            >  2:859984/1‑150 (MQ=255)
                                                                                                           aCAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCA‑‑GCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCg                                          >  1:218803/1‑150 (MQ=255)
                                                                                                           aCAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCA‑‑GCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCg                                          >  1:211978/1‑150 (MQ=255)
                                                                                                            cAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTc                                           <  2:393091/150‑1 (MQ=255)
                                                                                                            cAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCAGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGAGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTc                                           >  2:430200/1‑150 (MQ=255)
                                                                                                             aaaaTCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCg                                          >  1:1853181/1‑150 (MQ=255)
                                                                                                             aaaaTCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCg                                          >  1:1272179/1‑150 (MQ=255)
                                                                                                              aaaTCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAACCACGCCGCAAGGGTCGa                                         >  2:1481527/1‑150 (MQ=255)
                                                                                                              aaaTCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCA‑‑GCGCGCGCGTGCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAgc                                       <  1:50596/150‑1 (MQ=255)
                                                                                                                  cctcAGTGAACCCCGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCGGCATCTCTTCAGTTAATTGGGGTTGCATAGCATGCTCCTCAAAGTTCCGCGTACATGCGTCCGGGGGTAAACACGCCGCAAGGGTCGGgctg                                     >  2:100804/4‑150 (MQ=255)
                                                                                                                  ccGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAgctg                                     <  1:178922/150‑1 (MQ=255)
                                                                                                                      aGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGt                                 >  2:464/1‑150 (MQ=255)
                                                                                                                       ttGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGAATCTCTTCAGTTAATTGGGTTTGCTTAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGtt                                <  2:874376/149‑1 (MQ=255)
                                                                                                                       gTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGtt                                <  2:889779/150‑1 (MQ=255)
                                                                                                                        ttaaCACAGGCGCGGAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCGGCATCTCTTCAGTTATTTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACGGGGGTTTATCACGACGCAAGGGTCGAGCTGCTGttt                               <  2:23284/148‑1 (MQ=255)
                                                                                                                              caGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAgctg                         >  2:223686/1‑150 (MQ=255)
                                                                                                                               aGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAgctgc                        <  1:1880563/150‑1 (MQ=255)
                                                                                                                                ggCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAgctgct                       <  1:1116260/150‑1 (MQ=255)
                                                                                                                                  cgcgcTGGAAGCCGACAGCATCCTGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTgg                     >  1:409548/1‑150 (MQ=255)
                                                                                                                                   gcgcAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGTCCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCGGTTTAAGCTGCTGGt                    >  1:1201863/1‑150 (MQ=255)
                                                                                                                                    cgcAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGTGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTg                   <  2:1204446/150‑1 (MQ=255)
                                                                                                                                     gcAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTGTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGa                  <  2:725942/150‑1 (MQ=255)
                                                                                                                                        ggATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAg               >  2:447088/1‑150 (MQ=255)
                                                                                                                                            gCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGa           <  2:1533892/150‑1 (MQ=255)
                                                                                                                                             cttgcGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTTGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGaa          <  2:251717/146‑1 (MQ=255)
                                                                                                                                             cTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGaa          >  1:1944964/1‑150 (MQ=255)
                                                                                                                                               gTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATa        <  1:570861/150‑1 (MQ=255)
                                                                                                                                                cctgCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTAGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATaa       >  1:27599/4‑150 (MQ=255)
                                                                                                                                                  ggCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAACCACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTAATAGCGGAATAAAg     >  1:675172/1‑150 (MQ=255)
                                                                                                                                                  ggCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAg     <  1:192770/150‑1 (MQ=255)
                                                                                                                                                   gCTTCCAGCGCGCGCGCGTTCAGCCGCATCTCTTCCGTTAATTGGGGTTGCAGAGCCTGCTCCTCAAAGTTGCGCTTACATGCTACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAgg    <  1:315399/150‑1 (MQ=255)
                                                                                                                                                     ttCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAg  <  2:687829/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTAAGCGTGACTTCGTTGCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAG  >  NC_000913/3125260‑3125558

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