Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,862,488 G→A A471A (GCC→GCT yidE ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,4880GA99.7% 1311.0 / ‑5.6 361A471A (GCC→GCTyidEputative transport protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (180/180);  total (180/181)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

Only 100 of 368 total aligned reads displayed.

TGCCAGCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGGGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCGAGAGGA  >  NC_000913/3862398‑3862576
                                                                                          |                                                                                        
tGCCAGCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGc                                                                                        <  2:7702062/93‑1 (MQ=255)
 gCCAGCATTCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCa                                                                                       <  1:1215918/93‑1 (MQ=255)
 gCCAGCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCa                                                                                       <  1:4384797/93‑1 (MQ=255)
  ccAGCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCAc                                                                                      <  1:8075751/93‑1 (MQ=255)
   cAGCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCAcc                                                                                     <  2:1802839/93‑1 (MQ=255)
      cATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATa                                                                                  <  1:5673790/93‑1 (MQ=255)
           cGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCCTAACCaa                                                                             >  2:931802/1‑93 (MQ=255)
            gCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAAt                                                                            <  2:1517638/93‑1 (MQ=255)
            gCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAAt                                                                            <  1:3039186/93‑1 (MQ=255)
              acaTGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATcc                                                                          <  2:10032720/93‑1 (MQ=255)
              acaTGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATcc                                                                          >  1:7601414/1‑93 (MQ=255)
              acaTGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATcc                                                                          >  1:6528684/1‑93 (MQ=255)
              acaTGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATcc                                                                          <  2:5318950/93‑1 (MQ=255)
               caTGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCa                                                                         <  1:7734809/93‑1 (MQ=255)
                 aggTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGGGCACCATAACCAATCCAGc                                                                       <  2:5848515/92‑1 (MQ=255)
                    tCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTa                                                                    <  1:467852/93‑1 (MQ=255)
                    gcAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTa                                                                    <  1:5490620/92‑1 (MQ=255)
                      aGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGc                                                                  >  1:606776/1‑93 (MQ=255)
                       ggTAATTCATTTTTGCTAACAACCCCCCCCGAAAGCCAACAGTAATCAGCGGGACGGCGGGGGTCCAGGCACCATAAACAAACCACCTTAGcc                                                                 >  2:8763942/1‑93 (MQ=255)
                       ggTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGcc                                                                 >  2:5596085/1‑93 (MQ=255)
                       ggTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGcc                                                                 <  1:2129949/93‑1 (MQ=255)
                       ggTAATTCATGTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGcc                                                                 <  1:9035007/93‑1 (MQ=255)
                        gTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGccc                                                                >  1:8584975/1‑93 (MQ=255)
                        gTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGccc                                                                <  2:4967562/93‑1 (MQ=255)
                         tAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGTCCt                                                               >  2:2219214/1‑93 (MQ=255)
                         tAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCt                                                               <  1:3586293/93‑1 (MQ=255)
                         tAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCt                                                               >  1:7183232/1‑93 (MQ=255)
                          aaTTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCtt                                                              <  1:8389095/93‑1 (MQ=255)
                           aTTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTc                                                             >  2:3227226/1‑93 (MQ=255)
                           aTTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTc                                                             >  2:4375126/1‑93 (MQ=255)
                           aTTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTc                                                             <  2:513123/93‑1 (MQ=255)
                             tCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGc                                                           >  2:6197553/1‑93 (MQ=255)
                               aTTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCa                                                         <  2:121675/93‑1 (MQ=255)
                               aTTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCa                                                         >  1:3854487/1‑93 (MQ=255)
                               aTTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCa                                                         <  1:9288732/93‑1 (MQ=255)
                                ttttGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCAt                                                        <  1:5891770/93‑1 (MQ=255)
                                ttttGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCAt                                                        <  2:643970/93‑1 (MQ=255)
                                 tttGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCAtt                                                       <  2:5411899/93‑1 (MQ=255)
                                  ttGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTg                                                      <  1:135251/93‑1 (MQ=255)
                                  ttGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTg                                                      >  1:4201247/1‑93 (MQ=255)
                                   tGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGa                                                     >  2:2522266/1‑93 (MQ=255)
                                   tGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGa                                                     >  1:8676214/1‑93 (MQ=255)
                                    ggCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGAc                                                    <  1:6757321/93‑1 (MQ=255)
                                    ggCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGAc                                                    >  2:2619901/1‑93 (MQ=255)
                                    ggCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGAc                                                    <  2:2122817/93‑1 (MQ=255)
                                      cTAACAGCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCa                                                  <  2:7546723/93‑1 (MQ=255)
                                       tAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAg                                                 >  2:8264422/1‑93 (MQ=255)
                                         aCATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAggg                                               <  1:9442942/93‑1 (MQ=255)
                                         aCATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAggg                                               <  1:2818760/93‑1 (MQ=255)
                                          cATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCCCCATTGACCAggga                                              >  2:9866411/1‑92 (MQ=255)
                                           aTCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTa                                             >  2:4331572/1‑93 (MQ=255)
                                           aTCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTa                                             >  2:3972245/1‑93 (MQ=255)
                                           aTCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTa                                             >  1:4907634/1‑93 (MQ=255)
                                            tCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTAt                                            <  2:1995121/93‑1 (MQ=255)
                                            tCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGGGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTAt                                            >  1:1147216/1‑93 (MQ=255)
                                             ccGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTAtt                                           >  1:9397813/1‑93 (MQ=255)
                                             ccGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTAtt                                           <  1:7372299/93‑1 (MQ=255)
                                                  ccAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACaa                                      <  2:9104707/93‑1 (MQ=255)
                                                    aGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACaaaa                                    >  2:3727056/1‑93 (MQ=255)
                                                    aGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACaaaa                                    >  2:2930523/1‑93 (MQ=255)
                                                    aGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACaaaa                                    <  2:3568065/93‑1 (MQ=255)
                                                     gAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCGGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAAt                                   >  2:6565259/1‑93 (MQ=255)
                                                     gAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAAt                                   >  1:2839064/1‑93 (MQ=255)
                                                     gAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAAt                                   <  1:716752/93‑1 (MQ=255)
                                                     gAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAAt                                   >  1:2513473/1‑93 (MQ=255)
                                                     gAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAAt                                   >  2:9568624/1‑93 (MQ=255)
                                                       aTGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACCAAATcc                                 >  2:3156703/1‑93 (MQ=255)
                                                       aTGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATcc                                 >  1:2271788/1‑93 (MQ=255)
                                                        tGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATccc                                >  2:137320/1‑93 (MQ=255)
                                                        tGCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATccc                                >  1:5985209/1‑93 (MQ=255)
                                                         gCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATcccc                               >  2:7814014/1‑93 (MQ=255)
                                                         gCCAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATcccc                               >  1:7520544/1‑93 (MQ=255)
                                                          ccAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCcca                              >  2:4912290/1‑93 (MQ=255)
                                                           cAACAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCccac                             >  1:1588003/1‑93 (MQ=255)
                                                            aaCAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCccacc                            >  1:1883811/1‑93 (MQ=255)
                                                            aaCAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCccacc                            <  2:7263907/93‑1 (MQ=255)
                                                             aCAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCccacca                           >  2:8441825/1‑93 (MQ=255)
                                                              cAGTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGAGTATTCACAAAATCCCCACCAg                          >  1:6876825/1‑93 (MQ=255)
                                                                gTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCCATCCAGCTTAGCCCTTAGCCATTGACCAGGGTATTCAAAACAGCCCCACCAGAt                        >  1:1024398/1‑93 (MQ=255)
                                                                gTAATCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGAt                        >  2:7404301/1‑93 (MQ=255)
                                                                   aTCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTc                     <  2:6693125/93‑1 (MQ=255)
                                                                    tCAGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCa                    <  1:3227226/93‑1 (MQ=255)
                                                                      aGCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGa                  >  1:3441387/1‑93 (MQ=255)
                                                                       gCGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGAc                 >  2:328088/1‑93 (MQ=255)
                                                                        cGGAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGAcc                <  2:7076250/93‑1 (MQ=255)
                                                                          gAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCaa              <  2:1658433/93‑1 (MQ=255)
                                                                          gAACGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCaa              <  1:5145672/93‑1 (MQ=255)
                                                                           aaCGGCGGTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAAc             <  1:785933/93‑1 (MQ=255)
                                                                               gcggTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGAcc         >  2:8160483/1‑93 (MQ=255)
                                                                                 ggTGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCga       <  1:4028711/93‑1 (MQ=255)
                                                                                   tGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCgaga     >  2:680465/1‑93 (MQ=255)
                                                                                   tGATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCgaga     <  1:8293865/93‑1 (MQ=255)
                                                                                    gATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCACCGACCGagag    >  2:2602612/1‑93 (MQ=255)
                                                                                    gATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCGagag    <  2:6099916/93‑1 (MQ=255)
                                                                                    gATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCGagag    <  2:4652687/93‑1 (MQ=255)
                                                                                    gATCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCGagag    <  1:3721199/93‑1 (MQ=255)
                                                                                     aTCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCGAGAgg   <  2:6819576/93‑1 (MQ=255)
                                                                                     aTCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCGAGAgg   >  1:8894483/1‑93 (MQ=255)
                                                                                     aTCAGAGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCGAGAgg   >  1:9759718/1‑93 (MQ=255)
                                                                                      tCAGAGCACCATAACCAATCCAGATTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCGAGAGGa  <  2:5546803/93‑1 (MQ=255)
                                                                                          |                                                                                        
TGCCAGCATCCCGCACATGGTCAGGTAATTCATTTTGGCTAACATCCGCGCCAGAATGCCAACAGTAATCAGCGGAACGGCGGTGATCAGGGCACCATAACCAATCCAGCTTAGCCCTTCGCCATTGACCAGGGTATTCACAAAATCCCCACCAGATTTCAGACCAACGACCGAGAGGA  >  NC_000913/3862398‑3862576

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: