Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,213,360 A→G M102V (ATG→GTG)  aceB → malate synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,213,3600AG99.7% 1010.4 / NA 324M102V (ATG→GTG) aceBmalate synthase A
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (148/174);  total (150/174)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 339 total aligned reads displayed.

TCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGAC  >  NC_000913/4213270‑4213449
                                                                                          |                                                                                         
 ccATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTgg                                                                                        <  1:5765297/93‑1 (MQ=255)
  cATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAggtggt                                                                                       >  2:2667656/1‑93 (MQ=255)
  cATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAggtggt                                                                                       <  1:4492985/93‑1 (MQ=255)
  cATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAggtggt                                                                                       >  2:5150602/1‑93 (MQ=255)
   aTTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGgtggtg                                                                                      >  2:1651681/1‑93 (MQ=255)
    ttCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGa                                                                                     <  1:592422/93‑1 (MQ=255)
     tCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGAt                                                                                    <  1:5591502/93‑1 (MQ=255)
      cgcgATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATc                                                                                   <  1:3941072/93‑1 (MQ=255)
       gcgATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCa                                                                                  >  1:6995369/1‑93 (MQ=255)
          aTGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAAcg                                                                               <  2:7526016/93‑1 (MQ=255)
          aTGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAAcg                                                                               >  1:7597152/1‑93 (MQ=255)
           tGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAAcgc                                                                              >  2:4949138/1‑93 (MQ=255)
           tGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAAcgc                                                                              >  2:6637761/1‑93 (MQ=255)
            gCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAAcgcg                                                                             >  1:2731703/1‑93 (MQ=255)
               gATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCACGGTAGAGCGCAAGGTGGTGATCAACGCGCTc                                                                          <  2:6242443/93‑1 (MQ=255)
                aTTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCa                                                                         <  1:6908522/93‑1 (MQ=255)
                 ttGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCaa                                                                        >  1:5012575/1‑93 (MQ=255)
                 ttGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCaa                                                                        >  2:4549196/1‑93 (MQ=255)
                  tGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAAc                                                                       <  1:146524/93‑1 (MQ=255)
                   ggAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGTTCAACg                                                                      <  1:1840773/93‑1 (MQ=255)
                    gAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGc                                                                     >  2:5761989/1‑93 (MQ=255)
                       aaTTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCaa                                                                  <  1:1774483/93‑1 (MQ=255)
                       aaTTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACACGCTCAACGCCaa                                                                  >  1:7881190/1‑93 (MQ=255)
                        aTTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAAt                                                                 <  2:5577050/93‑1 (MQ=255)
                         ttCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAAtg                                                                >  2:6112472/1‑93 (MQ=255)
                          tCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAAtgt                                                               <  2:438142/93‑1 (MQ=255)
                            gtgggATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGa                                                             <  2:6513216/91‑1 (MQ=255)
                             cgGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGaa                                                            <  1:938180/93‑1 (MQ=255)
                             cgGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGaa                                                            >  2:4668506/1‑93 (MQ=255)
                               ggATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAg                                                          <  2:7361899/93‑1 (MQ=255)
                               ggATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAg                                                          >  1:4792649/1‑93 (MQ=255)
                               ggATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAg                                                          >  2:4095428/1‑93 (MQ=255)
                               ggATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAg                                                          >  2:3960635/1‑93 (MQ=255)
                                 aTTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTc                                                        >  2:1959696/1‑93 (MQ=255)
                                 aTTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTc                                                        >  2:1348124/1‑93 (MQ=255)
                                  ttCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCt                                                       <  2:3544706/93‑1 (MQ=255)
                                  ttCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCt                                                       >  1:7670623/1‑93 (MQ=255)
                                  ttCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCt                                                       >  1:7418012/1‑93 (MQ=255)
                                     cTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTa                                                    <  1:3234618/93‑1 (MQ=255)
                                     cTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTa                                                    <  2:6336743/93‑1 (MQ=255)
                                     cTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTa                                                    <  2:6766600/93‑1 (MQ=255)
                                      tGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTAt                                                   >  1:7514161/1‑93 (MQ=255)
                                      tGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTAt                                                   <  2:2814621/93‑1 (MQ=255)
                                       gCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATg                                                  <  2:5329901/93‑1 (MQ=255)
                                       gCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATg                                                  >  2:1044153/1‑93 (MQ=255)
                                        cGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATgg                                                 <  1:5355766/93‑1 (MQ=255)
                                         ggACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGc                                                <  1:246828/93‑1 (MQ=255)
                                          gACTTAGAAGACCGGCGCGTAGAGATAACTGGCCAGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGcc                                               <  2:320688/93‑1 (MQ=255)
                                             ttAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTGAACGCCAATGTGAAAGTCTTTATGGCCgct                                            >  2:2191893/1‑91 (MQ=255)
                                             ttAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGAt                                            <  2:6965490/93‑1 (MQ=255)
                                             ttAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGAt                                            >  2:6020301/1‑93 (MQ=255)
                                               agaagaCCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGAttt                                          <  1:4184210/93‑1 (MQ=255)
                                                gaagaCCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGGGATCAACGCGCTCACCGCCAATGTGAAAGTCTTTATGGCCGGtttg                                         >  1:2340126/1‑92 (MQ=255)
                                                  agaCCGCCGCGGAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGa                                       <  1:1200737/93‑1 (MQ=255)
                                                   gaCCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGaa                                      >  2:5809169/1‑93 (MQ=255)
                                                   gaCCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGaa                                      >  2:232436/1‑93 (MQ=255)
                                                    aCCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAg                                     >  2:1843336/1‑93 (MQ=255)
                                                        ccgcGTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTc                                 <  2:3230025/93‑1 (MQ=255)
                                                         cgcgTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCa                                <  1:6776461/93‑1 (MQ=255)
                                                         cgcgTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCa                                <  1:2988311/93‑1 (MQ=255)
                                                          gcgTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCAc                               <  1:7712407/93‑1 (MQ=255)
                                                          gcgTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCAc                               >  1:6571495/1‑93 (MQ=255)
                                                          gcgTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCAc                               >  1:3071029/1‑93 (MQ=255)
                                                          gcgTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCAc                               <  2:2159737/93‑1 (MQ=255)
                                                           cgTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACt                              >  1:1038558/1‑93 (MQ=255)
                                                            gTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTg                             <  2:3012299/93‑1 (MQ=255)
                                                            gTAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTg                             <  2:1154927/93‑1 (MQ=255)
                                                             tAGAGATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTgg                            >  1:4036907/1‑93 (MQ=255)
                                                             tAGAGATAACTGGACCGGTAGAGCGCAAGCTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTgg                            >  2:6844758/1‑93 (MQ=255)
                                                              agagATAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGc                           >  1:3176464/1‑93 (MQ=255)
                                                               gagaTAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCa                          >  1:3499279/1‑93 (MQ=255)
                                                                agaTAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCAc                         <  1:7123422/93‑1 (MQ=255)
                                                                 gaTAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCAcc                        >  2:3720545/1‑93 (MQ=255)
                                                                 gaTAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCAcc                        >  2:7548150/1‑93 (MQ=255)
                                                                  aTAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCa                       <  2:1888106/93‑1 (MQ=255)
                                                                   tAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAg                      >  2:5202009/1‑93 (MQ=255)
                                                                   tAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAg                      <  1:7355553/93‑1 (MQ=255)
                                                                   tAACTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAg                      >  2:7262024/1‑93 (MQ=255)
                                                                    aaCTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGa                     <  1:1295556/93‑1 (MQ=255)
                                                                     aCTGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGAc                    >  2:1607548/1‑93 (MQ=255)
                                                                       tGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGAcag                  >  1:6151056/1‑91 (MQ=255)
                                                                       tGGCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTg                  <  1:3831156/93‑1 (MQ=255)
                                                                        ggCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTgg                 >  1:2079243/1‑93 (MQ=255)
                                                                        ggCCCGGTAGAGCGAAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTgg                 <  2:7880236/93‑1 (MQ=255)
                                                                         gCCCGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGa                <  2:3470292/93‑1 (MQ=255)
                                                                          cccGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGaa               <  2:2894456/93‑1 (MQ=255)
                                                                          cccGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGaa               >  2:1469112/1‑93 (MQ=255)
                                                                            cGGTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACa             >  2:4781831/1‑93 (MQ=255)
                                                                             ggTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACaa            >  1:6051878/1‑93 (MQ=255)
                                                                             ggTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACaa            >  1:2249740/1‑93 (MQ=255)
                                                                              gTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACaaa           <  1:5247825/93‑1 (MQ=255)
                                                                              gTAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACaaa           >  2:1255636/1‑93 (MQ=255)
                                                                               tAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAg          <  2:6461120/93‑1 (MQ=255)
                                                                               tAGAGCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAg          <  1:3859209/93‑1 (MQ=255)
                                                                               tAGAGCGCAAGATGGAGGTCAACGTGCTCAACGCCAATGTGAAAGACTTTATGGCCGAGTTAGAAGATTCACTGCCACCTGAATAGCACaaaa          >  2:7651319/1‑92 (MQ=255)
                                                                                agagCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGt         <  1:3273967/93‑1 (MQ=255)
                                                                                  agCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGa       <  1:4394772/93‑1 (MQ=255)
                                                                                  agCGCAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGa       <  2:3330438/93‑1 (MQ=255)
                                                                                   gcgcAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGAt      <  1:3744766/93‑1 (MQ=255)
                                                                                      cAAGGTGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGa   <  2:843294/93‑1 (MQ=255)
                                                                                          |                                                                                         
TCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGAC  >  NC_000913/4213270‑4213449

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: