Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,001,421 2 bp→GC coding (381‑382/903 nt) rarD ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,001,4210CG100.0% 891.5 / NA 243G128R (GGA→CGA) rarDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (123/120);  total (123/120)
*NC_0009134,001,4220GC99.6% 889.2 / NA 243D127E (GAC→GAG) rarDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (122/120);  minor base A (1/0);  total (123/120)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 253 total aligned reads displayed.

TAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGTTCCAGCATATGGTGATTGTT  >  NC_000913/4001334‑4001508
                                                                                       ||                                                                                      
taataGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAAtcgctcgc                                                                                    >  2:4533991/1‑93 (MQ=255)
 aataGGTAGCGAACCAAAAGTCCACAGCTGGCCTAACCCGCCACATACCGCCAGAATCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGcc                                                                                   <  1:300253/93‑1 (MQ=255)
 aataGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGcc                                                                                   >  1:6139826/1‑93 (MQ=255)
   taGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGa                                                                                 <  2:6456955/93‑1 (MQ=255)
   taGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGa                                                                                 <  1:4327009/93‑1 (MQ=255)
    aGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAg                                                                                >  1:4890701/1‑93 (MQ=255)
    aGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAg                                                                                >  2:1533857/1‑93 (MQ=255)
      gTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGa                                                                              >  1:1338271/1‑93 (MQ=255)
       tAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGaa                                                                             <  1:4033848/93‑1 (MQ=255)
        aGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGaaa                                                                            <  2:3119615/93‑1 (MQ=255)
         gCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGaaaa                                                                           >  1:2749744/1‑93 (MQ=255)
         gCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGaaaa                                                                           <  1:6576372/93‑1 (MQ=255)
             aCCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAatcat                                                                       >  1:2635145/1‑93 (MQ=255)
             aCCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCAAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAatcat                                                                       <  1:3253785/93‑1 (MQ=255)
              ccAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAatcatc                                                                      >  1:2359088/1‑93 (MQ=255)
              ccAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAatcatc                                                                      <  2:188451/93‑1 (MQ=255)
               cAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATcc                                                                     >  2:5063384/1‑93 (MQ=255)
               cAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATcc                                                                     >  1:1146184/1‑93 (MQ=255)
                aaaaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATccc                                                                    <  1:2330297/93‑1 (MQ=255)
                aaaaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATccc                                                                    >  1:2916630/1‑93 (MQ=255)
                aaaaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATccc                                                                    <  2:5624770/93‑1 (MQ=255)
                 aaaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATcccc                                                                   >  1:226984/1‑93 (MQ=255)
                  aaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCa                                                                  <  1:831703/93‑1 (MQ=255)
                  aaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCa                                                                  <  1:2374895/93‑1 (MQ=255)
                   aGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAg                                                                 <  2:3637228/93‑1 (MQ=255)
                     tCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGca                                                               >  2:1565696/1‑93 (MQ=255)
                      ccACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGcac                                                              <  1:503320/93‑1 (MQ=255)
                       cacaGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGcaca                                                             >  2:199218/1‑93 (MQ=255)
                           gCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGt                                                         >  1:5847516/1‑93 (MQ=255)
                           gCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGt                                                         >  1:3424608/1‑93 (MQ=255)
                              ggACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTcc                                                      >  2:2189664/1‑92 (MQ=255)
                              ggACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCa                                                      >  1:4113482/1‑93 (MQ=255)
                              ggACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCa                                                      >  2:1756922/1‑93 (MQ=255)
                               gACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCAc                                                     >  1:4223102/1‑93 (MQ=255)
                               gACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCAc                                                     <  1:2579999/93‑1 (MQ=255)
                               gACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCAc                                                     <  2:5298071/93‑1 (MQ=255)
                               gACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGAACAATGTTCAc                                                     >  2:553133/1‑93 (MQ=255)
                                 cTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCa                                                   >  2:2322799/1‑93 (MQ=255)
                                 cTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCa                                                   <  2:844319/93‑1 (MQ=255)
                                  tAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAg                                                  >  2:3709580/1‑93 (MQ=255)
                                  tAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAg                                                  >  2:931404/1‑93 (MQ=255)
                                   aaCACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGc                                                 >  1:1483304/1‑93 (MQ=255)
                                     cacGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCgg                                               <  2:2059102/93‑1 (MQ=255)
                                        gCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGtt                                            <  1:358751/93‑1 (MQ=255)
                                        gCCACATATCCCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGtt                                            <  1:6211817/93‑1 (MQ=255)
                                         ccACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTta                                           >  2:4846828/1‑93 (MQ=255)
                                          cccgtatCACCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCCCAATGTTCACCAGCGGGTtaa                                          <  2:3306813/89‑1 (MQ=255)
                                          cacaTATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaa                                          >  2:6510713/1‑93 (MQ=255)
                                            caTATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGataata                                        >  1:3934020/1‑93 (MQ=255)
                                            caTATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaata                                        <  2:5372007/93‑1 (MQ=255)
                                              tatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATaaa                                      <  1:3353781/93‑1 (MQ=255)
                                              tatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATaaa                                      >  1:4793710/1‑93 (MQ=255)
                                              tatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATaaa                                      >  2:5087713/1‑93 (MQ=255)
                                              tatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATaaa                                      <  1:5616165/93‑1 (MQ=255)
                                               ttCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAg                                     <  2:6009186/92‑1 (MQ=255)
                                               atCGCCAGACCCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAg                                     <  1:2402213/93‑1 (MQ=255)
                                               atCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAg                                     <  1:6374786/93‑1 (MQ=255)
                                               atCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAg                                     >  2:3374443/1‑93 (MQ=255)
                                                 cGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTa                                   <  1:6452877/93‑1 (MQ=255)
                                                 cGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTa                                   <  1:1636494/93‑1 (MQ=255)
                                                  gCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTaa                                  >  1:430759/1‑93 (MQ=255)
                                                   ccAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAAc                                 <  2:3835009/93‑1 (MQ=255)
                                                     aGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCa                               >  2:1127567/1‑93 (MQ=255)
                                                     aGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCa                               <  1:2534421/93‑1 (MQ=255)
                                                      gAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCaa                              <  2:5878184/93‑1 (MQ=255)
                                                        aTCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAgg                            >  1:4969180/1‑93 (MQ=255)
                                                         tCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGc                           <  1:3612879/93‑1 (MQ=255)
                                                         gccccGCCAGCCGTTGCATCCGGCCGAATCGCTCGCCGACGAAAATCATCCCCAGCACCATGTTCACCAGCGGGTTAAAAAAGTAACCAAGGc                           <  1:5744524/90‑1 (MQ=255)
                                                          cACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCt                          <  2:3589902/93‑1 (MQ=255)
                                                          cACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCt                          <  1:5331633/93‑1 (MQ=255)
                                                           cccGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTc                         <  1:5915817/92‑1 (MQ=255)
                                                           aCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTc                         >  1:5117838/1‑93 (MQ=255)
                                                           aCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTc                         >  2:68126/1‑93 (MQ=255)
                                                           aCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTc                         >  1:3534359/1‑93 (MQ=255)
                                                            ccGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTcg                        <  2:4151167/93‑2 (MQ=255)
                                                            ccGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTcg                        <  2:3414215/93‑2 (MQ=255)
                                                             cGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTcgc                       <  1:1013139/93‑3 (MQ=255)
                                                              gccagccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTcgct                      <  2:2997876/93‑4 (MQ=255)
                                                                 agccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACGAAGGCTCGCTTcc                   >  1:189062/1‑93 (MQ=255)
                                                                 agccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTcc                   >  2:4498137/1‑93 (MQ=255)
                                                                 agccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTcc                   >  2:2646266/1‑93 (MQ=255)
                                                                 agccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTcc                   <  1:4929158/93‑1 (MQ=255)
                                                                 agccaTCGCATCGGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTcc                   <  2:3991485/93‑1 (MQ=255)
                                                                  gccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCa                  >  1:5487255/1‑93 (MQ=255)
                                                                  gccaTTGCATCCGGCGGAATCGCTCGCCGACGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCa                  >  2:6195507/1‑93 (MQ=255)
                                                                   ccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAg                 >  2:3085821/1‑93 (MQ=255)
                                                                   ccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAg                 <  2:4280719/93‑1 (MQ=255)
                                                                    caTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGc                >  1:4543192/1‑93 (MQ=255)
                                                                     aTTGCGTCCGGCGGAATCGCTAGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCa               <  2:6241391/93‑1 (MQ=255)
                                                                     aTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCa               <  1:5530955/93‑1 (MQ=255)
                                                                        gCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCatat            <  1:1971723/93‑1 (MQ=255)
                                                                          ttcaggcgcAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAAGCTCGCTACCAGCATATgg          <  1:1688065/89‑1 (MQ=255)
                                                                          aTCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATgg          <  2:4909649/93‑1 (MQ=255)
                                                                           tCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGt         >  1:5290552/1‑93 (MQ=255)
                                                                            ccGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTg        >  1:907612/1‑93 (MQ=255)
                                                                              ggcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAt      >  1:4521923/1‑93 (MQ=255)
                                                                               gcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTGACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAtt     >  2:4676901/1‑93 (MQ=255)
                                                                               gcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAtt     >  2:2376170/1‑93 (MQ=255)
                                                                               gcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAtt     <  2:1353020/93‑1 (MQ=255)
                                                                                  gAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGtt  <  1:5173511/93‑1 (MQ=255)
                                                                                       ||                                                                                      
TAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGTTCCAGCATATGGTGATTGTT  >  NC_000913/4001334‑4001508

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