Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 624,133 T→C V9A (GTA→GCA)  entC → isochorismate hydroxymutase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913624,1330TC99.6% 968.5 / ‑5.7 279V9A (GTA→GCA) entCisochorismate hydroxymutase 2
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (140/138);  total (141/138)
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.94e-01

Only 100 of 283 total aligned reads displayed.

GGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGTACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGTTTCGCCCGCT  >  NC_000913/624043‑624222
                                                                                          |                                                                                         
  ttAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCAcag                                                                                       <  1:6350204/93‑1 (MQ=255)
  ttAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCAcag                                                                                       <  2:6012654/93‑1 (MQ=255)
  ttAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCAcag                                                                                       <  1:4138299/93‑1 (MQ=255)
    aGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCAcagta                                                                                     <  1:5020445/93‑3 (MQ=255)
     gCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCAcagcag                                                                                    <  1:114764/93‑1 (MQ=255)
     gCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCAcagcag                                                                                    >  2:2894273/1‑93 (MQ=255)
      cGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGa                                                                                   <  2:4437682/93‑1 (MQ=255)
        aCCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGAcc                                                                                 >  2:6355314/1‑93 (MQ=255)
         ccGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCa                                                                                >  1:4535436/1‑93 (MQ=255)
          cGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCAt                                                                               >  2:4863500/1‑93 (MQ=255)
           gAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATg                                                                              >  2:5353482/1‑93 (MQ=255)
            aaaaTATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATgg                                                                             >  1:3241478/1‑93 (MQ=255)
             aaaTATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGc                                                                            <  1:1862059/93‑1 (MQ=255)
             aaaTATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGc                                                                            >  2:371567/1‑93 (MQ=255)
              aaTATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCa                                                                           <  2:417114/93‑1 (MQ=255)
               atatAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCaa                                                                          <  1:5921467/93‑1 (MQ=255)
                tataAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAac                                                                         >  2:3830678/1‑93 (MQ=255)
                 ataAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAaca                                                                        >  1:871998/1‑93 (MQ=255)
                  taAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAacac                                                                       <  1:6264795/93‑1 (MQ=255)
                  taAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAacac                                                                       >  1:4447138/1‑93 (MQ=255)
                  taAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAacac                                                                       >  2:109443/1‑93 (MQ=255)
                   aaaTGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACt                                                                      <  1:2472212/93‑1 (MQ=255)
                    aaTGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACtt                                                                     >  1:3284442/1‑93 (MQ=255)
                     aTGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTg                                                                    >  2:2498183/1‑93 (MQ=255)
                     aTGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTg                                                                    <  1:1188154/93‑1 (MQ=255)
                       gATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTgcc                                                                  >  1:3575135/1‑92 (MQ=255)
                         tAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGcc                                                                <  1:1309185/93‑1 (MQ=255)
                          aaTCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGccc                                                               <  2:1938700/93‑1 (MQ=255)
                          aaTCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGccc                                                               <  1:4167474/93‑1 (MQ=255)
                          aaTCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGccc                                                               <  2:4565198/93‑1 (MQ=255)
                           aTCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCa                                                              <  1:4616211/93‑1 (MQ=255)
                           aTCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCa                                                              >  1:3194168/1‑93 (MQ=255)
                             cATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAAt                                                            >  2:4222688/1‑93 (MQ=255)
                             cATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAAt                                                            >  2:3416245/1‑93 (MQ=255)
                              attattAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATc                                                           >  2:4338641/1‑93 (MQ=255)
                                 attaAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCt                                                        <  2:550372/93‑1 (MQ=255)
                                 attaAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCt                                                        >  1:283678/1‑93 (MQ=255)
                                  ttaAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCCATCGCtt                                                       >  1:566247/1‑93 (MQ=255)
                                  ttaAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCtt                                                       >  2:2058022/1‑93 (MQ=255)
                                  ttaAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCtt                                                       >  2:3269533/1‑93 (MQ=255)
                                  ttaAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCtt                                                       <  2:5389700/93‑1 (MQ=255)
                                    aaaGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCtttt                                                     >  2:2591545/1‑93 (MQ=255)
                                     aaGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCttttt                                                    >  1:2160909/1‑93 (MQ=255)
                                      aGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTc                                                   <  1:5875464/93‑1 (MQ=255)
                                      aGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTc                                                   >  1:966422/1‑93 (MQ=255)
                                      aGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTc                                                   <  2:3791684/93‑1 (MQ=255)
                                      aGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTc                                                   <  1:101875/93‑1 (MQ=255)
                                      aGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTc                                                   <  2:2296195/93‑1 (MQ=255)
                                      aGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTc                                                   >  2:4042966/1‑93 (MQ=255)
                                      aGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTc                                                   <  2:1911551/93‑1 (MQ=255)
                                       gCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCt                                                  <  1:5271091/93‑1 (MQ=255)
                                         cTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCttt                                                <  2:933986/93‑1 (MQ=255)
                                         cTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCttt                                                >  2:1374757/1‑93 (MQ=255)
                                          tttATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTa                                               <  2:2021912/93‑1 (MQ=255)
                                            tATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATg                                             <  2:3883194/93‑1 (MQ=255)
                                             aTCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGt                                            >  2:1467885/1‑93 (MQ=255)
                                              tCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTc                                           <  1:5522616/93‑1 (MQ=255)
                                              tCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTc                                           >  1:5144752/1‑93 (MQ=255)
                                               cATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCg                                          >  2:4762964/1‑93 (MQ=255)
                                                aTTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGc                                         <  2:564668/93‑1 (MQ=255)
                                                 ttttGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGcc                                        >  2:5264686/1‑93 (MQ=255)
                                                   ttGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGt                                      <  1:2396463/93‑1 (MQ=255)
                                                    tgtgGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTa                                     >  1:4988296/1‑93 (MQ=255)
                                                     gtgGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTAc                                    <  2:6401567/93‑1 (MQ=255)
                                                     gtgGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTAc                                    >  1:1624853/1‑93 (MQ=255)
                                                     gtgGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCCCTTTTTCTTTTTGTCCCCGTAc                                    >  1:3934438/1‑93 (MQ=255)
                                                       ggaggaTGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCg                                  <  1:4894182/93‑1 (MQ=255)
                                                        gaggaTGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGc                                 <  1:4868440/93‑1 (MQ=255)
                                                        gaggaTGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGc                                 <  2:1834080/93‑1 (MQ=255)
                                                         aggaTGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCa                                <  2:4404786/93‑1 (MQ=255)
                                                         aggaTGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCa                                >  2:139223/1‑93 (MQ=255)
                                                           gatgatATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGt                              <  2:1204787/93‑1 (MQ=255)
                                                             tgatATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGttt                            <  1:1940342/93‑1 (MQ=255)
                                                              gatATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGtttt                           <  2:1585628/93‑1 (MQ=255)
                                                              gatATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGtttt                           <  2:307028/93‑1 (MQ=255)
                                                              gatATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGtttt                           <  1:5439243/93‑1 (MQ=255)
                                                                tatGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTac                         >  1:378110/1‑93 (MQ=255)
                                                                  tGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTacgc                       >  1:2070416/1‑92 (MQ=255)
                                                                  tGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTacga                       <  2:697919/93‑1 (MQ=255)
                                                                   ggATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTacgac                      >  1:1735795/1‑93 (MQ=255)
                                                                   ggATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTacgac                      >  2:5211661/1‑93 (MQ=255)
                                                                    gATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTacgacg                     >  2:5099514/1‑93 (MQ=255)
                                                                    gATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTacgacg                     >  1:3358797/1‑93 (MQ=255)
                                                                      tACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTc                   >  2:3947748/1‑93 (MQ=255)
                                                                       aCGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCa                  >  2:424884/1‑93 (MQ=255)
                                                                       aCGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCa                  <  2:3613034/93‑1 (MQ=255)
                                                                         gTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAgg                >  2:3304767/1‑93 (MQ=255)
                                                                         gTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAgg                <  2:5982028/93‑1 (MQ=255)
                                                                          tCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGGTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGa               <  2:3434083/93‑1 (MQ=255)
                                                                             cTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGt            <  2:2408839/93‑1 (MQ=255)
                                                                               ggCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGttt          <  1:541579/93‑1 (MQ=255)
                                                                               ggCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGttt          <  2:5872843/93‑1 (MQ=255)
                                                                                gCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGTTTc         <  2:5843444/93‑1 (MQ=255)
                                                                                    aGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGTTTCGccc     >  2:2149032/1‑93 (MQ=255)
                                                                                    aGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGTTTCGccc     >  2:5030279/1‑93 (MQ=255)
                                                                                     ggAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGTTTCGCCCg    <  1:1558587/93‑1 (MQ=255)
                                                                                      gAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGTTTCGCCCGc   >  2:3104599/1‑93 (MQ=255)
                                                                                      gAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGTTTCGCCCGc   <  1:2905390/93‑1 (MQ=255)
                                                                                       aaGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGTTTCGCCCGCt  >  1:2713287/1‑93 (MQ=255)
                                                                                       aaGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGTTTCGCCCGCt  <  1:5307865/93‑1 (MQ=255)
                                                                                          |                                                                                         
GGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGTACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGTTTCGCCCGCT  >  NC_000913/624043‑624222

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