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,1330TC86.8% 711.0 / 84.2 287V9A (GTA→GCA) entCisochorismate hydroxymutase 2
Reads supporting (aligned to +/- strand):  ref base T (22/16);  new base C (132/117);  total (154/133)
Fisher's exact test for biased strand distribution p-value = 6.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 291 total aligned reads displayed.

GGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGTACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGTTTCGCCCGCTTCG  >  NC_000913/624043‑624225
                                                                                          |                                                                                            
ggTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGTAc                                                                                            >  1:1246763/1‑93 (MQ=255)
ggTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCAc                                                                                            >  1:2266682/1‑93 (MQ=255)
ggTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCAc                                                                                            >  2:5785206/1‑93 (MQ=255)
 gTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGcaca                                                                                           <  1:2573540/93‑1 (MQ=255)
 gTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGcaca                                                                                           <  2:5948587/93‑1 (MQ=255)
   tAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCAcagc                                                                                         <  2:4405351/93‑1 (MQ=255)
      cGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGa                                                                                      <  2:3614408/93‑1 (MQ=255)
       gACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGAc                                                                                     >  1:1373209/1‑93 (MQ=255)
       gACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGAc                                                                                     >  2:698498/1‑93 (MQ=255)
         ccGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGTACAGCAGACCa                                                                                   >  1:447070/1‑93 (MQ=255)
         ccGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCa                                                                                   <  1:59035/93‑1 (MQ=255)
         ccGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCa                                                                                   >  2:2015322/1‑93 (MQ=255)
         ccGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCa                                                                                   <  1:5218697/93‑1 (MQ=255)
          cGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCAt                                                                                  >  2:2375184/1‑93 (MQ=255)
           gAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATg                                                                                 >  1:1291952/1‑93 (MQ=255)
           gAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATg                                                                                 >  2:5767637/1‑93 (MQ=255)
           gAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATg                                                                                 >  1:6441796/1‑93 (MQ=255)
            aaaaTATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATgg                                                                                >  1:1449399/1‑93 (MQ=255)
            aaaaTATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATgg                                                                                <  2:4468157/93‑1 (MQ=255)
            aaaaTATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATgg                                                                                <  1:4834003/93‑1 (MQ=255)
             aaaTATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGc                                                                               >  2:4221187/1‑93 (MQ=255)
               atatAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCaa                                                                             >  1:513988/1‑93 (MQ=255)
               atatAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCaa                                                                             <  2:6611078/93‑1 (MQ=255)
               atatAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCaa                                                                             <  1:1852588/93‑1 (MQ=255)
                tataAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAac                                                                            >  1:1547093/1‑93 (MQ=255)
                tataAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAac                                                                            <  2:2515167/93‑1 (MQ=255)
                 ataAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGTACAGCAGACCATGGCAaca                                                                           <  2:382102/93‑1 (MQ=255)
                 ataAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAaca                                                                           >  1:2367728/1‑93 (MQ=255)
                  taAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGTACAGCAGACCATGGCAacac                                                                          <  1:4117561/93‑1 (MQ=255)
                    aaTGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGTACAGCAGACCATGGCAACACtt                                                                        >  1:6322827/1‑93 (MQ=255)
                    aaTGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACtt                                                                        <  2:2826983/93‑1 (MQ=255)
                    aaTGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACtt                                                                        <  2:1383887/93‑1 (MQ=255)
                     aTGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGTACAGCAGACCATGGCAACACTTg                                                                       >  2:755904/1‑93 (MQ=255)
                     aTGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTg                                                                       >  1:2917276/1‑93 (MQ=255)
                       gATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTgcg                                                                     >  2:1965490/1‑93 (MQ=255)
                       gATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTgcg                                                                     <  1:6208252/93‑1 (MQ=255)
                         tAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGcc                                                                   <  1:6878621/93‑1 (MQ=255)
                         tAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGcc                                                                   <  1:2162409/93‑1 (MQ=255)
                          aaTCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGccc                                                                  >  2:4844394/1‑93 (MQ=255)
                          aaTCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGccc                                                                  >  1:5826543/1‑93 (MQ=255)
                            tCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCaa                                                                >  2:2181596/1‑93 (MQ=255)
                              attattAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGTACAGCAGACCATGGCAACACTTGCGCCCAATc                                                              <  1:3362863/93‑1 (MQ=255)
                              attattAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATc                                                              >  2:6546510/1‑93 (MQ=255)
                              attattAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATc                                                              >  1:6557864/1‑93 (MQ=255)
                               ttattaAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCg                                                             >  1:1891266/1‑93 (MQ=255)
                               ttattaAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCg                                                             <  1:3024847/93‑1 (MQ=255)
                                tattaAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGTACAGCAGACCATGGCAACACTTGCGCCCAATCGc                                                            <  1:1508642/93‑1 (MQ=255)
                                 attaAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCt                                                           >  1:2858074/1‑93 (MQ=255)
                                   taAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCttt                                                         <  1:5814752/93‑1 (MQ=255)
                                   taAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCttt                                                         <  2:3403974/93‑1 (MQ=255)
                                    aaaGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCtttt                                                        <  2:1157260/93‑1 (MQ=255)
                                     aaGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCttttt                                                       >  1:5444537/1‑93 (MQ=255)
                                     aaGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCttttt                                                       <  2:7414364/93‑1 (MQ=255)
                                     aaGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCttttt                                                       >  2:7312554/1‑93 (MQ=255)
                                       gCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCt                                                     <  1:1038746/93‑1 (MQ=255)
                                       gCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCt                                                     <  2:3985213/93‑1 (MQ=255)
                                          tttATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTa                                                  >  2:5049205/1‑93 (MQ=255)
                                          tttATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTa                                                  <  2:1627694/93‑1 (MQ=255)
                                          tttATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTa                                                  <  1:2420348/93‑1 (MQ=255)
                                            tATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAGCACTTGCGCCCAATCGCTTTTTCTTTATg                                                >  2:1067203/1‑93 (MQ=255)
                                            tATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATg                                                >  1:7489728/1‑93 (MQ=255)
                                            tATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATg                                                <  1:4274318/93‑1 (MQ=255)
                                            tATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATg                                                >  1:3633442/1‑93 (MQ=255)
                                              tCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTc                                              <  2:6745124/93‑1 (MQ=255)
                                              tCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTc                                              <  2:5696950/93‑1 (MQ=255)
                                               cATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCg                                             >  1:5064026/1‑93 (MQ=255)
                                                aTTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGc                                            <  1:7242131/93‑1 (MQ=255)
                                                aTTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGc                                            <  1:3474025/93‑1 (MQ=255)
                                                aTTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGc                                            >  1:4353679/1‑93 (MQ=255)
                                                aTTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGc                                            >  1:4505596/1‑93 (MQ=255)
                                                 ttttGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGcc                                           >  2:609373/1‑93 (MQ=255)
                                                 ttttGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGcc                                           >  1:3958649/1‑93 (MQ=255)
                                                       ggaggaTGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCg                                     >  2:3257039/1‑93 (MQ=255)
                                                        gaggaTGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGc                                    >  2:6282007/1‑93 (MQ=255)
                                                         aggaTGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCa                                   <  2:4693618/93‑1 (MQ=255)
                                                           gatgatATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGt                                 <  2:5353789/93‑1 (MQ=255)
                                                             tgatATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGttt                               >  1:4220599/1‑93 (MQ=255)
                                                              gatATGGATACGTCACTGGCTGAGGAAGTACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGtttt                              >  2:5249662/1‑93 (MQ=255)
                                                                tatGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTac                            >  1:1667693/1‑93 (MQ=255)
                                                                  tGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTacga                          >  1:6089665/1‑93 (MQ=255)
                                                                  tGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTacga                          <  1:4694441/93‑1 (MQ=255)
                                                                     aTACGTCACTGGCTGAGGAAGTACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGt                       <  1:6695693/93‑1 (MQ=255)
                                                                     aTACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGt                       <  1:4387636/93‑1 (MQ=255)
                                                                     aTACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGt                       >  1:7463156/1‑93 (MQ=255)
                                                                      tACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTc                      >  1:881271/1‑93 (MQ=255)
                                                                      tACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTc                      >  1:1628364/1‑93 (MQ=255)
                                                                       aCGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCa                     <  2:4944040/93‑1 (MQ=255)
                                                                        cGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAg                    <  1:1424580/93‑1 (MQ=255)
                                                                        cGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAg                    <  2:684349/93‑1 (MQ=255)
                                                                          tCACTGGCTGAGGAAGTACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGa                  <  1:2997839/93‑1 (MQ=255)
                                                                          tCACTGGCTGAGGAAGTACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGa                  <  1:3208156/93‑1 (MQ=255)
                                                                          tCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGa                  <  2:3039749/93‑1 (MQ=255)
                                                                          tCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGa                  <  2:152630/93‑1 (MQ=255)
                                                                           cACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGAt                 <  1:4575101/93‑1 (MQ=255)
                                                                            aCTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATg                <  1:5374390/93‑1 (MQ=255)
                                                                             cTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGt               <  2:260906/93‑1 (MQ=255)
                                                                              tGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTATGACGTCAGGATGtt              <  2:7259745/93‑1 (MQ=255)
                                                                              tGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGtt              <  1:2743059/93‑1 (MQ=255)
                                                                               ggCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGttt             >  1:100494/1‑93 (MQ=255)
                                                                                   gAGGAAGTACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGTTTCGcc         <  2:5501105/93‑1 (MQ=255)
                                                                                          |                                                                                            
GGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGTACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGTTTCGCCCGCTTCG  >  NC_000913/624043‑624225

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