Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,806,424 G→A T86T (ACG→ACA)  kdtA → 3‑deoxy‑D‑manno‑octulosonic‑acid transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,806,4240GA99.7% 1035.5 / NA 291T86T (ACG→ACA) kdtA3‑deoxy‑D‑manno‑octulosonic‑acid transferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (132/158);  total (132/159)
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 299 total aligned reads displayed.

CTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACC  >  NC_000913/3806334‑3806512
                                                                                          |                                                                                        
ctccgtcGGGAAAACTCTGGCGGCAACCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTCCCGATTCCCGTAACAACCATGAcccc                                                                                        <  2:1642766/93‑1 (MQ=255)
 tccgtcGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCa                                                                                       >  2:4385985/1‑93 (MQ=255)
  cggtcGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCaa                                                                                      >  2:6029968/3‑93 (MQ=255)
  ccgtcGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCaa                                                                                      <  2:265356/93‑1 (MQ=255)
   cgtcGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAAc                                                                                     >  2:4126225/1‑93 (MQ=255)
    gtcGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAAcc                                                                                    <  2:905145/93‑1 (MQ=255)
    gtcGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAAcc                                                                                    <  2:5030859/93‑1 (MQ=255)
     tcGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCg                                                                                   <  1:4250612/93‑1 (MQ=255)
       ggTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGt                                                                                 <  1:4227808/93‑1 (MQ=255)
       ggTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGt                                                                                 <  1:5138616/93‑1 (MQ=255)
        gTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGtt                                                                                >  1:3755906/1‑93 (MQ=255)
         tGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTc                                                                               >  1:3823008/1‑93 (MQ=255)
         ggAAACTGTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTc                                                                               <  1:6228411/92‑1 (MQ=255)
          gAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCg                                                                              <  1:675446/93‑1 (MQ=255)
             aCTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAATCGGTTCGGAg                                                                           >  1:5555305/1‑93 (MQ=255)
                 tGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGcgcg                                                                       >  2:3592898/1‑93 (MQ=255)
                   gcggcAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTa                                                                     <  2:4218833/93‑1 (MQ=255)
                   gcggcAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTa                                                                     >  1:4352703/1‑93 (MQ=255)
                    cggcAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTAc                                                                    <  1:4541707/93‑1 (MQ=255)
                     ggcAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACa                                                                   >  1:4040982/1‑93 (MQ=255)
                     ggcAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACa                                                                   >  2:3539572/1‑93 (MQ=255)
                      gcAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACaa                                                                  <  1:2164121/93‑1 (MQ=255)
                      gcAACCCCGTTGGTGCGCGCTCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACaa                                                                  <  2:2612099/93‑1 (MQ=255)
                       cAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAAt                                                                 <  1:498366/93‑1 (MQ=255)
                       cAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAAt                                                                 <  1:1852211/93‑1 (MQ=255)
                        aaTCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATc                                                                <  2:4826572/93‑1 (MQ=255)
                          tCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCgg                                                              >  1:991388/1‑93 (MQ=255)
                          tCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCgg                                                              >  1:1080247/1‑93 (MQ=255)
                              gTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCttt                                                          <  1:2252488/93‑1 (MQ=255)
                              gTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCttt                                                          >  1:1344541/1‑93 (MQ=255)
                               ttGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTc                                                         >  1:3879656/1‑93 (MQ=255)
                               ttGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTc                                                         >  1:6152269/1‑93 (MQ=255)
                               ttGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTc                                                         <  2:3961085/93‑1 (MQ=255)
                                 tgTGCGCTCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCgg                                                       <  2:1382837/92‑1 (MQ=255)
                                 ggTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCgg                                                       >  1:1992192/1‑93 (MQ=255)
                                 ggTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCgg                                                       <  2:3758796/93‑1 (MQ=255)
                                  gTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCggg                                                      <  2:3249003/93‑1 (MQ=255)
                                   tGCGCGCGCTGCTTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGa                                                     <  1:1379688/93‑1 (MQ=255)
                                   tGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGa                                                     <  1:5238772/93‑1 (MQ=255)
                                    gcgcgcgcTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGaa                                                    >  2:3332719/1‑93 (MQ=255)
                                    gcgcgcgcTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGaa                                                    <  2:738756/93‑1 (MQ=255)
                                     cgcgcgcTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAg                                                   <  2:2990188/93‑1 (MQ=255)
                                      gcgcgcTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCACCCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAgg                                                  <  1:602756/93‑1 (MQ=255)
                                      gcgcgcTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAgg                                                  >  1:590764/1‑93 (MQ=255)
                                       cgcgcTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGa                                                 <  1:4655449/93‑1 (MQ=255)
                                       cgcgcTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGa                                                 <  1:2896516/93‑1 (MQ=255)
                                       cgcgcTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGa                                                 <  1:427493/93‑1 (MQ=255)
                                        gcgcTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGAt                                                >  1:3680558/1‑93 (MQ=255)
                                        gcgcTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGAt                                                >  2:5141123/1‑93 (MQ=255)
                                          gcTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGt                                              >  2:4814838/1‑93 (MQ=255)
                                          gcTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGt                                              <  1:5657209/93‑1 (MQ=255)
                                           cTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGtt                                             >  2:4906992/1‑93 (MQ=255)
                                            tGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTc                                            <  2:2510503/93‑1 (MQ=255)
                                            tGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTc                                            <  2:2041072/93‑1 (MQ=255)
                                            tGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTc                                            <  2:596065/93‑1 (MQ=255)
                                             gCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCa                                           <  1:1610624/93‑1 (MQ=255)
                                             gCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCa                                           >  2:2170169/1‑93 (MQ=255)
                                               gTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGc                                         >  1:240870/1‑93 (MQ=255)
                                                tCATCGTTCTCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCa                                        >  2:6118370/1‑93 (MQ=255)
                                                tCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCa                                        >  1:5097340/1‑93 (MQ=255)
                                                tCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGGTCGGAGCGCGTACCATCGGCTTTCGTGAAGGATGTTCAGcc                                        >  1:3138065/1‑92 (MQ=255)
                                                 cATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCAc                                       >  2:5711231/1‑93 (MQ=255)
                                                 cATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCAc                                       >  2:5943055/1‑93 (MQ=255)
                                                  aTCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACg                                      <  1:2781851/93‑1 (MQ=255)
                                                  aTCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACg                                      >  1:1212977/1‑93 (MQ=255)
                                                    cGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGtt                                    >  2:4436978/1‑93 (MQ=255)
                                                      ttATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTa                                  >  2:1811144/1‑93 (MQ=255)
                                                      ttATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTa                                  >  1:1167838/1‑93 (MQ=255)
                                                        aTCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATc                                >  1:1736290/1‑93 (MQ=255)
                                                        aTCCTGATTTACCGATTACCCTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATc                                <  2:5480553/93‑1 (MQ=255)
                                                         tCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCt                               <  2:6246378/93‑1 (MQ=255)
                                                          ccTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTg                              >  1:5332385/1‑93 (MQ=255)
                                                          ccTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTg                              <  2:3998577/93‑1 (MQ=255)
                                                          ccTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTg                              >  1:5711785/1‑93 (MQ=255)
                                                            tGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGcc                            >  2:5348377/1‑93 (MQ=255)
                                                             taTTTTCCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCg                           <  1:3449854/92‑1 (MQ=255)
                                                              aTTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGt                          <  1:4714367/93‑1 (MQ=255)
                                                               tttACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTa                         <  2:1269527/93‑1 (MQ=255)
                                                                ttACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTAt                        <  2:4717061/93‑1 (MQ=255)
                                                                 tACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATg                       <  1:867976/93‑1 (MQ=255)
                                                                 tACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATg                       <  1:2047755/93‑1 (MQ=255)
                                                                  aCCGATTACCGTAACGACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGa                      >  1:2507797/1‑93 (MQ=255)
                                                                  aCCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGa                      >  2:4299026/1‑93 (MQ=255)
                                                                  aCCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGa                      >  1:1528409/1‑93 (MQ=255)
                                                                   acGATTACCGTAACAACCATGAGACCACCCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGTATGTTCAGCACGTTTATCTGCCGTATGAt                     <  1:5786445/92‑1 (MQ=255)
                                                                    cGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATc                    <  1:3625061/93‑1 (MQ=255)
                                                                    cGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATc                    >  2:2088888/1‑93 (MQ=255)
                                                                     gATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCt                   >  2:2900439/1‑93 (MQ=255)
                                                                      aTTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTg                  >  1:5022341/1‑93 (MQ=255)
                                                                       ttACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGc                 >  2:1564608/1‑93 (MQ=255)
                                                                         tccGTAACAACCATGACACCAACCGGTTCGGAACGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGccc               <  2:5193657/92‑1 (MQ=255)
                                                                           cGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGACCGa             >  1:4574314/1‑93 (MQ=255)
                                                                              aacaacCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGc          <  1:4922831/93‑1 (MQ=255)
                                                                                 aacCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACt       <  2:3851931/93‑1 (MQ=255)
                                                                                  acCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTc      >  2:4617232/1‑93 (MQ=255)
                                                                                    cATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCaa    >  1:2288146/1‑93 (MQ=255)
                                                                                     aTGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAAc   <  2:726780/93‑1 (MQ=255)
                                                                                     aTGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAAc   >  1:10315/1‑93 (MQ=255)
                                                                                      tGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAAcc  >  1:3183858/1‑93 (MQ=255)
                                                                                      tGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAAcc  >  1:2683704/1‑93 (MQ=255)
                                                                                          |                                                                                        
CTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACC  >  NC_000913/3806334‑3806512

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