Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,792,768 A→G 8.1% intergenic (‑217/+58) kbl ← / ← yibB glycine C‑acetyltransferase/YibB family protein, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,792,7680AG 8.1% 557.2 / 22.5 171intergenic (‑217/+58)kbl/yibBglycine C‑acetyltransferase/YibB family protein, function unknown
Reads supporting (aligned to +/- strand):  ref base A (77/79);  new base G (7/7);  total (85/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 185 total aligned reads displayed.

TTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTT  >  NC_000913/3792622‑3792913
                                                                                                                                                  |                                                                                                                                                 
ttAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGc                                                                                                                                                 >  1:1272224/1‑149 (MQ=255)
      cTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCAtt                                                                                                                                           <  2:1249529/149‑1 (MQ=255)
       ttCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAttt                                                                                                                                          <  1:828530/149‑1 (MQ=255)
          tGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTc                                                                                                                                       <  2:1272224/149‑1 (MQ=255)
            cAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCt                                                                                                                                     >  4:837571/1‑149 (MQ=255)
                  tCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTCTTCCTATCTCa                                                                                                                               <  3:837571/149‑1 (MQ=255)
                     ggATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCAtttt                                                                                                                                         >  4:690072/1‑136 (MQ=255)
                     ggATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCAtttt                                                                                                                                         <  3:690072/136‑1 (MQ=255)
                      gATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   <  4:295250/125‑1 (MQ=255)
                      gATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   <  4:1437046/125‑1 (MQ=255)
                      gATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   >  3:1070572/1‑125 (MQ=255)
                      gATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   <  4:1070572/125‑1 (MQ=255)
                      gATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   >  3:1437046/1‑125 (MQ=255)
                      gATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   >  3:295250/1‑125 (MQ=255)
                      gATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTg                                                                                                                           >  1:1025933/1‑149 (MQ=255)
                            taaCGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTc                                                                                                                                       <  1:638921/131‑1 (MQ=255)
                            taaCGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTa                                                                                                                                    <  3:988620/134‑1 (MQ=255)
                            taaCGTGTAACTCGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCAtttt                                                                                                                                         <  2:1330391/129‑1 (MQ=255)
                            taaCGTGTAACTCGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCAtttt                                                                                                                                         >  1:1330391/1‑129 (MQ=255)
                                 tgtAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   <  3:714861/114‑1 (MQ=255)
                                 tgtAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   >  4:714861/1‑114 (MQ=255)
                                               tctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGcatcat                                                                                                                                            <  1:308424/107‑1 (MQ=255)
                                               tctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGcatcat                                                                                                                                            >  2:308424/1‑107 (MQ=255)
                                               tctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATgg                                                                                                  >  1:1077247/1‑149 (MQ=255)
                                                 tcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAttt                                                                                                                                          <  1:895130/107‑1 (MQ=255)
                                                   aaCACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGc                                                                                                                                                 <  3:899290/98‑1 (MQ=255)
                                                    acacCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGa                                                                                             <  2:1025933/149‑1 (MQ=255)
                                                    acacCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGa                                                                                             <  1:1446620/149‑1 (MQ=255)
                                                        cTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAtt                                                                                                                                           >  1:883281/1‑99 (MQ=255)
                                                        cTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGcatca                                                                                                                                             >  1:18572/1‑97 (MQ=255)
                                                        cTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCAttttt                                                                                                                                        <  1:710969/102‑1 (MQ=255)
                                                        cTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCAttttt                                                                                                                                        >  2:710969/1‑102 (MQ=255)
                                                        cTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCAttt                                                                                                                                          >  2:1111867/1‑100 (MQ=255)
                                                        cTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTg                                                                                                                           <  4:1403778/115‑1 (MQ=255)
                                                        cTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATc                                                                                                                      <  1:624290/120‑1 (MQ=255)
                                                        cTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAAcc                                                                                         >  3:517097/1‑149 (MQ=255)
                                                        cTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAAc                                                                                          >  4:396438/1‑148 (MQ=255)
                                                        cTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAAc                                                                                          <  3:396438/148‑1 (MQ=255)
                                                         tACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGca                                                                                                                                                <  3:654488/93‑1 (MQ=255)
                                                         tACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGca                                                                                                                                                >  4:654488/1‑93 (MQ=255)
                                                            aGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGca                                                                                                                                                <  3:743230/90‑1 (MQ=255)
                                                            aGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGca                                                                                                                                                >  4:743230/1‑90 (MQ=255)
                                                             gATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGaaa                                                                                           >  2:406176/1‑142 (MQ=255)
                                                                  tagAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGt                                                                                                                         <  2:567825/107‑1 (MQ=255)
                                                                    gAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCAttttt                                                                                                                                        <  4:26749/90‑1 (MQ=255)
                                                                    gAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCAttttt                                                                                                                                        >  3:26749/1‑90 (MQ=255)
                                                                    gAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGt                                                                                                                             >  3:1011304/1‑101 (MQ=255)
                                                                    gAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCtt                                                                                >  4:911510/1‑146 (MQ=255)
                                                                    gAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACGGCATTGATATAAAAAGAAAAAATGATGCCAAAGCGTCATTTTTCCTATCTCAGt                                                                                                                             <  4:1011304/101‑1 (MQ=255)
                                                                     aaaaaCTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaaa                                                                            <  4:517097/149‑1 (MQ=255)
                                                                               cAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAg                                                                                                <  1:129019/119‑1 (MQ=255)
                                                                               cAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAg                                                                                                >  2:129019/1‑119 (MQ=255)
                                                                               cAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaaa                                                                            <  3:115958/139‑1 (MQ=255)
                                                                               cAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaaa                                                                            >  4:115958/1‑139 (MQ=255)
                                                                                         tGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGt                                                                                                                             >  2:828786/1‑80 (MQ=255)
                                                                                              aTAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCt                                                                                                                                     >  2:663126/1‑67 (MQ=255)
                                                                                              aTAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCt                                                                                                                                     <  1:663126/67‑1 (MQ=255)
                                                                                               tAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   >  3:1344225/1‑52 (MQ=255)
                                                                                               tAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   <  2:887958/52‑1 (MQ=255)
                                                                                               tAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaaa                                                                                                                                                   <  4:1344225/52‑1 (MQ=255)
                                                                                                aaTGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTc                                                                                  >  2:369620/1‑116 (MQ=255)
                                                                                                 atgatgATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAttt                                                                                                                                          >  1:560579/1‑59 (MQ=25)
                                                                                                 atgatgATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCAtttt                                                                                                                                         <  1:829311/60‑1 (MQ=255)
                                                                                                   gatgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAtt                                                                                                                                           <  3:416182/56‑1 (MQ=25)
                                                                                                   gatgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAtt                                                                                                                                           >  4:416182/1‑56 (MQ=25)
                                                                                                   gatgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaaa                                                                            <  1:1148314/119‑1 (MQ=255)
                                                                                                   gatgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaaa                                                                            >  2:1148314/1‑119 (MQ=255)
                                                                                                    atgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaaa                                                                            <  2:965816/118‑1 (MQ=255)
                                                                                                    atgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaaa                                                                            >  1:965816/1‑118 (MQ=255)
                                                                                                       attgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTaa                                          >  3:202915/1‑149 (MQ=255)
                                                                                                        ttgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAAc                                                                                                                 >  3:1202991/1‑77 (MQ=255)
                                                                                                        ttgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAAc                                                                                                                 <  4:1202991/77‑1 (MQ=255)
                                                                                                        ttgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTAtt                                                                                     >  1:608423/1‑105 (MQ=255)
                                                                                                        ttgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTAtt                                                                                     <  2:608423/105‑1 (MQ=255)
                                                                                                              ttACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCCGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTa                                   <  4:202915/149‑1 (MQ=255)
                                                                                                                     aTTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAg                                                                                                <  3:1055063/81‑1 (MQ=255)
                                                                                                                     aTTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCtt                                                                                <  2:1200269/97‑1 (MQ=255)
                                                                                                                     aTTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGAAg                                                                                                >  4:1055063/1‑81 (MQ=255)
                                                                                                                     aTTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGGTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGaaa                                                                       >  2:634487/1‑106 (MQ=255)
                                                                                                                       tGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAttttt                                                  >  2:642965/1‑125 (MQ=255)
                                                                                                                         atatAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGaaaa                                                                      >  3:602647/1‑103 (MQ=255)
                                                                                                                          tataAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCtt                                                                                <  3:320954/92‑1 (MQ=255)
                                                                                                                          tataAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCtt                                                                                >  4:320954/1‑92 (MQ=255)
                                                                                                                           ataAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgt                                       <  1:392519/132‑1 (MQ=255)
                                                                                                                            taAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAttt                                                    >  1:84259/1‑118 (MQ=255)
                                                                                                                             aaaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGcc                                               >  3:753028/1‑122 (MQ=255)
                                                                                                                             aaaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGcc                                               <  4:753028/122‑1 (MQ=255)
                                                                                                                              aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTAc                                                                                                            <  1:1237736/60‑1 (MQ=255)
                                                                                                                              aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTa                                   <  2:483937/133‑1 (MQ=255)
                                                                                                                              aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTa                                   >  1:483937/1‑133 (MQ=255)
                                                                                                                                aaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCa                                                                <  4:794797/102‑1 (MQ=255)
                                                                                                                                aaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCa                                                                >  3:794797/1‑102 (MQ=255)
                                                                                                                                 aggaaaaaTGATGCTTTGGCATCATTTTTTCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTc                                                                                  >  2:89632/4‑83 (MQ=25)
                                                                                                                                     aaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaaa                                                                            <  4:487532/85‑1 (MQ=255)
                                                                                                                                     aaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaaa                                                                            >  3:487532/1‑85 (MQ=255)
                                                                                                                                     aaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTa                                   >  2:833816/1‑126 (MQ=255)
                                                                                                                                      aaaTGATGCTTTGGCATCATTTTTTCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAAtt                                                                    <  4:626293/92‑1 (MQ=38)
                                                                                                                                            tgCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAttttt                                                  >  4:570703/1‑104 (MQ=255)
                                                                                                                                             gCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGaa         <  2:195330/144‑1 (MQ=255)
                                                                                                                                                 aaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATc            <  3:137325/137‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTT  >  NC_000913/3792622‑3792913

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