Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 459,477 T→C Y456H (TAC→CAC)  lon → ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913459,4770TC99.7% 918.7 / NA 282Y456H (TAC→CAC) lonATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (156/125);  minor base A (1/0);  total (157/125)
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 303 total aligned reads displayed.

GACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATG  >  NC_000913/459387‑459566
                                                                                          |                                                                                         
 aCAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGaccaccacc                                                                                        >  1:2226516/1‑93 (MQ=255)
 aCAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGaccaccacc                                                                                        >  1:792402/1‑93 (MQ=255)
   aaaaTGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTg                                                                                      >  2:395313/1‑93 (MQ=255)
     aaTGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGa                                                                                    >  1:5499534/1‑93 (MQ=255)
      aTGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGaa                                                                                   <  1:4464571/93‑1 (MQ=255)
      aTGTCTTCTGACATGCGTGGCGATCCAGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGaa                                                                                   >  1:5545437/1‑93 (MQ=255)
       tGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAg                                                                                  >  2:1883099/1‑93 (MQ=255)
         tcttctGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTg                                                                                <  1:4641413/93‑1 (MQ=255)
         tcttctGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTg                                                                                >  2:6292652/1‑93 (MQ=255)
          cttctGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTgg                                                                               <  1:657857/93‑1 (MQ=255)
             ctGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGAtt                                                                            <  2:6179372/93‑1 (MQ=255)
               gACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTAc                                                                          >  1:1722406/1‑93 (MQ=255)
               gACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTAc                                                                          <  2:3873016/93‑1 (MQ=255)
               gACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTAc                                                                          >  1:1095671/1‑93 (MQ=255)
               gACATGCGGGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTAc                                                                          <  1:610896/93‑1 (MQ=255)
                 cATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGa                                                                        >  1:1482800/1‑93 (MQ=255)
                  aTGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGAt                                                                       >  1:3464159/1‑93 (MQ=255)
                    gCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGAtct                                                                     <  2:1065535/93‑1 (MQ=255)
                      gTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCa                                                                   <  1:5008948/93‑1 (MQ=255)
                      gTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCa                                                                   >  2:2796678/1‑93 (MQ=255)
                      gTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCa                                                                   >  2:850048/1‑93 (MQ=255)
                      gTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCa                                                                   <  1:3160398/93‑1 (MQ=255)
                       tGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAg                                                                  >  1:6690188/1‑93 (MQ=255)
                       tGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAg                                                                  >  2:5625730/1‑93 (MQ=255)
                         gCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCg                                                                <  1:4180504/93‑1 (MQ=255)
                           gATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGAc                                                              >  1:4776972/1‑93 (MQ=255)
                           gATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGAc                                                              >  2:5569419/1‑93 (MQ=255)
                              ccGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTg                                                           <  1:3275782/93‑1 (MQ=255)
                              ccGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTg                                                           >  1:4006623/1‑93 (MQ=255)
                               cGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGa                                                          <  2:2413359/93‑1 (MQ=255)
                                ggCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGAt                                                         >  2:2705966/1‑93 (MQ=255)
                                ggCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGAt                                                         <  1:1798205/93‑1 (MQ=255)
                                 gCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATg                                                        >  1:413376/1‑93 (MQ=255)
                                 gCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATg                                                        >  2:4144826/1‑93 (MQ=255)
                                    tctGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGttt                                                     >  1:3363042/1‑93 (MQ=255)
                                    tctGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGttt                                                     >  2:3516133/1‑93 (MQ=255)
                                     ctGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTg                                                    <  2:2104722/93‑1 (MQ=255)
                                      tGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGt                                                   <  2:5479395/93‑1 (MQ=255)
                                      tGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGt                                                   >  1:2006809/1‑93 (MQ=255)
                                      tGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGt                                                   <  1:934619/93‑1 (MQ=255)
                                       gCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTc                                                  <  1:3001044/93‑1 (MQ=255)
                                       gCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTc                                                  <  1:3557554/93‑1 (MQ=255)
                                       gCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTc                                                  <  1:5958270/93‑1 (MQ=255)
                                         aCTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTcgc                                                >  1:1269111/1‑93 (MQ=255)
                                         aCTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTcgc                                                <  2:1772390/93‑1 (MQ=255)
                                          cTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTcgcg                                               >  2:4534589/1‑93 (MQ=255)
                                          cTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTcgcg                                               >  1:5145161/1‑93 (MQ=255)
                                          cTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTcgcg                                               >  2:1981646/1‑93 (MQ=255)
                                           tGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGa                                              >  1:1032996/1‑93 (MQ=255)
                                           tGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGa                                              >  2:3658938/1‑93 (MQ=255)
                                            gtttGAAGTGCTGGATCCACAGCATAACGTAGCGATCAGCGACCACCACCTGGAAGAGGATTACGATCTCTGCGATGGGATGTTTGTcgcgcc                                             >  1:636248/3‑91 (MQ=255)
                                            gCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGAc                                             >  2:5787375/1‑93 (MQ=255)
                                            gCTTGAAGTGCCGGAACCAAAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGAc                                             <  2:4430413/93‑1 (MQ=255)
                                             cTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACg                                            <  1:1348201/93‑1 (MQ=255)
                                              ttGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGt                                           <  2:6594806/93‑1 (MQ=255)
                                              ttGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGt                                           >  1:2591746/1‑93 (MQ=255)
                                               tGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTc                                          >  2:214395/1‑93 (MQ=255)
                                                 aaGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGa                                        <  2:6070597/93‑1 (MQ=255)
                                                 aaGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGa                                        >  2:3426305/1‑93 (MQ=255)
                                                   gTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAAc                                      <  1:2086980/93‑1 (MQ=255)
                                                    tGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACt                                     <  2:1331181/93‑1 (MQ=255)
                                                     gCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTc                                    <  2:2622971/93‑1 (MQ=255)
                                                        ggATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCAt                                 <  1:1872735/93‑1 (MQ=255)
                                                          aTCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGa                               >  1:289914/1‑93 (MQ=255)
                                                          aTCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGa                               <  2:4648329/93‑1 (MQ=255)
                                                           tCCAGAGCAGAACGTAGCGTTCAGCGACCCCCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGaa                              <  1:2758177/93‑1 (MQ=255)
                                                             cAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACa                            >  2:3696083/1‑93 (MQ=255)
                                                              agagCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACAt                           >  1:486082/1‑93 (MQ=255)
                                                              agagCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACAt                           <  2:3105446/93‑1 (MQ=255)
                                                               gagCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACAtt                          >  2:4867339/1‑93 (MQ=255)
                                                               gagCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACAtt                          >  2:1525278/1‑93 (MQ=255)
                                                               gagCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACATGATGTTTGTCGCGACGTCGAACTCAATGAACAtt                          >  2:1239940/1‑93 (MQ=255)
                                                                agCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTc                         <  1:4329690/93‑1 (MQ=255)
                                                                 gCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTcc                        >  2:3240451/1‑93 (MQ=255)
                                                                  cAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCg                       <  2:2516027/93‑1 (MQ=255)
                                                                  cAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCg                       >  1:1511723/1‑93 (MQ=255)
                                                                   aGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCgg                      >  1:5244375/1‑93 (MQ=255)
                                                                    gAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGc                     >  1:5736292/1‑93 (MQ=255)
                                                                    gAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGc                     <  2:1698301/93‑1 (MQ=255)
                                                                      aCGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCAc                   >  2:141659/1‑93 (MQ=255)
                                                                        gTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCg                 <  1:2513240/93‑1 (MQ=255)
                                                                         tAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCgc                <  1:2301976/93‑1 (MQ=255)
                                                                         tAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCgc                <  2:218709/93‑1 (MQ=255)
                                                                         tAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCgc                <  1:2367272/93‑1 (MQ=255)
                                                                           gCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCgctg              >  1:5634766/1‑93 (MQ=255)
                                                                           gCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCgctg              >  1:2166774/1‑93 (MQ=255)
                                                                           gCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCgctg              >  2:5318332/1‑93 (MQ=255)
                                                                            cGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCgctgc             <  1:6225110/93‑1 (MQ=255)
                                                                             gTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCgctgct            <  1:730502/93‑1 (MQ=255)
                                                                             gTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCgctgct            >  2:4990199/1‑93 (MQ=255)
                                                                               tCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTACGGCACCGCTGCTgg          >  2:1641232/1‑93 (MQ=255)
                                                                                cAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGa         <  2:764796/93‑1 (MQ=255)
                                                                                cAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGa         <  1:318990/93‑1 (MQ=255)
                                                                                 aGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGAt        <  1:5366290/93‑1 (MQ=255)
                                                                                  gCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATc       >  2:360436/1‑93 (MQ=255)
                                                                                  gCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATc       <  1:3509882/93‑1 (MQ=255)
                                                                                   cGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTTCGGCACCGCTGCTGGATCg      <  1:6472365/93‑1 (MQ=255)
                                                                                   cGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCg      >  2:5091212/1‑93 (MQ=255)
                                                                                     accaccaccTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTa    >  1:3781273/1‑93 (MQ=255)
                                                                                      ccaccaCCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTAt   >  2:5992692/1‑93 (MQ=255)
                                                                                          |                                                                                         
GACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATG  >  NC_000913/459387‑459566

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