Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,883,371 T→C Y223H (TAC→CAC)  yidC → 60 KD inner‑membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,883,3710TC99.4% 1560.7 / ‑3.6 491Y223H (TAC→CAC) yidC60 KD inner‑membrane protein
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (183/305);  total (184/307)
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.98e-01

Only 100 of 512 total aligned reads displayed.

TTTGGTCAGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGTACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCGAAAGGTGGTTGG  >  NC_000913/3883281‑3883460
                                                                                          |                                                                                         
tttGGTCAGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCAc                                                                                         >  1:9522948/1‑93 (MQ=255)
  tGGTCAGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTc                                                                                       <  2:144360/93‑1 (MQ=255)
   ggTCAGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTcc                                                                                      >  2:2687038/1‑93 (MQ=255)
    gTCAGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCa                                                                                     <  2:9688044/93‑1 (MQ=255)
    gTCAGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCa                                                                                     <  2:3190647/93‑1 (MQ=255)
    gTCAGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCa                                                                                     <  1:8651613/93‑1 (MQ=255)
     tCAGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCAc                                                                                    <  1:3531229/93‑1 (MQ=255)
       aGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGc                                                                                  <  2:5242836/93‑1 (MQ=255)
       aGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGc                                                                                  >  1:4611365/1‑93 (MQ=255)
        gTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGACGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGcc                                                                                 >  2:1760306/1‑93 (MQ=255)
        gTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGcc                                                                                 >  1:1673732/1‑93 (MQ=255)
         ttGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCt                                                                                >  2:5439745/1‑93 (MQ=255)
           gAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGa                                                                              <  1:2645171/93‑1 (MQ=255)
           gAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGa                                                                              >  2:9372812/1‑93 (MQ=255)
           gAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGa                                                                              <  2:9021897/93‑1 (MQ=255)
             aGCAATCCATCCCTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACg                                                                            <  2:2217213/93‑1 (MQ=255)
             aGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACATTCCGTGGCGCGACGCACTCCACGCCTGACg                                                                            >  2:1491395/1‑93 (MQ=255)
              gCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACga                                                                           <  1:1046516/93‑1 (MQ=255)
               cAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACgag                                                                          >  2:7123826/1‑93 (MQ=255)
                aaTCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACgaga                                                                         >  1:6832971/1‑93 (MQ=255)
                aaTCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACgaga                                                                         <  1:3723315/93‑1 (MQ=255)
                 aTCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGaa                                                                        <  1:199089/93‑1 (MQ=255)
                  tCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAg                                                                       >  1:4337251/1‑93 (MQ=255)
                  tCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAg                                                                       <  1:6388849/93‑1 (MQ=255)
                   ccATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGt                                                                      <  2:595923/93‑1 (MQ=255)
                   ccATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGt                                                                      <  2:3053159/93‑1 (MQ=255)
                   ccATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGt                                                                      <  2:8368849/93‑1 (MQ=255)
                    cATCCCTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTa                                                                     <  1:3041474/93‑1 (MQ=255)
                      tCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATg                                                                   <  2:2501212/93‑1 (MQ=255)
                       cACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATga                                                                  <  1:503060/93‑1 (MQ=255)
                       cACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATga                                                                  >  1:1191370/1‑93 (MQ=255)
                       cACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATga                                                                  >  1:2532018/1‑93 (MQ=255)
                        aCTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATgag                                                                 >  1:5976429/1‑93 (MQ=255)
                          tctGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGaa                                                               >  2:789415/1‑93 (MQ=255)
                           ctGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGaaa                                                              >  1:3945141/1‑93 (MQ=255)
                            tGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAAt                                                             >  2:7391272/1‑93 (MQ=255)
                             gCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATa                                                            >  2:3545936/1‑93 (MQ=255)
                             gCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATa                                                            <  1:2112495/93‑1 (MQ=255)
                              ccACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAAtgc                                                           >  1:4429454/1‑91 (MQ=255)
                                 ccGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAg                                                        >  2:9621591/1‑93 (MQ=255)
                                 ccGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAg                                                        >  1:4877889/1‑93 (MQ=255)
                                 ccGCATCTCGAGACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGGGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAg                                                        <  1:5096216/93‑1 (MQ=255)
                                    cATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTc                                                     <  1:1511136/93‑1 (MQ=255)
                                    cATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTc                                                     <  1:2795052/93‑1 (MQ=255)
                                     aTCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCg                                                    >  1:3174980/1‑93 (MQ=255)
                                      tctcGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGa                                                   >  2:4102576/1‑93 (MQ=255)
                                      actcGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGa                                                   <  1:7123582/92‑1 (MQ=255)
                                       ctcGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGAt                                                  <  1:5607855/93‑1 (MQ=255)
                                       ctcGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGAt                                                  <  2:8658641/93‑1 (MQ=255)
                                          gATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATAcc                                               <  1:2303882/93‑1 (MQ=255)
                                            tACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCAt                                             <  1:7695226/93‑1 (MQ=255)
                                             aCCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCAtt                                            >  1:9368316/1‑93 (MQ=255)
                                             aCCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCAtt                                            <  1:8065951/93‑1 (MQ=255)
                                               cGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGc                                          <  1:613429/93‑1 (MQ=255)
                                                ggAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGcc                                         <  2:6702695/93‑1 (MQ=255)
                                                 gAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCg                                        >  1:7256184/1‑93 (MQ=255)
                                                   agcagcAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGAt                                      >  2:7611753/1‑93 (MQ=255)
                                                   agcagcAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGAt                                      <  2:6544493/93‑1 (MQ=255)
                                                   agcagcAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGAt                                      >  1:7822527/1‑93 (MQ=255)
                                                   agcagcAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACACCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGAt                                      <  1:182883/93‑1 (MQ=255)
                                                    gcagcaACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATa                                     <  2:5607967/93‑1 (MQ=255)
                                                    gcagcaACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATa                                     <  2:3931308/93‑1 (MQ=255)
                                                     cagcaACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATaa                                    <  2:6863436/93‑1 (MQ=255)
                                                      agcaACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAAc                                   <  2:5075896/93‑1 (MQ=255)
                                                        caACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGa                                 <  2:1885525/93‑1 (MQ=255)
                                                         aaCTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGaa                                <  2:8127454/93‑1 (MQ=255)
                                                           cTTCGCACTGCACACCTTCCGTGGCGCGGCGCGCTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGaaaa                              <  2:3309296/93‑1 (MQ=255)
                                                           cTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGaaaa                              >  2:1219718/1‑93 (MQ=255)
                                                            ttCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAAc                             >  1:3347928/1‑93 (MQ=255)
                                                            ttCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAAc                             >  1:5742356/1‑93 (MQ=255)
                                                               gCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTg                          <  2:4225964/93‑1 (MQ=255)
                                                                cACTGCACACCTTCCGTGGCGCGGCGTACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGa                         <  1:9208179/93‑1 (MQ=255)
                                                                  cTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAAc                       <  2:952634/93‑1 (MQ=255)
                                                                   tGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACa                      >  2:4810046/1‑93 (MQ=255)
                                                                    gCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACAt                     >  1:5288833/1‑93 (MQ=255)
                                                                      acacCTGCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACAtct                   <  1:6993916/93‑1 (MQ=255)
                                                                       cacCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACAtctc                  <  1:2938013/93‑1 (MQ=255)
                                                                       cacCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACAtctc                  >  1:5416286/1‑93 (MQ=255)
                                                                        acCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATctct                 >  1:9716963/1‑93 (MQ=255)
                                                                        acCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATctct                 <  2:238171/93‑1 (MQ=255)
                                                                        acCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATctct                 <  1:6322829/93‑1 (MQ=255)
                                                                        acCTTCCGTGGCGCGGCGCACTCCACGCATGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATctct                 >  1:6612132/1‑93 (MQ=255)
                                                                         ccTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCtt                >  1:8304766/1‑93 (MQ=255)
                                                                         ccTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCtt                <  1:4178018/93‑1 (MQ=255)
                                                                         ccTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCtt                <  1:65205/93‑1 (MQ=255)
                                                                          cTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTc               <  1:9575922/93‑1 (MQ=255)
                                                                           ttCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCg              <  2:377912/93‑1 (MQ=255)
                                                                           ttCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCg              >  2:281265/1‑93 (MQ=255)
                                                                           ttCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCg              >  2:7703989/1‑93 (MQ=255)
                                                                              cGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCGaaa           <  1:6173712/93‑1 (MQ=255)
                                                                              cGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCGaaa           <  2:799065/93‑1 (MQ=255)
                                                                              cGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCGaaa           <  1:1396889/93‑1 (MQ=255)
                                                                               gTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATTTCTTCGAAAg          >  2:4401927/1‑93 (MQ=255)
                                                                               gTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCGAAAg          <  2:5351474/93‑1 (MQ=255)
                                                                                 ggcgcggcgcACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCGAAAggg        >  2:4327330/1‑92 (MQ=255)
                                                                                   cgcgGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCGAAAggtgg      <  1:9917238/93‑1 (MQ=255)
                                                                                    gcggcgCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCGAAAggtggt     <  2:3671248/93‑1 (MQ=255)
                                                                                     cggcgCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGACCATCTCTTCGAAAGGTGGtt    <  1:8867644/93‑1 (MQ=255)
                                                                                     cggcgCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCGAAAGGTGGtt    <  2:8493743/93‑1 (MQ=255)
                                                                                      ggcgCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCGAAAGGTGGTTg   >  1:6918176/1‑93 (MQ=255)
                                                                                          |                                                                                         
TTTGGTCAGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGTACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCGAAAGGTGGTTGG  >  NC_000913/3883281‑3883460

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