Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,703,251 A→T 7.2% intergenic (+17/+17) add → / ← ydgJ adenosine deaminase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,703,2510AT 7.2% 209.3 / 2.4 76intergenic (+17/+17)add/ydgJadenosine deaminase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (31/39);  new base T (2/4);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCGCA  >  NC_000913/1703104‑1703400
                                                                                                                                                   |                                                                                                                                                     
ataCCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCCCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAAt                                                                                                                                                     <  1:152872/150‑1 (MQ=255)
        gCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCAccatc                                                                                                                                             <  1:1706883/150‑1 (MQ=255)
        gCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCAccatc                                                                                                                                             <  1:367816/150‑1 (MQ=255)
         ccGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCAccatcc                                                                                                                                            >  2:1689681/1‑150 (MQ=255)
              ccAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCACCATCCAtttt                                                                                                                                       <  2:563360/150‑1 (MQ=255)
              ccAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCAtttt                                                                                                                                       <  2:1592053/150‑1 (MQ=255)
              ccAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGAGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCAtttt                                                                                                                                       <  2:716877/150‑1 (MQ=255)
               cAGCTGCTGGGGTATCCCGCGAGCAAATCCGCCAGGCAAAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCACCATCCAttttg                                                                                                                                      <  1:298391/150‑2 (MQ=255)
                aGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTc                                                                                                                                     <  1:1689681/150‑1 (MQ=255)
                aGCTGCTGGGTTATCCCGCGAGAAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGATGAGGAAAAACGCGCACTGCGAGATAAAGTCGCCGCGAAGTCACAAAATGGATGGTGCAAATGCACCTTCCATTTTTc                                                                                                                                     >  2:1083893/1‑150 (MQ=255)
                 gctgctGGGTTATCCCGCTAGCAAATCCGCCAGGCACAGATTAATGGTCTGTAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCa                                                                                                                                    <  2:585099/150‑1 (MQ=255)
                 gctgctGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCa                                                                                                                                    >  1:411080/1‑150 (MQ=255)
                 gctgctGGGTTATCCCGCGAGCAAATCCCCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCACCATCCAttttgtt                                                                                                                                    <  1:808954/150‑4 (MQ=255)
                   tgctgTGTTATCCCGCGAGAAATTCCTGCATGAACAGATTAATGTTCTGGAATTGGCTTTCCTCTGCGCTGAGGAAATACGCGCACTGCGAGATAAAGTCGCCTCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATg                                                                                                                                  <  2:723022/150‑1 (MQ=255)
                   tgctgGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATg                                                                                                                                  >  1:748733/1‑150 (MQ=255)
                    gctgGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAATTGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGATGAAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGc                                                                                                                                 >  1:1217100/1‑150 (MQ=255)
                         gTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGc                                                                                                                            <  1:1352785/150‑1 (MQ=255)
                           tATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGcac                                                                                                                          >  2:1555895/1‑150 (MQ=255)
                                cgcgAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGt                                                                                                                     >  1:1193753/1‑150 (MQ=255)
                                         aTCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATg                                                                                                            <  2:1428595/150‑1 (MQ=255)
                                            cGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGttt                                                                                                         <  1:1599954/150‑1 (MQ=255)
                                              ccAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTgg                                                                                                       <  1:1703462/150‑1 (MQ=255)
                                               cAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGc                                                                                                      >  1:399412/1‑150 (MQ=255)
                                                      aGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCg                                                                                               >  1:1177162/1‑150 (MQ=255)
                                                         ttAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATg                                                                                            <  1:1231520/150‑1 (MQ=255)
                                                         ttAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATg                                                                                            <  1:1214524/150‑1 (MQ=255)
                                                         atAATGGTCTGGAAATCGGTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAATTCGCCGCGAAGTAACAAAATGGATTGTGCAAATGCACCATCCATTTTTCATGCAAGGCTCATAGTCGCGCGATGTTTGGCGGATTCGATg                                                                                            <  1:405914/149‑1 (MQ=255)
                                                          tAATGGTCTTGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCTCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGc                                                                                           <  2:1307096/150‑1 (MQ=255)
                                                           aaTGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGcc                                                                                          >  2:1374075/1‑150 (MQ=255)
                                                             tGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCa                                                                                        <  2:661621/150‑1 (MQ=255)
                                                                 cTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtc                                                                                    <  1:17657/150‑1 (MQ=255)
                                                                 cTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtc                                                                                    <  2:1193753/150‑1 (MQ=255)
                                                                   ggAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtcaa                                                                                  >  1:968123/1‑150 (MQ=255)
                                                                      aaTGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtcaatca                                                                               <  2:1573603/150‑1 (MQ=255)
                                                                      aaTGGCTTTCCTCAGCGCGGAGGAAAAACGCGCACTGCGAGAAAAAGTCGACGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtcaatca                                                                               <  1:1509614/150‑1 (MQ=255)
                                                                      aaTGGCGTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGGGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGGTTGGCGGATTCGATGCCCAGCtcaatca                                                                               <  1:1559543/150‑1 (MQ=255)
                                                                        cgtttttCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAAc                                                                             >  1:930559/5‑150 (MQ=255)
                                                                        cgtttttCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAAc                                                                             >  1:890759/5‑150 (MQ=255)
                                                                           ttttCCTCAGCGCTGAGGAAAACCGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTcc                                                                          >  1:858538/2‑150 (MQ=255)
                                                                            tttCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCa                                                                         <  1:83874/150‑1 (MQ=255)
                                                                               ccTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACCCAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAATGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTa                                                                      >  2:1678759/1‑150 (MQ=255)
                                                                               ccTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCa                                                                         <  1:1517141/147‑1 (MQ=255)
                                                                               ccTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCa                                                                         >  2:1517141/1‑147 (MQ=255)
                                                                                  cAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCt                                                                   >  2:1157046/1‑150 (MQ=255)
                                                                                  cAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCt                                                                   >  1:1272736/1‑150 (MQ=255)
                                                                                       cTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGAtt                                                              >  2:1604194/1‑150 (MQ=255)
                                                                                            gAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTg                                                         >  2:1130136/1‑150 (MQ=255)
                                                                                            gAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTg                                                         >  1:113606/1‑150 (MQ=255)
                                                                                              aaaaCGCGCACTGCGAGAAATAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGc                                                       <  1:1596929/150‑1 (MQ=255)
                                                                                              aaaaCGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAAAGCACAATCCATTTTTCAAGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGc                                                       <  1:121286/150‑1 (MQ=255)
                                                                                               aaaCGCGCACTGCGAGAAAAAGTCGCCGAGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATAACCTGGATTGCCTGGCt                                                      >  2:1461425/1‑150 (MQ=255)
                                                                                                  cgcgCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCCTGCAAGGCACAAAGTCGCGCGAGGTTTGGCGGATTCGATGACCAGCTCAATCAACTCCATTACCTGGAGTCCCTGGCTTGc                                                   >  1:295007/1‑150 (MQ=255)
                                                                                                   gcgcACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAATGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGcc                                                  <  1:865875/150‑1 (MQ=255)
                                                                                                    cgcACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCg                                                 <  2:748733/150‑1 (MQ=255)
                                                                                                      cACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGGGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGa                                               <  2:460000/150‑1 (MQ=255)
                                                                                                                  aaGTCGCCGCGAAGTAACAAAATGGATGGTGCATATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCa                                   <  1:64423/150‑1 (MQ=255)
                                                                                                                    gTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCAcc                                 >  2:471564/1‑150 (MQ=255)
                                                                                                                    gTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCAcc                                 <  1:1071527/150‑1 (MQ=255)
                                                                                                                        ccgcGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATcc                             >  2:800497/1‑149 (MQ=255)
                                                                                                                        ccgcGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCg                             >  1:639181/1‑150 (MQ=255)
                                                                                                                            gAAGTAACATAATGGAAGGAGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGAAGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCCCCATCGCCAt                         <  1:1083893/150‑1 (MQ=255)
                                                                                                                            gAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCAGTACCTGGATTGCCTGGCTTGCCGGCACCGGATTTTCACCATCGCCAt                         >  1:825554/1‑150 (MQ=255)
                                                                                                                            gAAGTAACAAAATGGATGGGGCAAATGCACCATCCATTTTTCATGCATGCCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCAGCATCGCCAt                         <  2:1217100/150‑1 (MQ=255)
                                                                                                                                 aaCAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAAc                    <  2:1334274/150‑1 (MQ=255)
                                                                                                                                   cAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGc                  >  1:1291007/1‑150 (MQ=255)
                                                                                                                                    aaaaTGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCAcc                                 <  1:974953/134‑1 (MQ=255)
                                                                                                                                    aaaaTGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCAcc                                 >  2:974953/1‑134 (MQ=255)
                                                                                                                                    aaaaTGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCa                 <  1:1555895/150‑1 (MQ=255)
                                                                                                                                     aaaTGGATGGTGCATTTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCAt                >  2:275081/1‑150 (MQ=255)
                                                                                                                                     aaaTGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCGGGATTGCCTGGCTTTCCGGAACCGGATTTTCACCATCGCCATTTAACGCAt                >  1:920938/1‑150 (MQ=255)
                                                                                                                                        tggatggTGCATTTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCAc             >  1:703727/1‑150 (MQ=255)
                                                                                                                                         ggatggTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACg            <  2:719694/150‑1 (MQ=255)
                                                                                                                                         ggatggTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACg            <  2:756606/150‑1 (MQ=255)
                                                                                                                                           atggTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGaa          <  2:1177162/150‑1 (MQ=255)
                                                                                                                                             ggTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATa        >  1:669298/1‑150 (MQ=255)
                                                                                                                                              gTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAg       <  2:858538/150‑1 (MQ=255)
                                                                                                                                                   aaTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCGCa  <  2:645307/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCGCA  >  NC_000913/1703104‑1703400

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: