Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,566,541 T→C I103V (ATA→GTA)  b1491 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,566,5410TC99.5% 664.3 / ‑7.0 206I103V (ATA→GTA) b1491hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (95/110);  total (96/110)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 216 total aligned reads displayed.

TCTTACCGGTCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTATGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGAGG  >  NC_000913/1566451‑1566630
                                                                                          |                                                                                         
tCTTACCGGTCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGc                                                                                         >  2:1793290/1‑93 (MQ=255)
 cTTACCGGTCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGcc                                                                                        <  1:1065415/93‑1 (MQ=255)
   tACCGGTCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCga                                                                                      <  1:4646587/93‑1 (MQ=255)
      cGGTCATAAGCTCGGACCACGGCATACTTTTCGATGGCCACAGGGCGGTCCCGTCCGGCTTGACCTGGACAAAGACCGTGTTTACGCCGAGAc                                                                                   <  1:6299512/93‑1 (MQ=255)
       ggTCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACg                                                                                  >  2:5194886/1‑93 (MQ=255)
        gTCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGt                                                                                 <  1:5637797/93‑1 (MQ=255)
        gTCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGt                                                                                 <  2:4479/93‑1 (MQ=255)
         tCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGtt                                                                                <  2:1558993/93‑1 (MQ=255)
         tCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGtt                                                                                >  1:1896664/1‑93 (MQ=255)
          cATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTg                                                                               <  2:4211750/93‑1 (MQ=255)
          cATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGCACCGTCCGGCTTGAACTGGAAAAAGACCGTGTTTACGCCGAGACGTTg                                                                               >  2:4646919/1‑93 (MQ=255)
           aTAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACGGTGTTTACGCCGAGACGTTGc                                                                              >  1:2365807/1‑93 (MQ=255)
            tAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCa                                                                             >  2:1194080/1‑93 (MQ=255)
            tAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCa                                                                             >  1:579290/1‑93 (MQ=255)
              aGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAga                                                                           >  1:3512496/1‑93 (MQ=255)
              aGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAga                                                                           <  1:2959989/93‑1 (MQ=255)
                aTCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAgatg                                                                         <  1:1139538/93‑1 (MQ=255)
                aTCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAgatg                                                                         >  1:1817050/1‑93 (MQ=255)
                aTCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAgatg                                                                         >  2:4557284/1‑93 (MQ=255)
                 tCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAgatga                                                                        <  1:6261664/93‑1 (MQ=255)
                 tCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAgatga                                                                        <  1:1334444/93‑1 (MQ=255)
                  cGGACCACGGCAAAATTTTCGATGTCCACACGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAgatgat                                                                       >  2:2617613/1‑93 (MQ=255)
                  cGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAgatgat                                                                       >  2:3068984/1‑93 (MQ=255)
                  cGGACAACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGCAAAAGACGGTGTTTACGCCGAGAAGTTGGAgatgat                                                                       >  2:3991015/1‑93 (MQ=255)
                   ggACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATc                                                                      <  2:4688134/93‑1 (MQ=255)
                   ggACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATc                                                                      >  1:992141/1‑93 (MQ=255)
                    gACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACATGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATcc                                                                     <  1:4447159/93‑1 (MQ=255)
                     aCCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCa                                                                    <  1:3285863/93‑1 (MQ=255)
                        aCGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATTCAGtt                                                                 >  2:5704191/1‑93 (MQ=255)
                        aCGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGtt                                                                 <  2:6211505/93‑1 (MQ=255)
                           gCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTTGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGt                                                              >  2:3526169/1‑93 (MQ=255)
                             aaaaTTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCg                                                            <  2:2044814/93‑1 (MQ=255)
                              aaaTTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGa                                                           <  1:4351792/93‑1 (MQ=255)
                              aaaTTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGa                                                           >  1:1405696/1‑93 (MQ=255)
                              aaaTTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGa                                                           >  1:390120/1‑93 (MQ=255)
                               aaTTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGat                                                          >  2:303929/1‑93 (MQ=255)
                               aaTTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGat                                                          >  1:1197725/1‑93 (MQ=255)
                                 ttttCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGatca                                                        <  2:847537/93‑1 (MQ=255)
                                 ttttCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGatca                                                        <  1:5045129/93‑1 (MQ=255)
                                  tttCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGatcat                                                       <  2:2654468/93‑1 (MQ=255)
                                  tttCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGatcat                                                       >  1:3677844/1‑93 (MQ=255)
                                  tttCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGatcat                                                       <  1:2063597/93‑1 (MQ=255)
                                      gATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGcc                                                   <  1:4836141/93‑1 (MQ=255)
                                       aTGGCCACAGGGCGTTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCt                                                  <  1:2237343/93‑1 (MQ=255)
                                       aTGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCt                                                  >  1:5117843/1‑93 (MQ=255)
                                         ggCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCtgt                                                >  1:4181523/1‑93 (MQ=255)
                                            cacaGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCtgttgt                                             <  2:4506146/93‑1 (MQ=255)
                                             acaGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCtgttgtt                                            >  2:3695175/1‑93 (MQ=255)
                                              caGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCtgttgttt                                           >  2:5998263/1‑92 (MQ=255)
                                               aGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCtgttgttgt                                          <  1:6019626/93‑1 (MQ=255)
                                               aGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCtgttgttgt                                          <  2:4433583/93‑1 (MQ=255)
                                                gggCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTa                                         <  2:4849876/93‑1 (MQ=255)
                                                 ggCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTac                                        >  1:39558/1‑93 (MQ=255)
                                                   cGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTacac                                      <  2:644323/93‑1 (MQ=255)
                                                    ggTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACg                                     <  1:921327/93‑1 (MQ=255)
                                                    ggTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACg                                     >  1:2174647/1‑93 (MQ=255)
                                                     gaaCCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGCCGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACgg                                    >  2:4192703/3‑93 (MQ=255)
                                                     gTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACgg                                    <  1:2589352/93‑1 (MQ=255)
                                                     gTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACgg                                    >  1:4871981/1‑93 (MQ=255)
                                                     gTACCGTCCGGCTTGACCTGGAAAAAAACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACgg                                    >  2:5003019/1‑93 (MQ=255)
                                                       aCCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGc                                  >  2:1860306/1‑93 (MQ=255)
                                                           tCCGGCTTGGCCTGGTAAAAGACCGTGTTTACGCCCAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGATCTACACGGGCACgg                              >  2:489572/1‑93 (MQ=255)
                                                            ccGGCTTGAGCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGGCTGTTGTTGTACACGGGCCCGGc                             <  1:5254567/93‑1 (MQ=255)
                                                              ggCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCtg                           >  2:5079888/1‑93 (MQ=255)
                                                              ggCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCtg                           >  1:827396/1‑93 (MQ=255)
                                                               gCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCtgg                          >  2:2548493/1‑93 (MQ=255)
                                                               gCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCtgg                          >  2:2217741/1‑93 (MQ=255)
                                                               gCTTGACCTCGAAAAAGACCGTGTTTACGACGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCtgg                          <  2:5865730/93‑1 (MQ=255)
                                                                cTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCtggt                         >  2:226418/1‑93 (MQ=255)
                                                                 ttGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCtggtg                        >  1:2095111/1‑93 (MQ=255)
                                                                  tGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCtggtgg                       <  1:3435141/93‑1 (MQ=255)
                                                                  tGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCtggtgg                       <  1:5555384/93‑1 (MQ=255)
                                                                  tGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCtggtgg                       <  1:2458929/93‑1 (MQ=255)
                                                                   gACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTggg                      <  1:6205093/93‑1 (MQ=255)
                                                                    aCCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTgggg                     >  2:736030/1‑93 (MQ=255)
                                                                    aCCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTgggg                     <  2:1135292/93‑1 (MQ=255)
                                                                     ccTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGt                    <  1:3620694/93‑1 (MQ=255)
                                                                     ccTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGt                    <  2:5420349/93‑1 (MQ=255)
                                                                     ccTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGt                    <  1:725580/93‑1 (MQ=255)
                                                                     ccTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGt                    >  1:3047087/1‑93 (MQ=255)
                                                                     ccTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGt                    <  1:3092478/93‑1 (MQ=255)
                                                                       tGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTa                  <  2:149418/93‑1 (MQ=255)
                                                                         gAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACt                >  2:3801013/1‑93 (MQ=255)
                                                                          aaaaaGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTa               >  1:2083984/1‑93 (MQ=255)
                                                                           aaaaGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTaa              >  2:4251598/1‑93 (MQ=255)
                                                                            aaaGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAAt             <  2:5125418/93‑1 (MQ=255)
                                                                            aaaGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAAt             <  2:4932625/93‑1 (MQ=255)
                                                                             aaGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATg            <  1:4682729/93‑1 (MQ=255)
                                                                             aaGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATg            >  2:824604/1‑93 (MQ=255)
                                                                              aGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGt           >  2:5554354/1‑93 (MQ=255)
                                                                              aGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGt           >  1:5868219/1‑93 (MQ=255)
                                                                               gACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGtt          >  2:4356427/1‑93 (MQ=255)
                                                                                aCCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTa         >  1:4779013/1‑93 (MQ=255)
                                                                                 ccGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTaa        >  1:691856/1‑93 (MQ=255)
                                                                                 ccGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTaa        <  1:4015979/93‑1 (MQ=255)
                                                                                  cGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAAc       <  1:1517958/93‑1 (MQ=255)
                                                                                   gtgtTTACTCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAAcc      <  2:1660279/93‑1 (MQ=255)
                                                                                   gtgtTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAAcc      >  1:4553470/1‑93 (MQ=255)
                                                                                    tgtTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAACCg     <  2:5393339/93‑1 (MQ=255)
                                                                                      tttACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGAg   >  1:4066303/1‑93 (MQ=255)
                                                                                          |                                                                                         
TCTTACCGGTCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTATGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGAGG  >  NC_000913/1566451‑1566630

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