Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,890,225 C→T R345R (CGG→CGA b1808 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,890,2250CT99.5% 712.6 / NA 197R345R (CGG→CGAb1808hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (108/88);  minor base A (0/1);  total (108/89)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 203 total aligned reads displayed.

GAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGCCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTC  >  NC_000913/1890135‑1890314
                                                                                          |                                                                                         
gAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCg                                                                                         >  1:7205122/1‑93 (MQ=255)
 aGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgt                                                                                        <  1:5737923/93‑1 (MQ=255)
 aGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgt                                                                                        <  2:6568791/93‑1 (MQ=255)
  gcggcgTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgtg                                                                                       >  1:1551755/1‑93 (MQ=255)
  gcggcgTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgtg                                                                                       <  1:2951742/93‑1 (MQ=255)
       gTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGa                                                                                  <  2:2121822/93‑1 (MQ=255)
        tGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGAt                                                                                 <  1:3973013/93‑1 (MQ=255)
        tGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGAt                                                                                 >  2:7451634/1‑93 (MQ=255)
         gagaGCCAGAGTAAAGTGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATa                                                                                >  1:1195678/1‑93 (MQ=255)
         gagaGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATa                                                                                >  2:6093689/1‑93 (MQ=255)
         gagaGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATa                                                                                >  2:7794163/1‑93 (MQ=255)
         gagaGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATa                                                                                >  1:6673497/1‑93 (MQ=255)
          agagCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATAc                                                                               <  2:5846722/93‑1 (MQ=255)
          agagCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATAc                                                                               >  1:765226/1‑93 (MQ=255)
           gagCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACa                                                                              >  1:2608405/1‑93 (MQ=255)
           gagCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACa                                                                              >  2:947548/1‑93 (MQ=255)
            agCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACaa                                                                             <  2:3412966/93‑1 (MQ=255)
             gCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAAc                                                                            >  1:5526122/1‑93 (MQ=255)
             gCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAAc                                                                            <  1:127038/93‑1 (MQ=255)
               cAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGCTACAATGt                                                                          >  2:4688415/1‑93 (MQ=255)
               cAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGt                                                                          <  2:4137520/93‑1 (MQ=255)
                agagTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTg                                                                         <  2:5113639/93‑1 (MQ=255)
                agagTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTg                                                                         <  2:5275909/93‑1 (MQ=255)
                agagTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTg                                                                         <  2:5737731/93‑1 (MQ=255)
                    tAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgcgc                                                                     <  2:6521444/93‑1 (MQ=255)
                        aTGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTc                                                                 <  1:3446033/93‑1 (MQ=255)
                           gcgcGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaaa                                                              >  1:910855/1‑93 (MQ=255)
                           gcgcGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaaa                                                              <  2:5235202/93‑1 (MQ=255)
                           gcgcGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaaa                                                              >  2:1703642/1‑93 (MQ=255)
                           gcgcGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaaa                                                              <  1:6059149/93‑1 (MQ=255)
                           gcccgAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaaa                                                              <  1:6508618/90‑1 (MQ=255)
                               gAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCTAATGcc                                                          >  1:783943/1‑93 (MQ=255)
                               gAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGcc                                                          >  2:2648951/1‑93 (MQ=255)
                                aaGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCg                                                         >  2:1742899/1‑93 (MQ=255)
                                  gTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCa                                                       >  1:7405921/1‑93 (MQ=255)
                                    gCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATc                                                     <  1:4854732/93‑1 (MQ=255)
                                    gCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATc                                                     >  1:5183666/1‑93 (MQ=255)
                                     cATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATcc                                                    <  1:1434481/93‑1 (MQ=255)
                                     cATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATcc                                                    <  1:225393/93‑1 (MQ=255)
                                       ttCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAg                                                  >  1:5223712/1‑93 (MQ=255)
                                       ttCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAg                                                  >  2:2899348/1‑93 (MQ=255)
                                        tCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGc                                                 <  2:4312491/93‑1 (MQ=255)
                                        tCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGc                                                 >  1:7358781/1‑93 (MQ=255)
                                          gTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCaa                                               <  1:3758027/93‑1 (MQ=255)
                                          gTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCaa                                               <  2:257363/93‑1 (MQ=255)
                                          gTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCaa                                               >  2:7673851/1‑93 (MQ=255)
                                          gTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCaa                                               >  1:7810005/1‑93 (MQ=255)
                                          gTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCaa                                               >  2:7544739/1‑93 (MQ=255)
                                           tACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAg                                              >  2:1255366/1‑93 (MQ=255)
                                           tACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAg                                              >  2:1550312/1‑93 (MQ=255)
                                             ccAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGc                                            >  1:2559335/1‑93 (MQ=255)
                                              cAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCg                                           <  2:3816142/93‑1 (MQ=255)
                                              cAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCg                                           <  1:6233161/93‑1 (MQ=255)
                                              cAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCg                                           >  1:3584564/1‑93 (MQ=255)
                                              cAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCg                                           >  2:6752137/1‑93 (MQ=255)
                                               aGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCgg                                          <  1:5598317/93‑1 (MQ=255)
                                                gTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGa                                         >  2:6870500/1‑93 (MQ=255)
                                                gTAGCTGTAACACGGCAGTTTGATCTCTTTTAGCCGCTACGGACGTGTACGTTCCAACGTAGCGCGCTCAAATGCCACATCCCGCCAGGCGGa                                         <  2:6000263/93‑1 (MQ=255)
                                                 tAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGaa                                        <  2:1176268/93‑1 (MQ=255)
                                                 tAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGaa                                        <  2:2722586/93‑1 (MQ=255)
                                                  aGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGaag                                       >  2:1933254/1‑92 (MQ=255)
                                                  aGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAAc                                       <  2:5415613/93‑1 (MQ=255)
                                                   gCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACg                                      >  1:4034989/1‑93 (MQ=255)
                                                    cTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGc                                     >  1:5030826/1‑93 (MQ=255)
                                                     tGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGcc                                    >  1:4703088/1‑93 (MQ=255)
                                                      gTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGccc                                   >  2:304262/1‑93 (MQ=255)
                                                      gTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGccc                                   >  1:3311819/1‑93 (MQ=255)
                                                      gTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGccc                                   >  2:2197559/1‑93 (MQ=255)
                                                         gCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCACATGCCGCATCCAGCAAGGCGGAACGCCCCAg                                >  2:3967991/1‑93 (MQ=255)
                                                         gCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAg                                >  2:705108/1‑93 (MQ=255)
                                                          cccGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGt                               <  1:7587264/93‑1 (MQ=255)
                                                          cccGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGt                               >  1:3175627/1‑93 (MQ=255)
                                                           ccGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTg                              >  1:2825821/1‑93 (MQ=255)
                                                           ccGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTg                              >  1:3782838/1‑93 (MQ=255)
                                                           ccGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGACGCATCCAGCACGGCGGAACGACCCAGTg                              >  2:4076030/1‑93 (MQ=255)
                                                              gCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACa                           >  2:2537952/1‑93 (MQ=255)
                                                               cTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAg                          <  2:803197/93‑1 (MQ=255)
                                                                tgattgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGt                         >  2:5297329/1‑93 (MQ=255)
                                                                tgattgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGt                         >  1:4396104/1‑93 (MQ=255)
                                                                tgattgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGt                         >  2:7224315/1‑93 (MQ=255)
                                                                tgattgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGt                         >  2:7189134/1‑93 (MQ=255)
                                                                 gattgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGtt                        <  1:3517496/93‑1 (MQ=255)
                                                                 gattgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGtt                        >  1:7156299/1‑93 (MQ=255)
                                                                  attgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGttt                       <  2:7320962/93‑1 (MQ=255)
                                                                  attgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGAATCCAGCAAGGCGGAACGCCCCAGTGACAGttt                       <  1:6965252/93‑1 (MQ=255)
                                                                    tgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCg                     >  1:6573549/1‑93 (MQ=255)
                                                                    tgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCg                     >  1:3710807/1‑93 (MQ=255)
                                                                    tgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCg                     >  1:775945/1‑93 (MQ=255)
                                                                     gatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGc                    >  1:7056374/1‑93 (MQ=255)
                                                                       tctcTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCa                  >  2:1614657/1‑93 (MQ=255)
                                                                       tctcTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCa                  >  2:6159543/1‑93 (MQ=255)
                                                                        ctctTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCACGGCGGAACGCCCCAGTGACAGTTTCGCCAc                 >  2:2446494/1‑93 (MQ=255)
                                                                         tctTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACg                >  2:5600705/1‑93 (MQ=255)
                                                                           ttttAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGCGACAGTTTCGCCACGTc              >  1:728159/1‑93 (MQ=255)
                                                                             ttAGCCGCTTCAGTCGTGTGCGGTACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCAt            >  2:3907756/1‑93 (MQ=255)
                                                                             ttAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCAt            <  2:2392215/93‑1 (MQ=255)
                                                                              tAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATa           >  1:888945/1‑93 (MQ=255)
                                                                                gCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAAc         >  2:6424535/1‑93 (MQ=255)
                                                                                   gCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACaaa      >  2:3060532/1‑93 (MQ=255)
                                                                                     ttCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGt    >  1:2479291/1‑93 (MQ=255)
                                                                                          |                                                                                         
GAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGCCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTC  >  NC_000913/1890135‑1890314

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