Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,703,251 A→T 9.3% 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 9.3% 226.6 / 8.2 71intergenic (+17/+17)add/ydgJadenosine deaminase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (33/31);  new base T (3/4);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ACCAGCACCGTAGCAGAGCTGGCTGCACATCCGCTGAAAACGTTCCTTGAGCATGGCATTCGTGCCAGCATTAACACTGACGATCCCGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCGCATAGTAAGCCGGATAATTCCCAGGCACCGTCAACAGCGTTTCTTCGACACGTTCCTCACCTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACG  >  NC_000913/1702986‑1703519
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                            
aCCAGCACCGTAGCAGAGCTGGCTGCACATCCGCTGAAAACGTTCCTTGAGCATGGCATTCGTGCCAGCATTAACACTGACGATCCCGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGCGGAAAAACGCCCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCAccat                                                                                                                                                                                                                                                                     >  2:117119/1‑275 (MQ=255)
     cACCGTAGCAGAGCTGGCTGCACATCCGCTGAAAACGTTCCTTGAGCATGGCATTCGTTCCAGCATTACCACTGACGATCCCGGCGTACAGGGAGTGGATATCATTCAGGAATATACCGTTGCCGCTACAGCTGATGGTTTATCACGCGTGCAAATACGCCAGGCACAGGTTAATGGCTTGGAAATGGCTTTCCTCAGCGCTGCGGAAAAACGCTCACATCGCCACAAATTCGCCACGAATTAACACAAGCGAGGGTGCaacctgaccactcat                                                                                                                                                                                                                                                                 >  4:414385/1‑261 (MQ=255)
          tAGCAGAGCTGGCTGCACATCCGCTGAAAACGTTCCTTGAGCATGGCATTCGTGCCAGCATTAACACTGACGATCCCGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCa                                                                                                                                                                                                                                                           >  2:204040/1‑275 (MQ=255)
               aagCTGGCTGCACATCCGCTGAAAACGTTCCTTGAGCATGGCATTCGTGCCAGCATTAACACTGACGATCCCGGCGTACAGGGAGTGGATATCATTCACGAATATCCCGTTGCCGCGCCAGCTGCTGAGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCa                                                                                                                                                                                                                                                       <  1:93202/273‑1 (MQ=255)
                  ctggctgCACATCCGCTGAAAACGTTCCTTGAGCATGGCATTCGTGCCAGCATTAACACTGACGATCCCGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAATGc                                                                                                                                                                                                                                                   >  2:152535/1‑275 (MQ=255)
                                                          gtcggtccAGCATTTACGCTGACGATCCCGGCGGAAAGGGAGTGTGTATCATTAACGAATAGACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATAGGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtc                                                                                                                                                                                                           <  4:117118/269‑1 (MQ=255)
                                                                ccAGAATGAACACTGACCATCCGAGCGTACAGGGAGTTGATATCTTTCATGACGATAGCGTGGCCGAGGCCGTCGCTGGGTGTGCCCGCGGGCAAATCCGCCAGGCACAGATTAATTGTCTGGAAAATGCTGTCGTAAGGGCTGAGGAAAAACGTGCACAGCGAGAAAAAGTCGCCGCGGAGTAACAAAATGGATTGTTCAAATGCACCAGCCATTTTTCATGCAAGGCACAAAGTCTCGCGATGTTTGGCGGATTCGATGCCCAGCtcaatcaa                                                                                                                                                                                                     <  4:492147/275‑1 (MQ=255)
                                                                  aGCATTAACACTGACGATCCCGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACt                                                                                                                                                                                                   >  2:6269/1‑275 (MQ=255)
                                                                               aCGATCCCGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAAGGGTCTGGAAATGGCTTTCCTCAGTGCTGAGGAAAAACGAGCACTGCGAGCAAAAGTAGCCGAGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTATTGCACGGAAACAAGCCGCGCGATGCATGGCGGACTCGAGCCCCACCCCCATCCACCCCAATACCTGGAt                                                                                                                                                                                      >  4:273211/1‑275 (MQ=255)
                                                                                   tCCCGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGc                                                                                                                                                                                                                  <  2:311527/243‑1 (MQ=255)
                                                                                           tACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGc                                                                                                                                                                          <  2:414390/275‑1 (MQ=255)
                                                                                                         atCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGAtttt                                                                                                                                                            >  1:147796/1‑275 (MQ=255)
                                                                                                           cATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACACAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTCACTGGATTGCCTGGCTTGCCGGAACCGGATTTTc                                                                                                                                                           >  2:172126/1‑274 (MQ=255)
                                                                                                            aTTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCACCTCCATTACCTGGATTTCCTGGCTTCCCGGAACCGGAattttc                                                                                                                                                          >  4:374003/1‑272 (MQ=255)
                                                                                                                aCGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCAc                                                                                                                                                                                                                                                                        <  1:103560/160‑1 (MQ=255)
                                                                                                                   aaTATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCg                                                                                                                                                    <  2:374008/273‑1 (MQ=255)
                                                                                                                       taCCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAAt                                                                                                                                                                                                                                                                            <  7:337601/149‑1 (MQ=255)
                                                                                                                       taCCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAAt                                                                                                                                                                                                                                                                            >  8:35214/1‑149 (MQ=255)
                                                                                                                        aCCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTACTGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGANGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGAGGCCCAGCCAAATCAANCTCATTACCTGGATTGCCTGGCTTGNCGTTACCCGCCTTTCACCATCCCCATTTa                                                                                                                                             >  3:97908/1‑275 (MQ=255)
                                                                                                                         ccGTTGCGGCGCCAGTGGCTGGGTTACCCCCCGAGCAAATGCACCAGGCACAGATGAATGGTCTGGTAATGGCTTGACTCAGCACTGAGGGAAAACGCGCACTGCGAGACAAAGTAGCCGCGAAGGAACAAAATGGATGGTGCAAAGGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCAATGTGTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTaa                                                                                                                                            <  4:204037/275‑1 (MQ=255)
                                                                                                                          cGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCa                                                                                                                                                                                                                                                                         <  7:62894/149‑1 (MQ=255)
                                                                                                                            ttGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGAAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACg                                                                                                                                          >  2:341020/1‑274 (MQ=255)
                                                                                                                              gCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccat                                                                                                                                                                                                                                                                     <  6:363250/149‑1 (MQ=255)
                                                                                                                               ccGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATc                                                                                                                                      >  2:79042/1‑275 (MQ=255)
                                                                                                                                      aGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCAttttt                                                                                                                                                                                                                                                             >  5:74022/1‑149 (MQ=255)
                                                                                                                                           ctgGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGc                                                                                                                                                                                                                                                        >  7:132490/1‑149 (MQ=255)
                                                                                                                                            tgGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCa                                                                                                                                                                                                                                                       >  6:333129/1‑149 (MQ=255)
                                                                                                                                                ttATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGc                                                                                                                                                                                                                                                   <  6:290183/149‑1 (MQ=255)
                                                                                                                                                     ccGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTc                                                                                                                                                                                                                                           >  2:61030/1‑152 (MQ=255)
                                                                                                                                                        cgAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTg                                                                                                                                                                                                                               >  2:220339/1‑161 (MQ=255)
                                                                                                                                                         gAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTcgcgc                                                                                                                                                                                                                                       <  4:61029/152‑1 (MQ=255)
                                                                                                                                                          aGCAAAACCGCCAGGCAAAAGTTAAGGGTGTGGAAATGGCTGTCCCCAGCTCTGAGGAAAAAAGCGGCATGCGAGAAAAAGTCGCCTCGAAGTATCAAAATGGATGGTGCAAAGTCACCATCCATTTTTCATGCAAGGAACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCGCATAGTAAGCCGGAt                                                                                                            <  4:172124/274‑1 (MQ=255)
                                                                                                                                                            cAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTcgcgc                                                                                                                                                                                                                                       <  8:132490/149‑1 (MQ=255)
                                                                                                                                                                tCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATg                                                                                                                                                                                                                                   >  5:102931/1‑149 (MQ=255)
                                                                                                                                                                    ccAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTAGCGCGATGTTTg                                                                                                                                                                                                                               <  7:35214/149‑1 (MQ=255)
                                                                                                                                                                     cAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTgg                                                                                                                                                                                                                              <  7:340165/149‑1 (MQ=255)
                                                                                                                                                                      aGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGc                                                                                                                                                                                                                             >  5:332525/1‑149 (MQ=255)
                                                                                                                                                                          acaGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTc                                                                                                                                                                                                                                                            <  8:164585/114‑1 (MQ=255)
                                                                                                                                                                          acaGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTc                                                                                                                                                                                                                                                            >  7:164585/1‑114 (MQ=255)
                                                                                                                                                                                tAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATg                                                                                                                                                                                                                   <  6:332525/149‑1 (MQ=255)
                                                                                                                                                                                      tCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGGGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGc                                                                                                                                                                                                             <  5:333129/149‑1 (MQ=255)
                                                                                                                                                                                        tGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAGTGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCa                                                                                                                                     >  2:105873/1‑219 (MQ=255)
                                                                                                                                                                                                  tttCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTcc                                                                                                                                                                                                 >  7:364185/1‑149 (MQ=255)
                                                                                                                                                                                                    tCCTAAGCCCTGAGGAAAAACGCGCACTGCGCGACAACGTCGCCGCGAAGTAACAAAATGGATGGTGCATATTCCCCATCCATTTTTCATGCAAGGCACAAAGTCGCTCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCGCATAGTAATCCGGATAATTCCCAGGCACCGTCAACAGCGTTTCTTCGACACGTTCCt                                                                  <  1:55698/274‑1 (MQ=255)
                                                                                                                                                                                                     ccTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAAt                                                                                                                                                                                                                                                                            <  7:5624/71‑1 (MQ=255)
                                                                                                                                                                                                     ccTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAAt                                                                                                                                                                                                                                                                            >  8:5624/1‑71 (MQ=255)
                                                                                                                                                                                                          gcgcTGAGGAAAAACGCGCAGTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCCTCACGAATAGCCGCATAGTAAGCCGGATAATTCCCAGGCACCGTCCACCGCGTTTCGTCGAACCGGTCCTCAACTTc                                                           >  2:96795/1‑275 (MQ=255)
                                                                                                                                                                                                                aGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCGCATAGTAAGCCGGATAATTCCCAGGCACCGTCAACAGCGTTTCTTCGACACGTTCCTCACCTTCCAcgc                                                      >  2:24242/1‑274 (MQ=255)
                                                                                                                                                                                                                 agAAAAACGCGCAATTCGTGAAAAAGTCACCGCGAAGTAACACAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCGCATAGTAAGCCGGATAATTCCCAGGCACCGTCAACAGCGTTTCTTCGACACGTTCCTCACCTTCCACGCgg                                                    <  2:455065/274‑1 (MQ=255)
                                                                                                                                                                                                                  gAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCt                                                                                                                                                                                 >  6:396984/1‑149 (MQ=255)
                                                                                                                                                                                                                    aaaaCGCGAACTGCGAGAAAAAGTCGCCGCGGAGTAACAAAAGTGATGGTGCAAATGCGCCGTCCATTTTTCATGCAAGGCACAAAGTGGCGCGATGTTGGGCGGATTCGATGCCCAGCTCAATCAACTCCATTAACTGGATTGCCTGTCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCGCATAGTAAGCCGGATAATTCCCAGGCACCGTCAACAGCGTTTCTTCGAAACGTTCCTCACCTTCCACGCGGGt                                                  <  4:24242/274‑1 (MQ=255)
                                                                                                                                                                                                                         gcgcACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGTCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCGCATAGTAAGCCGGATAATTCCCAGGCACCGTCAACAGCGTTTCTTCGACACGTTCCTCACCTTCCACGCGGGTAAGTa                                             <  2:126089/274‑1 (MQ=255)
                                                                                                                                                                                                                         gcgcACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCGCATAGTAAGCCGGATAATTCCCAGGCACCGTCAACAGCGTTTCTTCGACACGTTCCTCACCTTCCACGCGGGTAAGTa                                             <  2:273215/274‑1 (MQ=255)
                                                                                                                                                                                                                          cgcACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGcc                                                                                                                                                                         <  6:74022/149‑1 (MQ=255)
                                                                                                                                                                                                                          cgcACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCGCATAGTAAGCCGGATAATTCCCAGGCACCGTCAACAGCGTTTCTTCGACACGTTCCTCACCTTCCACGCGGGTAAGTa                                             <  1:97908/273‑1 (MQ=255)
                                                                                                                                                                                                                                    gaAAAAGTAGACAGGAGGGAACAAAAAGGNNNGGTAAAATCAACCATCAATTTTTGAGGCAAGGGACAAAGGCGCGCGTTGTTTGGCGGATTCGATGCCCAGCTCAATCAAACCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCGTCACGAATAGCCGCATAGTAAGCCGGATAATTCCCAGGCACCGTCAACAGCGTTTCTTCGACACGTTCCTCACCTTCCACGCGGGTAAGTACGCCATCACGCa                                 <  3:147796/275‑1 (MQ=255)
                                                                                                                                                                                                                                         ggTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACGAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCa                                                                                                                                                          <  7:458999/148‑1 (MQ=255)
                                                                                                                                                                                                                                                      tAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATc                                                                                                                                                     >  7:88177/1‑141 (MQ=255)
                                                                                                                                                                                                                                                      tAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATc                                                                                                                                                     <  8:88177/141‑1 (MQ=255)
                                                                                                                                                                                                                                                      tAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTa                                                                                                                                             >  5:278622/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      atgaaaaaTGGATGGTGCATTTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGccc                                                                                                                                                                                                                >  8:197621/5‑82 (MQ=25)
                                                                                                                                                                                                                                                      atgaaaaaTGGATGGTGCATTTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGccc                                                                                                                                                                                                                <  7:197621/78‑1 (MQ=25)
                                                                                                                                                                                                                                                          aaaaTGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGc                                                                                                                                         <  5:396984/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            aaTGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCAt                                                                                                                                       >  5:273719/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               ggatggTGCATTTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCAc                                                                                                                                    >  5:310407/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                gatggTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCGCATAGTAAGCCGGATAATTCCCAGGCACCGTCAACAGCGTTTCTTCGACACGTTCCTCACCTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAg     >  2:7707/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                gatggTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATCGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCGCATAGTAAGCCGGATAATTCCCAGGCACCGTCAACAGCGTTTCTTCGACACGTTCCTCACCTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAg     >  2:213961/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                gatggTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGACGTTTCGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGTTTTTCACCATCGCCATTTAACGCATCACGAATAGCCGCATAg                                                                                                                      >  2:154849/1‑162 (MQ=255)
                                                                                                                                                                                                                                                                 atgtTGCATTTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTTCCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGa                                                                                                                                  >  8:325234/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  tggTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGaa                                                                                                                                 <  6:278622/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   ggTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCCTTTAACGCATCACGAATAGCCGCCTAGTAAGCCGGATAATTCCCAGGCACCGTCAACAGCGTTTCTTCGACACGTTCCTCACCTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGTTAGACg  >  2:422319/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                    ttGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCGCATAGTAAGCCGGATAATTCCCAGGCACCGTCAACAGCGTTTCTTCGACACGTTCCt                                                                  <  2:198374/209‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aaTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCGc                                                                                                                          >  6:318170/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                            
ACCAGCACCGTAGCAGAGCTGGCTGCACATCCGCTGAAAACGTTCCTTGAGCATGGCATTCGTGCCAGCATTAACACTGACGATCCCGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCGCATAGTAAGCCGGATAATTCCCAGGCACCGTCAACAGCGTTTCTTCGACACGTTCCTCACCTTCCACGCGGGTAAGTACGCCATCACGCATATCGTAGCCCCAGTCTTCCTGCGGTAGACG  >  NC_000913/1702986‑1703519

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: