Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,530,657 C→T G365D (GGT→GAT)  engA ← ribosome biogenesis GTPase Der

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,530,6570CT98.9% 243.5 / ‑5.2 89G365D (GGT→GAT) engAribosome biogenesis GTPase Der
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (1/87);  total (1/88)
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.75e-01

ACAACAAATGTCGGTGGTTCTATAGCAACTTGTGTTGCATAAAAGACATTCAAACGTCTACCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCACCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTAATCTTGTGC  >  NC_007793/1530561‑1530753
                                                                                                |                                                                                                
acaacaAATGTCGGTGGTTCTATAGCAACTTGTGTTGCATAAAAGACATTCAAACGTCTACCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAg                                                                                               <  1:3182139/100‑1 (MQ=255)
   acaaATGTCGGTGGTTCTATAGCAACTTGTGTTGCATAAAAGACATTCAAACGTCTACCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGtct                                                                                            <  1:4372705/100‑3 (MQ=255)
    caaATGTCGGTGGTTCTATAGCAACTTGTGTTGCATAAAAGACATTCAAACGTCTACCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCACCAGTaac                                                                                           <  1:1506693/100‑1 (MQ=255)
     aaaTGTCGGTGGTTCTATAGCAACTTGTGTTGCATAAAAGACATTCAAACGTCTACCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTaacg                                                                                          <  1:3503415/100‑2 (MQ=255)
     aaaGGTCGGTGGTTCTATAGCAACTTGTGTTGCATAAAAGACATTCAAACGTCTACCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTaaca                                                                                          <  1:5028568/100‑1 (MQ=255)
                             ttGTGTTGCATAAAAGACATCCAAACGTCTACCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGa                                                                  <  1:5674973/100‑1 (MQ=255)
                                    gCATAAAAGACATTCAAACGTCTACCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGta                                                           <  1:815414/100‑2 (MQ=255)
                                              cATTCAAACGTCTACCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGtt                                                 <  1:3193824/100‑1 (MQ=255)
                                               ctcttccgatCTACCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGttt                                                <  1:1351465/91‑1 (MQ=255)
                                               ctcttccgatCTACCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGttt                                                <  1:3738897/91‑1 (MQ=255)
                                               aTTCAAACGTCTGCCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGttt                                                <  1:5125891/100‑1 (MQ=255)
                                                ttCAAACGTCTACCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGttta                                               <  1:2438888/100‑2 (MQ=255)
                                                       gTCTACCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGc                                        <  1:4333363/100‑1 (MQ=255)
                                                        tCTACCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTgag                                       <  1:4738364/100‑3 (MQ=255)
                                                        tCTACCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCt                                       <  1:4396482/100‑1 (MQ=255)
                                                        gctcttccgatCTGTTGGGGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCt                                       <  1:4958348/90‑1 (MQ=255)
                                                         ctcttccgatCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCtt                                      <  1:380001/91‑1 (MQ=255)
                                                         ctcttccgatCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCtt                                      <  1:1454560/91‑1 (MQ=255)
                                                         ctcttccgatCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCtg                                      <  1:3621657/91‑2 (MQ=255)
                                                         ctcttccgatCTGTTGGAGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCtt                                      <  1:1754301/91‑1 (MQ=255)
                                                         cTACCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCtt                                      <  1:2941675/100‑1 (MQ=255)
                                                         cTACCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCta                                      <  1:1332171/100‑2 (MQ=255)
                                                          tcttccgatCTGTTGGCGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCttg                                     <  1:5211744/92‑2 (MQ=255)
                                                          tcttccgatCTGTCGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTc                                     <  1:3873059/92‑1 (MQ=255)
                                                          tACCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCttg                                     <  1:2964658/100‑2 (MQ=255)
                                                          tACCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTc                                     <  1:4599901/100‑1 (MQ=255)
                                                          tACCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTc                                     <  1:5675529/100‑1 (MQ=255)
                                                          tACCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACGTCATTTAGAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTc                                     <  1:3610905/100‑1 (MQ=255)
                                                           cttccgatCTGTTGGGGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCa                                    <  1:1291167/93‑1 (MQ=255)
                                                            ttccgatCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTTAt                                   <  1:3590157/94‑1 (MQ=255)
                                                            ttccgatCTGTTGGAGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCag                                   <  1:5701303/94‑2 (MQ=255)
                                                             cTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCAtt                                  <  1:2152095/100‑1 (MQ=255)
                                                             cTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCAtt                                  <  1:4914944/100‑1 (MQ=255)
                                                              tttGTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTg                                 <  1:3397577/100‑1 (MQ=255)
                                                              cgatctgtttGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTg                                 <  1:3486340/92‑1 (MQ=255)
                                                              ccgatCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACCCGTTTTTTATGGTTTTCACTTGCTTCATTg                                 <  1:4140469/96‑1 (MQ=255)
                                                              ccgatCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTg                                 <  1:4696845/96‑1 (MQ=255)
                                                              ccgatCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTg                                 <  1:2367756/96‑1 (MQ=255)
                                                              ccgatCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAATTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTg                                 <  1:2435228/96‑1 (MQ=255)
                                                              ccgatCTGGTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTg                                 <  1:3205345/96‑1 (MQ=255)
                                                                 gTCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGAta                              <  1:3362919/100‑2 (MQ=255)
                                                                 atCTGTTGGCGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATg                              <  1:1307115/99‑1 (MQ=255)
                                                                  tCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGt                             <  1:3354073/100‑1 (MQ=255)
                                                                  tCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGt                             <  1:5601038/100‑1 (MQ=255)
                                                                  tCTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGt                             <  1:898227/100‑1 (MQ=255)
                                                                  tCTGTTGGAGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGt                             <  1:1256605/100‑1 (MQ=255)
                                                                   cTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTa                            <  1:112482/100‑1 (MQ=255)
                                                                   cTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTa                            <  1:622499/100‑1 (MQ=255)
                                                                   cTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTa                            <  1:3090285/100‑1 (MQ=255)
                                                                   cTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTa                            <  1:1913500/100‑1 (MQ=255)
                                                                   cTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTa                            <  1:5082996/100‑1 (MQ=255)
                                                                   cTGTTGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTa                            >  1:3011450/1‑100 (MQ=255)
                                                                   cTGTTGGGGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTa                            <  1:1743993/100‑1 (MQ=255)
                                                                   cTGTTGGCGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTa                            <  1:5301529/100‑1 (MQ=255)
                                                                   cTGTTGGCGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGCTTTCACTTGCTTCATTGATGTa                            <  1:4496994/100‑1 (MQ=255)
                                                                   cTGTTGCTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTa                            <  1:1668555/100‑1 (MQ=255)
                                                                    ctgtGGTGTAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTag                           <  1:454723/97‑2 (MQ=255)
                                                                         gtgtAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAggggg                      <  1:3774969/100‑3 (MQ=255)
                                                                           gtAGGGTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGatac                    <  1:1308775/100‑4 (MQ=255)
                                                                                gTTCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTa               <  1:3014320/100‑1 (MQ=255)
                                                                                  tCATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACg             <  1:1793836/100‑1 (MQ=255)
                                                                                   cATGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGt            <  1:94567/100‑1 (MQ=255)
                                                                                    aTGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTa           <  1:3906149/100‑1 (MQ=255)
                                                                                    aTGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTa           <  1:4025535/100‑1 (MQ=255)
                                                                                    aTGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTa           <  1:2313428/100‑1 (MQ=255)
                                                                                    aTGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTa           <  1:2043781/100‑1 (MQ=255)
                                                                                    aTGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTa           <  1:3853030/100‑1 (MQ=255)
                                                                                    aTGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGGTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTa           <  1:2151831/100‑1 (MQ=255)
                                                                                    aTGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAAATTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTa           <  1:693780/100‑1 (MQ=255)
                                                                                     tGGAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTag          <  1:2242587/100‑2 (MQ=255)
                                                                                       gAAATTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTAAta        <  1:527338/100‑2 (MQ=255)
                                                                                        aaaTTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTAATCt       <  1:4855208/100‑1 (MQ=255)
                                                                                          aTTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTAATCtta     <  1:4723011/100‑2 (MQ=255)
                                                                                          aTTGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTAATCTtg     <  1:3996598/100‑1 (MQ=255)
                                                                                           ttGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTAATCTtgt    <  1:542086/100‑1 (MQ=255)
                                                                                            tGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTAATCTtgtg   <  1:1591620/100‑1 (MQ=255)
                                                                                            tGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTAATCTtgtg   <  1:78624/100‑1 (MQ=255)
                                                                                            tGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTAATCTtgtg   <  1:3388160/100‑1 (MQ=255)
                                                                                            tGCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTAATCTtgtg   <  1:2684930/100‑1 (MQ=255)
                                                                                             gCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTAATCtagcac  <  1:3822079/100‑6 (MQ=255)
                                                                                             gCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTAATCTTGTGc  <  1:5479590/100‑1 (MQ=255)
                                                                                             gCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTAATCTTGTGc  <  1:5443729/100‑1 (MQ=255)
                                                                                             gCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTAATCTTGTGc  <  1:3869680/100‑1 (MQ=255)
                                                                                             gCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTAATCTTGTGc  <  1:4100314/100‑1 (MQ=255)
                                                                                             gCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTAATCTTGTGc  <  1:195264/100‑1 (MQ=255)
                                                                                             gCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTAATCTTGTGc  <  1:3055706/100‑1 (MQ=255)
                                                                                             gCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTAATCTTGTGc  <  1:2927170/100‑1 (MQ=255)
                                                                                             gCATCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTAATCTTGTGc  <  1:3674860/100‑1 (MQ=255)
                                                                                             gCATCAGTAACAACTTCATTTAAAGTTGAACTTTAAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTAATCTTGTGc  <  1:2403778/100‑1 (MQ=255)
                                                                                                |                                                                                                
ACAACAAATGTCGGTGGTTCTATAGCAACTTGTGTTGCATAAAAGACATTCAAACGTCTACCTTTGTCTGTTGGTGTAGGGTTCATGGAAATTGCACCAGTAACAACTTCATTTAAAGTTGAACTTTGAACACGTTTTTTATGGTTTTCACTTGCTTCATTGATGTAAGGGAATAATGTACGTAATCTTGTGC  >  NC_007793/1530561‑1530753

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