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.4% 502.3 / ‑5.4 158Y456H (TAC→CAC) lonATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (91/66);  minor base G (1/0);  total (92/66)
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 177 total aligned reads displayed.

GACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATG  >  NC_000913/459387‑459566
                                                                                          |                                                                                         
gACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACcaccac                                                                                         <  1:1536914/93‑1 (MQ=255)
 aCAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGaccaccacc                                                                                        <  1:3269964/93‑1 (MQ=255)
  cAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCt                                                                                       <  2:1294990/93‑1 (MQ=255)
   aaaaTGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTg                                                                                      <  2:664675/93‑1 (MQ=255)
   aaaaTGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTg                                                                                      >  2:1427228/1‑93 (MQ=255)
     aaTGTCTTCTTACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGa                                                                                    <  2:3373284/93‑1 (MQ=255)
     aaTGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGa                                                                                    >  2:2074350/1‑93 (MQ=255)
       tGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAg                                                                                  >  2:37050/1‑93 (MQ=255)
         tcttctGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTg                                                                                >  1:2448426/1‑93 (MQ=255)
           ttgtGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGa                                                                              <  2:3636996/90‑1 (MQ=255)
           ttctGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGa                                                                              >  1:2777985/1‑93 (MQ=255)
           ttctGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGa                                                                              >  1:645386/1‑93 (MQ=255)
            tctGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGAt                                                                             >  1:703029/1‑93 (MQ=255)
             ctGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGAtt                                                                            <  1:517660/93‑1 (MQ=255)
              tGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTa                                                                           >  2:2526232/1‑93 (MQ=255)
              tGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTa                                                                           <  2:1629032/93‑1 (MQ=255)
              tGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTa                                                                           >  2:1727662/1‑93 (MQ=255)
              tGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGGGGATTa                                                                           >  1:2863022/1‑93 (MQ=255)
                aCATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACg                                                                         >  1:2207166/1‑93 (MQ=255)
                 cATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGa                                                                        >  1:1712445/1‑93 (MQ=255)
                 cATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGa                                                                        <  2:2185459/93‑1 (MQ=255)
                  aTGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGAt                                                                       >  2:4484836/1‑93 (MQ=255)
                   tGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGAtc                                                                      >  1:4570401/1‑93 (MQ=255)
                   tGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGAtc                                                                      >  2:4407527/1‑93 (MQ=255)
                   tGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGAtc                                                                      >  2:342355/1‑93 (MQ=255)
                   tGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGAtc                                                                      >  1:1197373/1‑93 (MQ=255)
                      gTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCa                                                                   <  2:1742308/93‑1 (MQ=255)
                      gTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCa                                                                   >  1:3896459/1‑93 (MQ=255)
                       tGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAg                                                                  >  2:4187174/1‑93 (MQ=255)
                       tGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAg                                                                  <  1:1053648/93‑1 (MQ=255)
                       tGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAg                                                                  <  1:3919014/93‑1 (MQ=255)
                        ggCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGc                                                                 <  1:2895652/93‑1 (MQ=255)
                        ggCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGc                                                                 <  2:4320543/93‑1 (MQ=255)
                        ggCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGc                                                                 >  1:4091405/1‑93 (MQ=255)
                           gATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGAc                                                              >  1:8917/1‑93 (MQ=255)
                             tCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGt                                                            >  1:2666981/1‑93 (MQ=255)
                                ggCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGAt                                                         >  1:1655639/1‑93 (MQ=255)
                                 gCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATg                                                        >  1:1096718/1‑93 (MQ=255)
                                   ctctGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGtt                                                      >  1:1325490/1‑93 (MQ=255)
                                   ctctGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGtt                                                      <  2:3519364/93‑1 (MQ=255)
                                    tctGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGttt                                                     >  1:3449553/1‑93 (MQ=255)
                                     ctGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTg                                                    <  2:2937996/93‑1 (MQ=255)
                                      tGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGt                                                   <  2:1201201/93‑1 (MQ=255)
                                       gCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTc                                                  <  2:3010141/93‑1 (MQ=255)
                                       gCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTc                                                  >  2:4269349/1‑93 (MQ=255)
                                         aCTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTcgc                                                >  1:4150005/1‑93 (MQ=255)
                                            gCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGAc                                             >  1:4535774/1‑93 (MQ=255)
                                               tGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTc                                          <  2:3443134/93‑1 (MQ=255)
                                                gAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCg                                         >  2:4516300/1‑93 (MQ=255)
                                                 aaGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGa                                        >  1:488406/1‑93 (MQ=255)
                                                  aGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGaa                                       >  1:1507996/1‑93 (MQ=255)
                                                   gTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCTACGTGATGTTTGTCGCGACGTCGAAc                                      >  2:3796373/1‑93 (MQ=255)
                                                   gTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAAc                                      >  2:2676176/1‑93 (MQ=255)
                                                     gCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTc                                    <  2:734708/93‑1 (MQ=255)
                                                     gCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTc                                    >  1:2537564/1‑93 (MQ=255)
                                                     gCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTc                                    >  2:1816048/1‑93 (MQ=255)
                                                      cTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTcc                                   <  1:3054516/93‑1 (MQ=255)
                                                       tGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCa                                  <  2:1261922/93‑1 (MQ=255)
                                                       tGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCa                                  >  2:1113648/1‑93 (MQ=255)
                                                       tGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCa                                  >  2:2742266/1‑93 (MQ=255)
                                                        ggATCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCAt                                 >  2:1174772/1‑93 (MQ=255)
                                                           tCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGaa                              <  2:3959591/93‑1 (MQ=255)
                                                           tCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGaa                              <  2:4094886/93‑1 (MQ=255)
                                                           tCCAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGaa                              >  1:3431336/1‑93 (MQ=255)
                                                            ccAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAAc                             <  2:3330112/93‑1 (MQ=255)
                                                             cAGAGCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACa                            <  2:2286099/93‑1 (MQ=255)
                                                              cgagCAGATCGTAGCGTTCAGCGACCACCACGTGGAAGTGGATTACGAGCTCAGCGACGTGATGTTTGTCGCGACGTCGAGCTCCATGAACAt                           <  2:1928929/92‑1 (MQ=255)
                                                               gagCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACAtt                          >  2:4110496/1‑93 (MQ=255)
                                                               gagCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACAtt                          >  2:2305673/1‑93 (MQ=255)
                                                               gagCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACAtt                          >  1:1574890/1‑93 (MQ=255)
                                                                agCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTc                         >  1:986004/1‑93 (MQ=255)
                                                                 gCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGCACTCCATGAACATTcc                        >  2:92164/1‑93 (MQ=255)
                                                                 gCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTcc                        <  1:239527/93‑1 (MQ=255)
                                                                 gCAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTcc                        >  1:392606/1‑93 (MQ=255)
                                                                  cAGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCg                       <  1:1215702/93‑1 (MQ=255)
                                                                   aGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCgg                      <  1:537854/93‑1 (MQ=255)
                                                                   aGAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCgg                      >  1:4436487/1‑93 (MQ=255)
                                                                    gAACGTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGc                     >  1:1628066/1‑93 (MQ=255)
                                                                    gAACGTAGCGTCCAGCGTCCACCACCTGGAAGTTGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGc                     <  1:1006712/93‑1 (MQ=255)
                                                                        gTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCGGCGACGCGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCg                 >  1:2516518/1‑93 (MQ=255)
                                                                        gTAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCg                 <  1:1809628/93‑1 (MQ=255)
                                                                         tAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCgc                >  2:3008274/1‑93 (MQ=255)
                                                                         tAGCGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCgc                >  2:4044962/1‑93 (MQ=255)
                                                                            cGTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCgctgc             <  2:1589329/93‑1 (MQ=255)
                                                                             gTTCAGCGACCACCACCTGGAAGTGGATTTCGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCgctgct            <  1:3080778/93‑1 (MQ=255)
                                                                             gTTCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCgctgct            <  2:32596/93‑1 (MQ=255)
                                                                              ttCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGctgctg           <  1:3227363/93‑1 (MQ=255)
                                                                               tCAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTgg          <  2:1273456/93‑1 (MQ=255)
                                                                                cAGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGa         <  1:4322463/93‑1 (MQ=255)
                                                                                 aGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGAt        >  2:2824808/1‑93 (MQ=255)
                                                                                 aGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGAt        >  1:2648179/1‑93 (MQ=255)
                                                                                 aGCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGAt        <  2:4246117/93‑1 (MQ=255)
                                                                                  gCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATc       <  1:2084611/93‑1 (MQ=255)
                                                                                  gCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATc       >  2:4042955/1‑93 (MQ=255)
                                                                                  gCGACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATc       <  2:3071361/93‑1 (MQ=255)
                                                                                    gACCACCACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCAACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGt     <  1:2568835/93‑1 (MQ=255)
                                                                                     accaccaccTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTa    <  1:581339/93‑1 (MQ=255)
                                                                                     accaccaccTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTa    >  2:2364002/1‑93 (MQ=255)
                                                                                      ccaccaCCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTAt   >  1:1002234/1‑93 (MQ=255)
                                                                                       caccacCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATg  <  2:1582253/93‑1 (MQ=255)
                                                                                          |                                                                                         
GACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATG  >  NC_000913/459387‑459566

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