Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,581,138 C→T 6.5% G109S (GGT→AGT)  hsdS ← type I restriction enzyme EcoKI specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,581,1380CT 6.5% 178.8 / 2.5 62G109S (GGT→AGT) hsdStype I restriction enzyme EcoKI specificity protein
Reads supporting (aligned to +/- strand):  ref base C (32/25);  new base T (2/2);  total (35/27)
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.21e-01

CAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACA  >  NC_000913/4581001‑4581280
                                                                                                                                         |                                                                                                                                              
caaatcaaaGCTTGCCGGCTTAATATTATTTATATTTGCACCGGCAGAAAGTGATGAAATTGTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCaaa                                                                                                                                           >  1:945583/1‑143 (MQ=255)
caaatcaaaGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCaaa                                                                                                                                           >  1:1006075/1‑143 (MQ=255)
 aaatcaaaGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAACCACAGAGGCTATCAGTTTTTCAGGACGTAATACACCGCaaaa                                                                                                                                          >  2:899755/1‑143 (MQ=255)
  aatcaaaGCTTGCCGGCTTAATATTATTAATATTTGCAGCAGCAGAAAGTGATGAAATTTTGTTTCGAGAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAAt                                                                                                                                         <  2:1040310/143‑1 (MQ=255)
  aatcaaaGCCTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAAt                                                                                                                                         <  1:379017/143‑1 (MQ=255)
   atcaaaGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATg                                                                                                                                        <  1:975082/143‑1 (MQ=255)
    tcaaaGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATgc                                                                                                                                       >  1:76282/1‑143 (MQ=255)
             gCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGGAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACt                                                                                                                              >  2:864993/1‑143 (MQ=255)
               cGGCTTAATATTATGAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAACCCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTAc                                                                                                                            <  2:1207294/143‑1 (MQ=255)
                  cTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAGTATAAGTGTTTCAGGACGTAATACACCGCAAAATGCGCCGACACTACAtt                                                                                                                         >  1:723827/1‑143 (MQ=255)
                     aaTATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCaa                                                                                                                      >  1:1012323/1‑143 (MQ=255)
                       tattatTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAAt                                                                                                                    <  1:954223/143‑1 (MQ=255)
                            ttaATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAg                                                                                                               >  2:1146868/1‑143 (MQ=255)
                            ttaATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAg                                                                                                               <  2:732322/143‑1 (MQ=255)
                             taATATTTGCACCAGCAGAAAGTGATGAAATTCTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGa                                                                                                              <  1:495502/143‑1 (MQ=255)
                               atatTGGCACCAGTAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATg                                                                                                            <  2:744626/143‑1 (MQ=255)
                                         cagcagAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACTGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGAtgtgtg                                                                                                  >  2:1097681/1‑141 (MQ=255)
                                         cagcagAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACTGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGAtgtgtg                                                                                                  >  2:1096742/1‑141 (MQ=255)
                                         cagcagAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCg                                                                                                  <  2:788596/143‑1 (MQ=255)
                                           gcagAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGa                                                                                                <  2:293246/143‑1 (MQ=255)
                                           gcagAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGa                                                                                                <  2:338956/143‑1 (MQ=255)
                                               aaaGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTa                                                                                            <  1:125403/143‑1 (MQ=255)
                                               aaaGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCCGATGTGCGGAttga                                                                                            >  2:1259868/1‑141 (MQ=255)
                                                  gTGATGACATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCa                                                                                         >  1:908017/1‑143 (MQ=255)
                                                           ttttGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGAt                                                                                <  1:1008568/143‑1 (MQ=255)
                                                            tttGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGAtt                                                                               <  1:1072870/143‑1 (MQ=255)
                                                               gTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGc                                                                            <  1:516446/143‑1 (MQ=255)
                                                                 ttCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCCGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTc                                                                          >  2:634698/1‑143 (MQ=255)
                                                                  gcGATAAAGAGAAGATTTTGTGAAATGAGAAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTcc                                                                         <  2:1213887/142‑1 (MQ=255)
                                                                    gATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCt                                                                       >  1:643578/1‑143 (MQ=255)
                                                                    gATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCt                                                                       <  2:1034748/143‑1 (MQ=255)
                                                                     tctaaaagaagaTTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGAATTTACCAACTACGGATTTGCTCCCtg                                                                      >  2:141348/7‑143 (MQ=255)
                                                                      tAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCtga                                                                     >  2:14020/1‑143 (MQ=255)
                                                                       aaaGAGAAGAGTTCGTGAAATGAGCCATAGAACCAGAAAATATAAGTGTCTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGAGTTACCAACTACGGATTTGCTCCCtgat                                                                    <  2:950212/143‑1 (MQ=255)
                                                                         agagAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCtgatga                                                                  >  1:1019324/1‑143 (MQ=255)
                                                                           agaagaTTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACa                                                                >  1:582558/1‑143 (MQ=255)
                                                                             aagaTTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACAtt                                                              >  1:250368/1‑143 (MQ=255)
                                                                              agaTTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTg                                                             >  1:546665/1‑143 (MQ=255)
                                                                                aTTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTCGAAGCTGATGTGCGGATTTACCAACTACGGCTTTGCTCCCCAATGACAtcaca                                                           >  2:616588/1‑139 (MQ=255)
                                                                                 ttttGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCaa                                                          >  2:399799/1‑143 (MQ=255)
                                                                                   tggggAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAAGACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATa                                                        <  2:316913/139‑1 (MQ=255)
                                                                                    ggtgAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATaa                                                       <  1:490571/142‑1 (MQ=255)
                                                                                         aaTGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATCGCAATACCAata                                                  >  1:229878/1‑143 (MQ=255)
                                                                                               cAATAAATCCATAAAATATAAGTTGTGCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCa                                            <  2:303189/143‑1 (MQ=255)
                                                                                                aaTAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAg                                           <  2:1012323/143‑1 (MQ=255)
                                                                                                 aTAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAgg                                          >  1:1137690/1‑143 (MQ=255)
                                                                                                  tAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGga                                         <  1:919821/143‑1 (MQ=255)
                                                                                                  tAAAACCAGAAAATATAAGTTTTTCAGGAAACAAAAGACAGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGga                                         >  1:878929/1‑143 (MQ=255)
                                                                                                        cAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATAtt                                   >  1:298192/1‑143 (MQ=255)
                                                                                                           aaaaTATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGAGGGGTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATATTCAGGGGATAttttt                                >  1:654545/1‑143 (MQ=255)
                                                                                                              atatAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGCGACATTTTTTGa                             >  2:1083412/1‑143 (MQ=255)
                                                                                                                   attttttCAGGACGTAATACACCGCATATTGCGCCGAAAATACAGTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGActttc                        <  2:140687/141‑1 (MQ=255)
                                                                                                                     tttttCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGActttctt                      >  1:185962/1‑143 (MQ=255)
                                                                                                                     tttttCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGActttctt                      <  1:629411/143‑1 (MQ=255)
                                                                                                                       tttCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATGTTCAGGAGATCTCTTTTGACTTTCTTTa                    >  2:554675/1‑143 (MQ=255)
                                                                                                                           aGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACaa                >  1:1113251/1‑143 (MQ=255)
                                                                                                                            ggACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAg               >  1:1148981/1‑143 (MQ=255)
                                                                                                                               cGTAATACACCGCAAAATGCGACGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGAATTGCACCCTGATGACATTGCAATAACAATAACGTCAGGAGATATTTTTAGACTTTCGTTaaggagaat            >  2:1057318/1‑136 (MQ=255)
                                                                                                                                   aTACACTGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGTGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTa        <  1:1097681/143‑1 (MQ=255)
                                                                                                                                   aTACACTGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGTGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTa        <  1:1096742/143‑1 (MQ=255)
                                                                                                                                     acacCCCAAAATGCGCCGAAACTACCTTCAAATGGTAGCTGCTGATGTGCGGAATTACCAACTACAGATTTCCTCCCTGATGACATTGCAATAACAATATCCTCAAGAGATACTTTTTGACTTTCTTTTACAAAATTTTTAgg      >  2:1213470/1‑143 (MQ=255)
                                                                                                                                      cacCGCAAAATGCGCCGAAACTACATTCAAATGGTAGCTGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAGTATCTTCAGGAGATATTTTTTGACTTTCTTTCACAACATTTATAgta     >  2:593249/1‑141 (MQ=255)
                                                                                                                                         cGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGaa    <  1:468019/141‑1 (MQ=255)
                                                                                                                                         cGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGaa    >  2:468019/1‑141 (MQ=255)
                                                                                                                                         cGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACa  >  1:890100/1‑143 (MQ=255)
                                                                                                                                         cGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGCGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTCCAATAACAATATCCTCAGCAGACATTTTTTCACTTTCTTTAACAAGATCCTTAGGaccc  >  2:325675/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                              
CAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACA  >  NC_000913/4581001‑4581280

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