Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,499,661 T→C intergenic (+161/‑1599) insD6 → / → yjgZ IS2 insertion element protein InsB/uncharacterized protein YjgZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,499,6610TC95.1% 117.4 / ‑3.8 41intergenic (+161/‑1599)insD6/yjgZIS2 insertion element protein InsB/uncharacterized protein YjgZ
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (15/24);  total (16/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

CATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCC  >  NC_000913/4499411‑4499897
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                            
cATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCtt                                                                                                                                                                                                                                              <  1:15825/251‑1 (MQ=255)
       cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAAcc                                                                                                                                                                                                                                       <  1:217964/251‑1 (MQ=255)
        gcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCCGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATCTTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAAcc                                                                                                                                                                                                                                       >  2:268495/1‑250 (MQ=255)
                                          cAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGAtt                                                                                                                                                                                                    <  1:426912/251‑1 (MQ=255)
                                             cgagCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGa                                                                                                                                                                                                  >  2:41752/4‑250 (MQ=255)
                                              gggCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGAAGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACGGTCTGCCGGTATATGCCTCGCTCCCTCAAAAACCATAACAGGTTAGTTATCTTTGTCTGTCAGCGGATTGAAt                                                                                                                                                                                                <  2:425110/251‑1 (MQ=39)
                                               agCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATg                                                                                                                                                                                               >  1:27752/2‑251 (MQ=255)
                                                              aaGTGATAACAGATGACTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTct                                                                                                                                                                                 >  2:307488/1‑249 (MQ=255)
                                                                       aagacggcTGGAAAAATAGGGGCAAAACCACCTTGGTCTGATATGAATACGGATATTTTATACAGAATCCCTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGCAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCCCTCTCCCTTAAAAACCATAACAGGGTAGCTATCTTCGTCTGGCAGCGGATTGAATGGCGGTATATTTCCAGTACAAATAc                                                                                                                                                                       <  2:236782/244‑1 (MQ=255)
                                                                         gATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCg                                                                                                                                                                     >  2:123742/1‑251 (MQ=255)
                                                                                  aaaTATAGGGTAAAAGCCGCCTTGGGCTGATATGAATACGGATATTTTATTGAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCAc                                                                                                                                                            <  2:254978/251‑1 (MQ=255)
                                                                                      ataGGGGCAAATCCACCCTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGa                                                                                                                                                         <  1:261037/250‑1 (MQ=255)
                                                                                          gggaaaaTCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGCTTATCTGTGATCAGGAAGAATAAGAGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACAATAAAAGGTTAGTTATCTTCGTCTGTCAGCGTCTTGAATGGCGGGATATTTTTAATACAAATAACGGTCCTCCCCTGacaaa                                                                                                                                                     >  2:39859/5‑250 (MQ=255)
                                                                                          gggCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGAtt                                                                                                                                                                                                    <  1:433184/203‑1 (MQ=255)
                                                                                          gggCAAATCAACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGAtt                                                                                                                                                                                                    >  2:433184/1‑203 (MQ=255)
                                                                                                      cTTGTGCTGCTATGAATACGGATATTTTCTCCAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGCTTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTCTATTCCTCTCTCCCTCAAAAACCATACCAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGTTATATTTTCAGTACAAATTCCGGACATCCACTAAAAAATACCTgacaaatg                                                                                                                                        >  2:316169/1‑244 (MQ=255)
                                                                                                         gtgCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACg                                                                                                                                     <  1:41752/251‑1 (MQ=255)
                                                                                                          tgCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCg                                                                                                                                                                                                                  >  1:207935/1‑173 (MQ=255)
                                                                                                          tgCTGATATGAATACGGAAATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCGGCCTCAAAAACCATAACAGGTTAGTTATCTTCg                                                                                                                                                                                                                  <  2:207935/173‑1 (MQ=255)
                                                                                                                       aCGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTc                                                                                                                       <  1:56985/251‑1 (MQ=255)
                                                                                                                           atatTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGAATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACt                                                                                                                   >  1:320138/1‑251 (MQ=255)
                                                                                                                             atTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCAAAATGACGGGCACTGATTTTTCTACTg                                                                                                                  >  1:345151/1‑250 (MQ=255)
                                                                                                                                  ataCAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGAATATTTTCAGTACCAATACCGGTCATCCACTGAAAAATACCTGCGAACAGGCGGGCACTGATTTTTCTACTGAAGGGc                                                                                                            >  2:120762/1‑251 (MQ=255)
                                                                                                                                        aaTACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCtgatga                                                                                                      <  1:241019/251‑1 (MQ=255)
                                                                                                                                           aCTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGAc                                                                                                   <  1:285790/251‑1 (MQ=255)
                                                                                                                                                 cAGTACAAGAGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGTAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTAt                                                                                              <  2:116016/250‑1 (MQ=255)
                                                                                                                                                  ataaCAAGGGACTGTATTGCAGTCGGCCATAGATGAGTACCTGTGATCAGGAAGAATCAGTGGCTGAGGTTTCAAACTGTAGGACGGTATATTCCTCTCTCCCACAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCTGTATATTTTCAGTACAAATTCCGGTCATCCACGGAAAAATCCCTGCGAAATGACGGGCACTGATTTTTATCCTGACGTGCTGATGAGCCGTTATGt                                                                                            <  2:88704/248‑1 (MQ=255)
                                                                                                                                                                     cAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATc                                                                           >  1:75610/1‑249 (MQ=255)
                                                                                                                                                                                      ttACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATAAGTggagcactttaacctc                                                        >  1:379288/1‑237 (MQ=255)
                                                                                                                                                                                         ccTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTaaaaatgcgtggagcaatgtaggctaccc                                                     >  1:323650/1‑228 (MQ=255)
                                                                                                                                                                                                  aGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCAGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGAc                                            <  2:448463/251‑1 (MQ=255)
                                                                                                                                                                                                            aGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAatgatg                                  <  1:22091/251‑1 (MQ=255)
                                                                                                                                                                                                                      gTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTg                        <  1:120762/251‑1 (MQ=255)
                                                                                                                                                                                                                                gTCTGCCTGTATATTCCTCTCTCCCTCAATAACCATAACAGGTTAGTTATCTTCGTCTGTAAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACTCTGTAGCCTCCCTGCACATGACCATAATGATGAACAGTATTGGCACTGTc                <  2:130054/249‑1 (MQ=255)
                                                                                                                                                                                                                                    gacGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGATATCTTCGTCTGTCAGCGGATTGAAGGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCTCTGTCGCTGAc          <  1:363174/249‑1 (MQ=255)
                                                                                                                                                                                                                                     ccGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCt                                                     <  2:118780/207‑1 (MQ=255)
                                                                                                                                                                                                                                     ccGGTATATGCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAACGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCTCTGAcc         <  2:460664/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGaaa                                                                                                                                                       <  1:94708/108‑1 (MQ=255)
                                                                                                                                                                                                                                      cGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGaaa                                                                                                                                                       >  2:94708/1‑108 (MQ=255)
                                                                                                                                                                                                                                      cGGTAAATTCCTCCCTCCCTCAAAAACCATAACAGGTAAGGTGTCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACTGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATTAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACTCTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCg        <  2:132441/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gTATATTCCCCTCCCCCTCAAAAACCACAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGAc      <  2:323650/251‑1 (MQ=255)
                                                                                                                                                                                                                                            atTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGcc  <  1:266598/251‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                            
CATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCC  >  NC_000913/4499411‑4499897

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