Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 660,616 C→T A255T (GCC→ACC)  ybeF ← LysR family putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913660,6160CT98.9% 306.8 / ‑4.5 88A255T (GCC→ACC) ybeFLysR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (43/44);  total (43/45)
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.87e-01

CATTTAAGGACGGAAAAGGCAGTTTTTCCAGCGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGCAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTGAACTCTCTGCCGCAGACCGCTGAAAT  >  NC_000913/660520‑660711
                                                                                                |                                                                                               
cATTTAAGGACGGAAAAGGCAGTTTTTCCAGCGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAAc                                                                                              >  2:292559/1‑100 (MQ=255)
 aTTTAAGGACGGAAAAGGCAGTTTTTCCAGCGGCCAGCAGCGGCTAAACAGGTTATAAGAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAAcc                                                                                             <  2:2421480/100‑1 (MQ=255)
     aaGGACGGAAAAGGCAGTTTTTCCAGCGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTg                                                                                         >  1:714768/1‑100 (MQ=255)
     aaGGACGGAAAAGGCAGTTTTTCCAGCGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTg                                                                                         >  2:1373158/1‑100 (MQ=255)
      aGGACGGAAAAGGCAGTTTTTCCAGCGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGc                                                                                        <  1:2424667/100‑1 (MQ=255)
          cGGAAAAGGCAGTTTTTCCAGCGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCg                                                                                    <  2:2571936/100‑1 (MQ=255)
          cGGAAAAGGCAGTTTTTCCAGCGGCCAGCAGCGGCTAAAAAGAATATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCg                                                                                    <  2:123336/100‑1 (MQ=255)
            gAAAAGGCAGTTTTTCCAGCGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGAt                                                                                  >  2:1405162/1‑100 (MQ=255)
             aaaaGGCAGTTTTTCCAGCGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGAtt                                                                                 <  2:333103/100‑1 (MQ=255)
              aaaGGCAGTTTTTCCAGCGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTg                                                                                <  1:399731/100‑1 (MQ=255)
              aaaGGCAGTTTTTCCAGCGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTg                                                                                <  1:1754648/100‑1 (MQ=255)
                aGGCAGTTTTTCCAGCGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTc                                                                              >  2:1766883/1‑100 (MQ=255)
                 ggCAGTTTTTCCAGCGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCa                                                                             >  2:2257221/1‑100 (MQ=255)
                 ggCAGTTTTTCCAGCGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCa                                                                             >  1:1328105/1‑100 (MQ=255)
                 ggCAGTTTTTCCAGCGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCa                                                                             >  2:2271780/1‑100 (MQ=255)
                     gTTTTTCCAGCGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAAt                                                                         >  2:129535/1‑100 (MQ=255)
                     gTTTTTCCAGCGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAAt                                                                         >  2:1121092/1‑100 (MQ=255)
                        tttCCAGCGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAAtgtt                                                                      <  1:391040/100‑1 (MQ=255)
                         ttCCAGCGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAAtgttg                                                                     <  2:34277/100‑1 (MQ=255)
                           ccAGCGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTa                                                                   >  2:2877404/1‑100 (MQ=255)
                            cAGCGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAg                                                                  >  2:2537259/1‑100 (MQ=255)
                              gCGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCt                                                                >  2:2137411/1‑100 (MQ=255)
                               cGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTg                                                               <  2:1800297/100‑1 (MQ=255)
                                ggCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTgg                                                              <  1:831406/100‑1 (MQ=255)
                                 gCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGt                                                             <  2:1302146/100‑1 (MQ=255)
                                  ccAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTg                                                               >  1:1780248/1‑97 (MQ=255)
                                  ccAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTg                                                            >  2:2267627/1‑100 (MQ=255)
                                     gcagcGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGaaa                                                         <  2:1260588/100‑1 (MQ=255)
                                       agcGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAAtt                                                       >  1:2894275/1‑100 (MQ=255)
                                       agcGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAAtt                                                       <  2:2912410/100‑1 (MQ=255)
                                       agcGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAAtt                                                       <  2:647921/100‑1 (MQ=255)
                                          ggCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGAt                                                    <  1:313615/100‑1 (MQ=255)
                                          ggCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGAt                                                    >  2:685365/1‑100 (MQ=255)
                                          ggCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGAt                                                    <  2:2885667/100‑1 (MQ=255)
                                           tcTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGAtt                                                   <  1:2953528/99‑1 (MQ=255)
                                           gCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGAtt                                                   <  2:409904/100‑1 (MQ=255)
                                            cTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGAttt                                                  >  2:405954/1‑100 (MQ=255)
                                             tAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTg                                                 >  2:1069252/1‑100 (MQ=255)
                                             tAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTg                                                 >  2:1485842/1‑100 (MQ=255)
                                             tAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTg                                                 >  1:1605774/1‑100 (MQ=255)
                                              ccgatctttATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGc                                                <  2:1210528/93‑1 (MQ=255)
                                              aaaCAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGAttaga                                                >  2:2953512/1‑97 (MQ=255)
                                              aaaCAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGc                                                >  1:846693/1‑100 (MQ=255)
                                              aaaCAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGc                                                >  2:1685149/1‑100 (MQ=255)
                                               aaCAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGcc                                               <  1:2877420/100‑1 (MQ=255)
                                               aaCAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGcc                                               <  2:1701018/100‑1 (MQ=255)
                                               aaCAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGcc                                               <  2:1459834/100‑1 (MQ=255)
                                               aaCAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGcc                                               <  2:2894259/100‑1 (MQ=255)
                                               aaCAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGCAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGcc                                               <  2:2013231/100‑1 (MQ=255)
                                                aCAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTATCTGGTGAAATTGATTTGcc                                               >  1:2410472/1‑99 (MQ=255)
                                                 cAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCgg                                             <  2:1306158/100‑1 (MQ=255)
                                                 cAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCgg                                             >  2:1340295/1‑100 (MQ=255)
                                                   ggTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGCTTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTc                                           >  2:1283723/1‑100 (MQ=255)
                                                   ggTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTc                                           >  2:1274110/1‑100 (MQ=255)
                                                   ggTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTc                                           >  2:1988111/1‑100 (MQ=255)
                                                   ggTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTc                                           >  1:698933/1‑100 (MQ=255)
                                                    gTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTcc                                          >  2:1248168/1‑100 (MQ=255)
                                                     ttATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGcagatcgg                                         >  1:1210534/1‑93 (MQ=255)
                                                      tataAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCgc                                        <  2:721554/100‑2 (MQ=255)
                                                       ataAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGa                                       <  1:733098/100‑1 (MQ=255)
                                                        taAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGaa                                      <  1:1485851/100‑1 (MQ=255)
                                                         aaaaaCGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGaaa                                     >  2:1362011/1‑100 (MQ=255)
                                                         aaaaaCGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGaaa                                     <  1:196402/100‑1 (MQ=255)
                                                          aaaaCGGCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAAc                                    <  1:1401291/100‑1 (MQ=255)
                                                                gCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCt                              >  1:1580528/1‑100 (MQ=255)
                                                                gCTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCt                              <  1:2137423/100‑1 (MQ=255)
                                                                 cTGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCtt                             <  1:1638992/100‑1 (MQ=255)
                                                                  tGGGAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTg                            <  1:662362/100‑1 (MQ=255)
                                                                   gggAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTGa                           <  2:868585/100‑1 (MQ=255)
                                                                   gggAATAATCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTGa                           <  2:846690/100‑1 (MQ=255)
                                                                       ataatCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTGAActc                       <  2:1107835/100‑1 (MQ=255)
                                                                         aatCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTGAActctc                     <  2:2002334/100‑1 (MQ=255)
                                                                           tCGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTGAACTCTCTg                   <  1:129537/100‑1 (MQ=255)
                                                                            cGCTAACATGTCACTGTTTGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTGAACTCTCTg                   <  1:299472/99‑1 (MQ=255)
                                                                            cGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTGAACTCTCTGc                  >  2:2639650/1‑100 (MQ=255)
                                                                            cGCTAACATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTGAACTCTCTGc                  <  1:2484130/100‑1 (MQ=255)
                                                                                aaCATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTGAACTCTCTGCCGCa              <  1:2670088/100‑1 (MQ=255)
                                                                                 aCATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTGAACTCTCTGCCGCAg             <  1:1248174/100‑1 (MQ=255)
                                                                                  cATGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTGAACTCTCTGCCGCAGa            >  1:2970594/1‑100 (MQ=255)
                                                                                    tGTCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTGAACTCTCTGCCGCAGAcc          <  1:292561/100‑1 (MQ=255)
                                                                                      tCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTGAACTCTCTGCCGCAGACCGc        >  1:1494015/1‑100 (MQ=255)
                                                                                      tCACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTGAACTCTCTGCCGCAGACCGc        >  2:1438108/1‑100 (MQ=255)
                                                                                       cACTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTGAACTCTCTGCCGCAGACCGCt       >  2:2148153/1‑100 (MQ=255)
                                                                                         cTGTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTGAACTCTCTGCCGCAGACCGCTGa     >  1:1922707/1‑100 (MQ=255)
                                                                                           gTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTGAACTCTCTGCCGCAGACCGCTGaaa   >  1:2224427/1‑100 (MQ=255)
                                                                                           gTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTGAACTCTCTGCCGCAGACCGCTGaaa   >  1:1266105/1‑100 (MQ=255)
                                                                                           gTTGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTGAACTCTCTGCCGCAGACCGCTGaaa   >  1:123775/1‑100 (MQ=255)
                                                                                            ttGGTAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTGAACTCTCTGCCGCAGACCGCTGAAAt  <  2:1232209/100‑1 (MQ=255)
                                                                                                |                                                                                               
CATTTAAGGACGGAAAAGGCAGTTTTTCCAGCGGCCAGCAGCGGCTAAACAGGTTATAAAAACGGCTGGGAATAATCGCTAACATGTCACTGTTGGCAACCAGTGCAGCGATTGTCAAAATGTTGTAGCTGGTGAAATTGATTTGCCGGTCCGGAAACATTTCTTGAACTCTCTGCCGCAGACCGCTGAAAT  >  NC_000913/660520‑660711

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