Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 84,735 A→G D123G (GAC→GGC)  leuO → global transcription factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091384,7350AG85.5% 201.3 / 18.7 76D123G (GAC→GGC) leuOglobal transcription factor
Reads supporting (aligned to +/- strand):  ref base A (3/8);  new base G (30/35);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 3.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.26e-01

ACCGACTGCTCGCGCATTTCAACTTTTTGGTTCAGTTCGTCAGGCATTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGACAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTAAATCAGAACACTGAACATCAGCTGCGTTATCAGGAAACGGAGTTTGTGAT  >  NC_000913/84598‑84864
                                                                                                                                         |                                                                                                                                 
aCCGACTGCTCGCGCATTTCAACTTTTTGGTTCAGTTCGTCAGGCATTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGAc                                                                                                                                  <  1:140990/139‑1 (MQ=255)
    aCTGCTCGCGCATTTCAACTTTTTGGTTCAGTTCGTCAGGCATTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGgcagca                                                                                                                              >  2:1357115/1‑139 (MQ=255)
      tGCTCGCGCATTTCAACTTTTTGGTTCAGTTCGTCAGGCATTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCAtt                                                                                                                            >  2:1601350/1‑139 (MQ=255)
      tGCTCGCGCATTTCAACTTTTTGGTTCAGTTCGTCAGGCATTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCAtt                                                                                                                            <  1:1695377/139‑1 (MQ=255)
      tGCTCGCGCATTTCAACTTTTTGGTTCAGTTCGTCAGGCATTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCAtt                                                                                                                            <  1:1357115/139‑1 (MQ=255)
         tCGCGCATTTCAACTTTTTGGTTCAGTTCGTCAGGCATTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTg                                                                                                                         >  1:1624907/1‑139 (MQ=255)
               aTTTCAACTTTTTGGTTCAGTTCGTCAGGCATTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCg                                                                                                                   <  1:743086/139‑1 (MQ=255)
                tttCAACTTTTTGGTTCAGTTCGTCAGGCATTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGc                                                                                                                  <  2:1663801/139‑1 (MQ=255)
                tttCAACTTTTTGGTTCAGTTCGTCAGGCATTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGACAGCATTCTGACCTCGc                                                                                                                  <  1:1120026/139‑1 (MQ=255)
                tttCAACTTTTTGGTTCAGTTCGTCAGGCATTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGACAGCATTCTGACCTCGc                                                                                                                  <  2:728336/139‑1 (MQ=255)
                  tCAACTTTTTGGTTCAGTTCGTCAGGCATTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGACAGCATTCTGACCTCGCAg                                                                                                                >  1:1247880/1‑139 (MQ=255)
                       tttttGGTTCAGTTCGTCAGGCATTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTta                                                                                                           <  2:420131/139‑1 (MQ=255)
                          ttGGTTCAGTTCGTCAGGCATTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATaa                                                                                                        <  2:536742/139‑1 (MQ=255)
                           tGGTTCAGTTCGTCAGGCATTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGa                                                                                                               >  1:1367434/1‑131 (MQ=255)
                           tGGTTCAGTTCGTCAGGCATTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGa                                                                                                               <  2:1367434/131‑1 (MQ=255)
                                   ttCGTCAGGCATTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGa                                                                                               >  1:1089683/1‑139 (MQ=255)
                                    tCGTCAGGCATTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGACAGCATTCTGACCTCGCAGATTTATAATCACATTGAg                                                                                              <  2:637276/139‑1 (MQ=255)
                                       tCAGGCATTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAg                                                                                           >  1:624176/1‑139 (MQ=255)
                                            cATTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTgc                                                                                      <  2:1089683/139‑1 (MQ=255)
                                             aTTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTgcg                                                                                     >  2:816227/1‑139 (MQ=255)
                                             aTTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGACAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTgcg                                                                                     <  2:1247880/139‑1 (MQ=255)
                                               tGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGAtt                                                                                                             >  1:722254/1‑113 (MQ=255)
                                               tGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGAtt                                                                                                             <  2:722254/113‑1 (MQ=255)
                                                gCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCa                                                                                  <  2:624176/139‑1 (MQ=255)
                                                  aaCTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGACAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCaaa                                                                                <  1:484244/139‑1 (MQ=255)
                                                       gTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATAc                                                                           >  1:135292/1‑139 (MQ=255)
                                                       gTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATAc                                                                           <  1:833409/139‑1 (MQ=255)
                                                         aCAAAATGAATTGCCTGGTTCAGGTTTTAAACCCGCGAGCAGTGCACGTGTATTTCATACTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACat                                                                         >  2:1728303/1‑139 (MQ=255)
                                                             aaTGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACatgtta                                                                     >  1:1618736/1‑139 (MQ=255)
                                                                           ttCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGcc                                                                                   <  1:1477615/111‑1 (MQ=255)
                                                                           ttCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGcc                                                                                   >  2:1477615/1‑111 (MQ=255)
                                                                           ttCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGACAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCa                                                       >  1:212494/1‑139 (MQ=255)
                                                                              aGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTa                                                    <  2:242621/139‑1 (MQ=255)
                                                                                gTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATaa                                                                                                        >  1:750493/1‑85 (MQ=255)
                                                                                gTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATaa                                                                                                        <  2:750493/85‑1 (MQ=255)
                                                                                gTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGcc                                                                                   <  1:1070096/106‑1 (MQ=255)
                                                                                gTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGcc                                                                                   >  2:1070096/1‑106 (MQ=255)
                                                                                gTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTa                                                    <  2:1264653/137‑1 (MQ=255)
                                                                                gTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTa                                                    >  1:1264653/1‑137 (MQ=255)
                                                                                    tGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCg                                                                                                                   <  1:439245/70‑1 (MQ=255)
                                                                                    tGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCg                                                                                                                   >  2:439245/1‑70 (MQ=255)
                                                                                    tGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTa                                                    >  1:964277/1‑133 (MQ=255)
                                                                                    tGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTa                                                    <  2:964277/133‑1 (MQ=255)
                                                                                     gAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGAtt                                                                                        >  1:576758/1‑96 (MQ=255)
                                                                                     gAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGAtt                                                                                        <  2:576758/96‑1 (MQ=255)
                                                                                           gcgAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTAAATCAGAACACTg                                       <  2:1624907/139‑1 (MQ=255)
                                                                                            cgAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGACAGCATTCTGACCTCGCAGAtt                                                                                                             <  1:16137/68‑1 (MQ=255)
                                                                                            cgAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGACAGCATTCTGACCTCGCAGAtt                                                                                                             >  2:16137/1‑68 (MQ=255)
                                                                                                      aCGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTAAATCAGAACACTGAACATCAGCTg                            <  2:1003513/139‑1 (MQ=255)
                                                                                                           tATTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTAAATCAGAACACTGAACATCAGCTGCGTTa                       <  1:1601350/139‑1 (MQ=255)
                                                                                                            aTTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATAc                                                                           <  1:1064403/86‑1 (MQ=255)
                                                                                                            aTTTCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATAc                                                                           >  2:1064403/1‑86 (MQ=255)
                                                                                                             tttCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTAAATCAGAACACTGAACATCAGCTGCGTTATc                     >  2:190340/1‑139 (MQ=255)
                                                                                                               tCATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTAAATCAGAACACTGAACATCAGCTGCGTTATCAg                   >  2:1462717/1‑139 (MQ=255)
                                                                                                               aacgggctgtgTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCaaa                                                                                <  1:1727267/71‑1 (MQ=255)
                                                                                                               aacgggctgtgTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCaaa                                                                                >  2:1727267/8‑78 (MQ=255)
                                                                                                                cATCTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTAAATCAGAACACTGAACATCAGCTGCGTTATCAgg                  >  2:366294/1‑139 (MQ=255)
                                                                                                                   cTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACatg                                                                        <  1:445291/82‑1 (MQ=255)
                                                                                                                   cTTTGTGTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACatg                                                                        >  2:445291/1‑82 (MQ=255)
                                                                                                                    tttGTGTTTGCAGCCCGTTAGACAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTAAATCAGAACACTGAACATCAGCTGCGTTATCAGGAAAc              <  2:212494/139‑1 (MQ=255)
                                                                                                                      tgtgTTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTAAATCAGAACACTGAACATCAGCTGCGTTATCAGGAAACgg            >  1:1223961/1‑139 (MQ=255)
                                                                                                                        tgtTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTc                                                           >  1:2452/1‑90 (MQ=255)
                                                                                                                        tgtTTGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTc                                                           <  2:2452/90‑1 (MQ=255)
                                                                                                                           ttGCAGCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTAAATCAGAACACTGAACATCAGCTGCGTTATCAGGAAACGGAGttt       <  1:366294/139‑1 (MQ=255)
                                                                                                                                gCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTAAATCAGAACACTGAACATCAGCTGCGTTATCAGGAAACGGAGTTTGTGAt  <  1:816227/139‑1 (MQ=255)
                                                                                                                                gCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTAAATCAGAACACTGAACATCAGCTGCGTTATCAGGAAACGGAGTTTGTGAt  <  2:1618736/139‑1 (MQ=255)
                                                                                                                                gCCCGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTAAATCACAACACTGAACATCAGCTGCGTTATCAGGAAACGGAGTTTGTGAt  <  2:1223961/139‑1 (MQ=255)
                                                                                                                                 cccGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTAAATCa                                               <  2:1236357/93‑1 (MQ=255)
                                                                                                                                 cccGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTAAATCa                                               >  1:1236357/1‑93 (MQ=255)
                                                                                                                                  ccGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTc                                                           >  1:791766/1‑80 (MQ=255)
                                                                                                                                  ccGTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTc                                                           <  2:791766/80‑1 (MQ=255)
                                                                                                                                    gTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACatgtt                                                                      >  2:1582861/1‑67 (MQ=255)
                                                                                                                                    gTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACatgtt                                                                      <  1:1582861/67‑1 (MQ=255)
                                                                                                                                    gTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTa                                                    <  1:1381315/85‑1 (MQ=255)
                                                                                                                                    gTTAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTa                                                    >  2:1381315/1‑85 (MQ=255)
                                                                                                                                     ttAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTAAATCAGAACACTGAACATCAGCTGCGTTATCAGGAAACGGAGTTTGTGAt  >  2:1042190/1‑134 (MQ=255)
                                                                                                                                     ttAGGCAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTAAATCAGAACACTGAACATCAGCTGCGTTATCAGGAAACGGAGTTTGTGAt  <  1:1042190/134‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                 
ACCGACTGCTCGCGCATTTCAACTTTTTGGTTCAGTTCGTCAGGCATTGCAACTAGTACAAAATGAATTGCCTGGTTCAGGTTTTGAACCCGCGAGCAGTGAACGTGTATTTCATCTTTGTGTTTGCAGCCCGTTAGACAGCATTCTGACCTCGCAGATTTATAATCACATTGAGCAGATTGCGCCAAATATACATGTTATGTTCAAGTCTTCATTAAATCAGAACACTGAACATCAGCTGCGTTATCAGGAAACGGAGTTTGTGAT  >  NC_000913/84598‑84864

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