Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,287,413 G→A intergenic (+388/‑13) agaI → / → yraH putative deaminase AgaI/putative fimbrial protein YraH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,287,4130GA81.1% 63.8 / 11.8 39intergenic (+388/‑13)agaI/yraHputative deaminase AgaI/putative fimbrial protein YraH
Reads supporting (aligned to +/- strand):  ref base G (5/2);  new base A (18/12);  total (25/14)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

CTATTAATAAATCACCTCGTATTCTGACTCGCCCCCTGGCGCAGCAAAGTCAGGCATTTATACTCTTTTGAATACATAGAAAATTGATATCAATATAATGAAAATATCAAATATTTGCATATAAATATAATCTTAAAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGGATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTTTGATGGCGAGATTTTAAAATCCGCTTGTGAAATCAATGACTCTGATAAGAAAATTGAAGTTGCTCTTGGTCACTATAATGCCGAACAGTTTCGTAATATTGGCGAGCGTAGCCCAAAAATTCCATTCACTATTCCTTTAGTTAACTGCCCA  >  NC_000913/3287175‑3287662
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                         
cTATTAATAAATCACCTCGTATTCTGACTCGCCCCCTGGCGCAGCAAAGTCAGGCCTTTATACTCTTTTGAATACATAGAAAATTGATATCAATATAATGAAAATATCAAATATTTGCATATAAATATAATCTTAAAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGAATAAGGTtatat                                                                                                                                                                                                                                               >  2:462438/1‑251 (MQ=255)
  aTTAATAAATCACCTCGTATTCTGACTCGCCCCCTGGCGCAGCAAAGTCAGGCATTTATACTCTTTTGAATACATAGAAAATTGATATCAATATAATGAAAATATCAAATATTTGCATATAAATATAATCTTAAAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGGATAAGGTTatatat                                                                                                                                                                                                                                             >  2:533915/1‑251 (MQ=255)
  aTTAATAAATCACCTCGTATTCTGACTCGCCCCCTGGCGCAGCAAAGTCAGGCATTTATACTCTTTTGAATACATAGAAAATTGATATCAATATAATGAAAATATCAAATATTTGCATATAAATATAATCTTAAAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGAATAAGGTtata                                                                                                                                                                                                                                                <  1:77787/248‑1 (MQ=255)
    taataaATCACCTCGTATTCTGACTCGCCCCCTGGCGCCGCAAAGTCAGGCATTTATACTCTTTTGAATACATAGAACATTGATATCAATATAATGCAAATATCAAATATTTGCATATAAATATAATCTTACAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATACTCTTTCAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGAATAAGGTTATATATg                                                                                                                                                                                                                                            >  2:580344/1‑250 (MQ=255)
            cACCCCGTATTCTGACTCGCCCCCTGGCGCAGCAAAGTCAGGCATTTATACTCTTTTGAATACATAGAAAATTTATATCAATATAATGAAAATATCAAATATTTGCATATAAATATAATCTTAAAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGAGGTTTAACATATATTCACATGAAATATGATTAGGTACGTGTTACACGAATAAGGTTATATATGAATAAAGt                                                                                                                                                                                                                                    <  2:600949/250‑1 (MQ=255)
                              gCCCCCTGGCGCAGCAAAGTCAGGCATTTATACTCTTTTGAATACATAGAAAATTGATATCAATATAATGAAAATATCAAATATTTGCATATAAATATAATCTTAAAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACTAGAATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCg                                                                                                                                                                                                                 <  1:372346/251‑1 (MQ=255)
                                                atCAGGCATTTATACTCTTTTGAATACATAGAAAATTGATATCAATATAATGAAAATATCAAATATTTGCATATAAATATAATCTTAAATATCAGTCTATTTAATGTTCAATGAAATATTTCTGCATGTATAATCTTTAAAGATGTTAAACATATATTCACATTAAATATGATTATATACTTGTTACAATAATAAGGTTATATTTGATGAAAGTTACAAAAACACCTATTGCGGGCTTCCTTGCACttttt                                                                                                                                                                                               <  2:72880/250‑2 (MQ=255)
                                                     gCATTTATACTCTTTTGAATACATAGAAAATTGATATCAATATAATGAAAATATCAAATATTTGCATATAAATATAATCTTAAAGTTCAGTCTATTTAATGTTCGATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGGATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGccgg                                                                                                                                                                                          >  2:584186/1‑248 (MQ=255)
                                                          tataCTCTTTTGAATACATAGAAAATTGATATCAATATAATGAAAATATCAAATATTTGCATATAAATATAATCTTAAAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGAATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATg                                                                                                                                                                                     >  2:42313/1‑251 (MQ=255)
                                                                               aaaaTTGATATCAATATAATGAAAATATCAAATATTTGCATATAAATATAATCTTAAAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGAATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGaaa                                                                                                                                                                <  1:66432/251‑1 (MQ=255)
                                                                                aaaTTGATATCAATATAATGAAAATATCAAATATTTGCATATAAATATAATCTTAAAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGAATAAGGTTATATATGAATACAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCGTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAAt                                                                                                                                                               <  2:27357/251‑1 (MQ=255)
                                                                                 aaTTGATATCAATATAATGAAAATATCAAATATTTGCATATAAATATAATCTTAAAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGGATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATa                                                                                                                                                              <  2:170048/251‑1 (MQ=255)
                                                                                  aTTGATATCAATATAATGAAAATATCAAATATTTGCATATAAATATAATCTTAAAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGAATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAAt                                                                                                                                                               >  1:7720/1‑249 (MQ=255)
                                                                                     gATATCAATATAATGAAAATATCAAATATTTGCATATAAATATAATCTTAAAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGAGTATGGACTTGTTACAAGGAGAAGGTTATATATGAATAAAGTTACAAAAACAGCGATTGCGGGGTTACTTGCGCTGTTCGCAGGTAAGGCCGCTGCAACCGATGGTTAAATAGttt                                                                                                                                                          <  2:361668/251‑1 (MQ=255)
                                                                                       tatCAATATAATGAAAATATCAAATATTTGCATATAAATATAATCTTAAAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGGATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGAAGGTAATGCCGCTGCAACCGATGGTGAAATAGttttt                                                                                                                                                        >  1:306/1‑251 (MQ=255)
                                                                                           aaTATAATGAAAATATCAAATATTTGCATATAAATATAATCTTAAAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGGATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTTTGATg                                                                                                                                                    >  1:41137/1‑251 (MQ=255)
                                                                                                aaTGAAAATATCAAATATTTGCATATAAATATAATCTTAAAGTTAAGTCTATTTAATGTTCAATGCAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGAATAAGGTTATATATGACTAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTTTGATGGCg                                                                                                                                                 >  1:361591/1‑249 (MQ=255)
                                                                                                         atCAAATATTTGCATATAAATATAATCTTAAAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGAATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTTTGATGGCGAGATTTTaaaa                                                                                                                                      <  1:580344/251‑1 (MQ=255)
                                                                                                         atCAAATATTTGCATATAAATATAATCTTAAAGTTCAGTCTATTTAATGTTCAATGAAATATGTCTGCCTGTATAAGCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGAATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTCTGCTGGCGAGATTTTaaaa                                                                                                                                      <  2:27126/251‑1 (MQ=255)
                                                                                                          tCAAATATTTGCATATAAATATAATCTTACAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACCTATATTCACATTAAATATGATTATGTACTTGTTACAAGAATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATCGTTTTTGATGGCGAGATTTTAAAAt                                                                                                                                     >  2:589700/1‑251 (MQ=255)
                                                                                                               tatTTGCAACTACATATACTGTTAAAGTTCTGTCTATTTAATGTTCAATGCAATATCTATCCCTGTATAACCTTCAACGATGTTGAACCTATATTCACTTTCAATATGAATATGTAATTGTTACACGAATAAGGTTATATATGAATAACGTTACCAACACAGCTATGGAGGGCTTACTTGCGCTTTTCGAGGGTAATGCGCCAGCTACCGATGGTGATATAGTTTTTGATGGCGAGATTTTAAAATCCGCt                                                                                                                                >  2:296350/1‑251 (MQ=255)
                                                                                                                atTTGCATATAAATATAATCTTAAAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGAATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTTTGATGGCGAGATTTTAAAATCCGCtt                                                                                                                               >  1:114519/1‑251 (MQ=255)
                                                                                                                  ttGCATATAAATATAATCTTAAAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGAATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTTTGATGGCGAGCTTTTAAAATCCGCTtgt                                                                                                                             >  1:162905/1‑251 (MQ=255)
                                                                                                                            atataaTCTTAAAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGAATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTTTGATGGCGAGATTTTAAAATCCGCTTGTGAAATCAATg                                                                                                                   >  1:513603/1‑251 (MQ=255)
                                                                                                                               taaTCTTAAAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGAATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTTTGATGGCGAGATTTTAAAATCCGCTTGTGAAATCAATGAct                                                                                                                >  1:460758/1‑251 (MQ=255)
                                                                                                                                aaTCTTAAAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGAATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTTTGATGGCGAGATTTTAAAATCCGCTTGTGAAATCAATGActc                                                                                                               <  2:483144/251‑1 (MQ=255)
                                                                                                                                 aTCTTAAAGTTCAGTCTATTTAATGTTCAATGAACTATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGAATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGAGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTTTGATGGCGAGATTTTACAATCCGCTTGTGAAATCAATGActct                                                                                                              >  1:230382/1‑251 (MQ=255)
                                                                                                                                      aaaGTACACTCTATTTAATGTTCAATGAATTATTTCTGACTGTATACTCTATAAAGATGTTGAACATATATTCACATTAAGTCTGATTATGTACTTGTTACACGAATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGAGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTTTGATGGCGAGATTTTAAAATCCGCTTGTGAAATAAATGACTCTGATa                                                                                                          >  2:593641/1‑250 (MQ=255)
                                                                                                                                                 tATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGAATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTTTGATGGCGAGATTTTAAAATCCGCTTGTGAAATCAATGACTCTGATAAGAAAATTGaa                                                                                               <  2:201023/250‑1 (MQ=255)
                                                                                                                                                 tATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGAATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTTTGATGGCGAGATTTTAAAATCCGCTTGTGAAATCAATGACTCTGATAAGAAAATTGa                                                                                                >  1:68703/1‑249 (MQ=255)
                                                                                                                                                           tCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGCATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTTTGATGGCGAGATTTTAAAATCCGCTTGTGAAATCAATGACTCTGATAAGAAAATTGAAGTTGCTCTTgg                                                                                    >  2:74502/1‑251 (MQ=255)
                                                                                                                                                                                      cTTTAAAGCTGTTGAACATAAATTCACACTAATTATAATTATGTACTTGTTACAAGAATCAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGAGGGCTTAATTGCGCTTTTCGCAGGTAAAGCAGCTGCAAACGATGGTGAAATAGTTTTTGATGGCGCGCTTTTAAAATCCGCTTGTGAAATCAATGACTCTGCTAAGAAAATTGAAGTTGCTCTTGGTCACTATAATGCCGAACAGTTTCgaa                                                          >  2:628882/1‑248 (MQ=255)
                                                                                                                                                                                      atttAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGCATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGTCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTTTGATGACGCGATTTTAAAATCCGCTTGTGAAATCAATGACTCTGATAAGAAAATTGAAGTTGCTCTTGGTCACTATAATGCCGAACAGTTTCGTa                                                          >  2:22103/2‑250 (MQ=255)
                                                                                                                                                                                                gTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGAATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTTTGATGGCGAGATTTTAAAATCCGCTTGTGAAATCAATGACTCTGATAAGAAAATTGAAGTTGCTCTTGGTCACTATAATGCCGAACAGTTTCGTAATATTGGCg                                                 >  1:116173/1‑249 (MQ=255)
                                                                                                                                                                                                   gAACATATATTCACATTAAATATGATTATGGACTTGTTACAAGAATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTTTGATGGCGAGATTTTAAAATCCGCTTGTGAAATCAATGACTCTGATAAGGAAATTGAAGTTGCTCTTGGTCACTATAATGCCGCACAGTTTCGTAATATTGGCGAGCgg                                            >  2:42374/1‑250 (MQ=255)
                                                                                                                                                                                                                caTTAAATATGATTATGTACTTGTTACAAGAATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTGTTTGATGGCGAGATTTTAAAATCCGCTTGTGAAATCTATGACTCTGATAAGAAAATTGAAGTTGCTCTTGGTCACTATAATGCCGAACAGTTTCGTAATATTGGCGAGCGTAGCCCAAAAATTc                               <  1:416771/251‑1 (MQ=255)
                                                                                                                                                                                                                         tGATTATGTACTTGTTACAAGAATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTTTGATGGCGAGATTTTAAAATCCGCTTGTGAAATCAATGACTCTGATAAGAAAATTGAAGTTGCTCTTGGTCACTATAATGCCGAACAGTTTCGTAATATTGGCGAGCGTAGCCCAAAAATTCCATTCACt                       >  1:516639/1‑250 (MQ=255)
                                                                                                                                                                                                                                gTACGTGTTACAAGAATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTTTGATGGCGAGATTTTAAAATCCGCTTGTGAAATCAATGACTCTGATAAGAAAATTGAAGTTGCTCTTGGTCACTATAATGCCGAACAGTTTCGTAATATTGGCGAGCGTAGCCCAAAAATTCCATTCTCTATTCCtt                <  2:415562/250‑1 (MQ=255)
                                                                                                                                                                                                                                     tGTTAAAAGGATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTTTGATGGCGAGATTTTAAAATCCGCTTGTGAAATCAATGACTCTGATAAGAAAATTGAAGTTGCTCTTGGTCACTATAATGCCGAACAGTTTCGTAATATTGGCGAGCGTAGCCCAAAAATTCCATTCACTATTCCTTTAGtt           >  1:556401/1‑250 (MQ=255)
                                                                                                                                                                                                                                              gATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCATACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTTTGATGGCGAGATTTTAAAATCCGCTTGTGAAATCAATGACTCTGATAAGAAAATTGAAGTTGCTCTTGGTCACTATAATGCCGAACAGTTTCGTAATATTGGCGAGAGTAGCCCAAAAATTCCATTCACTATTCCTTTAGTTAACTGCCCa  >  1:576085/1‑250 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                         
CTATTAATAAATCACCTCGTATTCTGACTCGCCCCCTGGCGCAGCAAAGTCAGGCATTTATACTCTTTTGAATACATAGAAAATTGATATCAATATAATGAAAATATCAAATATTTGCATATAAATATAATCTTAAAGTTCAGTCTATTTAATGTTCAATGAAATATTTCTGCCTGTATAATCTTTAAAGATGTTGAACATATATTCACATTAAATATGATTATGTACTTGTTACAAGGATAAGGTTATATATGAATAAAGTTACAAAAACAGCTATTGCGGGCTTACTTGCGCTTTTCGCAGGTAATGCCGCTGCAACCGATGGTGAAATAGTTTTTGATGGCGAGATTTTAAAATCCGCTTGTGAAATCAATGACTCTGATAAGAAAATTGAAGTTGCTCTTGGTCACTATAATGCCGAACAGTTTCGTAATATTGGCGAGCGTAGCCCAAAAATTCCATTCACTATTCCTTTAGTTAACTGCCCA  >  NC_000913/3287175‑3287662

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