Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 317,566 G→A 7.5% pseudogene (128/129 nt) ykgQ → pseudogene, putative dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913317,5660GA 7.5% 258.1 / 10.1 80pseudogene (128/129 nt)ykgQpseudogene, putative dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (38/36);  new base A (3/3);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-02

GTTACGTTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTGCGTTGTTTAAGAATTTCCCGCGCCGAATCTGCCATTATCATGACTGCACCTGCCAGCATCAG  >  NC_000913/317335‑317827
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                     
gTTACGTTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTaa                                                                                                                                                                                                                            <  4:502820/275‑1 (MQ=255)
 ttACGTTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAAc                                                                                                                                                                                                                           >  4:147237/1‑275 (MQ=255)
   aCGTTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTTGTTGTCTTGCTTTAAAAATACATTATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAACTTTTaaacg                                                                                                                                                                                                                         >  4:412595/1‑272 (MQ=255)
     gTTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACa                                                                                                                                                                                                                          <  2:147238/272‑1 (MQ=255)
                        gTTCATGTCTAATGGAGGGACTGACAGTTNNATTTTCCGTTTAGTCTGCTTCAGAGTATTGNAGATGCCCATCAAAAATCGATTTTTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTANGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCAct                                                                                                                                                                                                    <  3:57722/275‑1 (MQ=255)
                                   aTGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACATAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCtt                                                                                                                                                                                                        >  1:57722/1‑260 (MQ=255)
                                    tGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATAtctc                                                                                                                                                                                        <  4:363672/275‑1 (MQ=255)
                                                                                  ttGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACCTTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATTTATTTTACCTAAGACGACAAAAACCTTTTAACCGGGAGTGATATTGACCTTCCCTCTGTCATATCTCCGGTAATATCGCGTCAGGCTTCATCATACGGAAGGGGGCCTGACGc                                                                                                                                          >  2:89582/1‑275 (MQ=255)
                                                                                      aGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACAtgattgat                                                                                                                                                                                                                                                                     <  8:852634/148‑1 (MQ=255)
                                                                                        aTGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATCTGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGaa                                                                                                                                     >  2:186276/1‑274 (MQ=255)
                                                                                        aTGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGaaa                                                                                                                                    >  4:286480/1‑275 (MQ=255)
                                                                                          gCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGNTTGATATTAATCATGTATTTTACCTAAGACGACAAAAACATTTTAACAGGGAGTGATATTGATCTTCACTCTGTNATATCTCCGGTAATATGTCGTAAGGCTTCANNATAAGAAAGGGGGCCTGACGCCTGaaa                                                                                                                                    >  3:68392/1‑273 (MQ=255)
                                                                                              aTCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGt                                                                                                                                <  1:68392/273‑1 (MQ=255)
                                                                                                        gATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTaa                                                                                                                                                                                                                                                                >  7:500854/1‑135 (MQ=255)
                                                                                                        gATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTaa                                                                                                                                                                                                                                                                <  8:500854/135‑1 (MQ=255)
                                                                                                         aTTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGCTCTTCACTCTGTCATATCTCCGGTAATATGCCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCGGAAAAAGGGAacatcaaagag                                                                                                                   >  4:526665/1‑270 (MQ=255)
                                                                                                           ttCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGatat                                                                                                                                                                                                                                                                   >  8:412511/1‑129 (MQ=255)
                                                                                                           ttCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGatat                                                                                                                                                                                                                                                                   <  7:412511/129‑1 (MQ=255)
                                                                                                                  tGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTaa                                                                                                                                                                                                                                                                <  6:415001/125‑1 (MQ=255)
                                                                                                                  tGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTaa                                                                                                                                                                                                                                                                >  5:415001/1‑125 (MQ=255)
                                                                                                                     cAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCa                                                                                                                                                                                                                                                             >  6:278195/1‑125 (MQ=255)
                                                                                                                     cAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCa                                                                                                                                                                                                                                                             <  5:278195/125‑1 (MQ=255)
                                                                                                                        aTAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGt                                                                                                                                                                                                                    >  4:490085/1‑163 (MQ=255)
                                                                                                                         tAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTa                                                                                                                                                                                                                                                <  5:850848/134‑1 (MQ=255)
                                                                                                                         tAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTa                                                                                                                                                                                                                                                >  6:850848/1‑134 (MQ=255)
                                                                                                                           aaGAAATGGCACATATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAttt                                                                                                                                                                                                                                                      <  8:741097/126‑1 (MQ=255)
                                                                                                                           aaGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTCGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAttt                                                                                                                                                                                                                                                      >  7:741097/1‑126 (MQ=255)
                                                                                                                            aGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGata                                                                                                                                                                                                                <  2:490086/163‑1 (MQ=255)
                                                                                                                             gAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCtttt                                                                                                                                                                                                                              >  7:900689/1‑148 (MQ=255)
                                                                                                                                    cacaAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAggg                                                                                                                                                                                                                       <  8:900689/148‑1 (MQ=255)
                                                                                                                                      caAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTa                                                                                                                                                                                                                                                <  6:922204/121‑1 (MQ=255)
                                                                                                                                      caAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTa                                                                                                                                                                                                                                                >  5:922204/1‑121 (MQ=255)
                                                                                                                                      caAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGt                                                                                                                                                                                                                                                          >  6:578140/1‑111 (MQ=255)
                                                                                                                                      caAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGt                                                                                                                                                                                                                                                          <  5:578140/111‑1 (MQ=255)
                                                                                                                                               tCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCa                                                                                                                                                                                                                                                             <  6:153618/99‑1 (MQ=255)
                                                                                                                                               tCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCa                                                                                                                                                                                                                                                             >  5:153618/1‑99 (MQ=255)
                                                                                                                                                  cTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgacg                                                                                                                                                                                                                                           >  6:748486/1‑114 (MQ=255)
                                                                                                                                                  cTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgacg                                                                                                                                                                                                                                           <  5:748486/114‑1 (MQ=255)
                                                                                                                                                   ttGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAggg                                                                                                                                                                                                                       <  7:350821/133‑1 (MQ=255)
                                                                                                                                                   ttGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAggg                                                                                                                                                                                                                       >  8:350821/1‑133 (MQ=255)
                                                                                                                                                           ggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGt                                                                                                                                                                                                >  5:1055336/1‑148 (MQ=255)
                                                                                                                                                             ttATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCttt                                                                                                                                                                                                                               >  6:1105218/1‑115 (MQ=255)
                                                                                                                                                             ttATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCttt                                                                                                                                                                                                                               <  5:1105218/115‑1 (MQ=255)
                                                                                                                                                                 gTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTa                                                                                                                                                                                                                                                >  8:48364/1‑94 (MQ=255)
                                                                                                                                                                 gTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTa                                                                                                                                                                                                                                                <  7:48364/94‑1 (MQ=255)
                                                                                                                                                                         gTCGGTAAAACATTTTAATCCTGAATTTGTTTGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGa                                                                                                                              <  1:159299/200‑1 (MQ=255)
                                                                                                                                                                              tAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATgg                                                                                                                                                                             >  6:866185/1‑148 (MQ=255)
                                                                                                                                                                               aaaaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCt                                                                                                                                                                                                                                 >  8:1115445/1‑95 (MQ=255)
                                                                                                                                                                               aaaaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCt                                                                                                                                                                                                                                 <  7:1115445/95‑1 (MQ=255)
                                                                                                                                                                                  aCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACAtgattga                                                                                                                                                                                                                                                                      <  6:44718/55‑1 (MQ=255)
                                                                                                                                                                                  aCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACAtgattga                                                                                                                                                                                                                                                                      >  5:44718/1‑55 (MQ=255)
                                                                                                                                                                                  aCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATAtt                                                                                                                                                                                                              <  7:474739/111‑1 (MQ=255)
                                                                                                                                                                                  aCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATAtt                                                                                                                                                                                                              >  8:474739/1‑111 (MQ=255)
                                                                                                                                                                                  aCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTc                                                                                                                                                                         <  5:866185/148‑1 (MQ=255)
                                                                                                                                                                                     ttttAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCtt                                                                                                                                                                                                                                <  8:165759/90‑1 (MQ=255)
                                                                                                                                                                                     ttttAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCtt                                                                                                                                                                                                                                >  7:165759/1‑90 (MQ=255)
                                                                                                                                                                                         aaTCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTtc                                                                                                                                                                  <  6:1055336/148‑1 (MQ=255)
                                                                                                                                                                                             cTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTGCGTTGTTTAAGAATTTCCCGCGCCGAATCTGcc                               <  2:412596/275‑1 (MQ=255)
                                                                                                                                                                                              tGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAg                                                                                                                                                                                                                         >  8:440404/1‑88 (MQ=255)
                                                                                                                                                                                              tGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAg                                                                                                                                                                                                                         <  7:440404/88‑1 (MQ=255)
                                                                                                                                                                                              tGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTtcatca                                                                                                                                                              >  7:716377/1‑147 (MQ=255)
                                                                                                                                                                                              tGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTtcatca                                                                                                                                                              <  8:716377/147‑1 (MQ=255)
                                                                                                                                                                                                 aTTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTc                                                                                                                                                                                               <  6:557436/111‑1 (MQ=255)
                                                                                                                                                                                                 aTTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTc                                                                                                                                                                                               >  5:557436/1‑111 (MQ=255)
                                                                                                                                                                                                       tCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGCTCTTCACTCTGTCATATCCCCGGTAATATGGCGTCAGGCTtcatcatc                                                                                                                                                            <  6:1074570/140‑1 (MQ=255)
                                                                                                                                                                                                       tCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTtcatcatc                                                                                                                                                            >  5:752629/1‑140 (MQ=255)
                                                                                                                                                                                                       tCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTtcatcatc                                                                                                                                                            <  6:752629/140‑1 (MQ=255)
                                                                                                                                                                                                       tCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCCCCGGTAATATGGCGTCAGGCTtcatcatc                                                                                                                                                            >  5:1074570/1‑140 (MQ=255)
                                                                                                                                                                                                          tGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATc                                                                                                                                                                                                          >  5:1056179/1‑91 (MQ=255)
                                                                                                                                                                                                          tGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATc                                                                                                                                                                                                          <  6:1056179/91‑1 (MQ=255)
                                                                                                                                                                                                                        tAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGata                                                                                                                                                                                                                <  6:555807/71‑1 (MQ=255)
                                                                                                                                                                                                                        tAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGata                                                                                                                                                                                                                >  5:555807/1‑71 (MQ=255)
                                                                                                                                                                                                                          aaaaTACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGa                                                                                                                                                                                                            >  7:1056216/1‑73 (MQ=255)
                                                                                                                                                                                                                          aaaaTACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGa                                                                                                                                                                                                            <  8:1056216/73‑1 (MQ=255)
                                                                                                                                                                                                                          aaaaTACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTCCTGCGTTGTTGAAGAATTTCTCGCGCCAAATCTGCCATTATCATGACTGCACCTGCCAACATCag  >  2:119748/1‑275 (MQ=255)
                                                                                                                                                                                                                               aCATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCa                                                                                                                                                                        <  6:797305/104‑1 (MQ=255)
                                                                                                                                                                                                                               aCATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCa                                                                                                                                                                        >  5:797305/1‑104 (MQ=255)
                                                                                                                                                                                                                               aCATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAgg                                                                                                                                                     <  5:596286/123‑1 (MQ=255)
                                                                                                                                                                                                                               aCATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAgg                                                                                                                                                     >  6:596286/1‑123 (MQ=255)
                                                                                                                                                                                                                               aCATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCAc                                                                                                                                                                                                     <  6:185024/75‑1 (MQ=38)
                                                                                                                                                                                                                               aCATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCAc                                                                                                                                                                                                     >  5:185024/1‑75 (MQ=38)
                                                                                                                                                                                                                               aCATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAgg                                                                                                                                                     <  5:875160/123‑1 (MQ=255)
                                                                                                                                                                                                                               aCATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAgg                                                                                                                                                     >  6:875160/1‑123 (MQ=255)
                                                                                                                                                                                                                                     ttgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAatat                                                                                                                                                                               <  6:1064739/91‑1 (MQ=255)
                                                                                                                                                                                                                                     ttgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAatat                                                                                                                                                                               >  5:1064739/1‑91 (MQ=255)
                                                                                                                                                                                                                                     ttgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATgg                                                                                                                                                                             <  8:653129/93‑1 (MQ=255)
                                                                                                                                                                                                                                     ttgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCa                                                                                                                                                                        <  7:921024/98‑1 (MQ=255)
                                                                                                                                                                                                                                     ttgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCa                                                                                                                                                                        >  8:921024/1‑98 (MQ=255)
                                                                                                                                                                                                                                     ttgatATTAAGCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATgg                                                                                                                                                                             >  7:653129/1‑93 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                     
GTTACGTTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTGCGTTGTTTAAGAATTTCCCGCGCCGAATCTGCCATTATCATGACTGCACCTGCCAGCATCAG  >  NC_000913/317335‑317827

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