Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,510,339 T→C V221A (GTG→GCG)  b1440 → putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,510,3390TC99.4% 548.3 / NA 177V221A (GTG→GCG) b1440putative transport protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (78/98);  minor base A (0/1);  total (78/99)
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 188 total aligned reads displayed.

TTGCGGACGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGTGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGAGGGT  >  NC_000913/1510249‑1510427
                                                                                          |                                                                                        
ttGCGGACGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGgcgc                                                                                        >  1:2048857/1‑93 (MQ=255)
 tGCGGACGCTGCTTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCt                                                                                       <  2:5963048/93‑1 (MQ=255)
    ggACGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGaa                                                                                    >  2:5099589/1‑93 (MQ=255)
     gACGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGaaa                                                                                   >  1:2212770/1‑93 (MQ=255)
     gACGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGaaa                                                                                   >  1:4299145/1‑93 (MQ=255)
         cTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGtgc                                                                               >  2:601689/1‑93 (MQ=255)
          tGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGtgct                                                                              >  2:4914860/1‑93 (MQ=255)
          tGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTCAAAGggcg                                                                              >  1:1575977/1‑92 (MQ=255)
           gCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGtgctg                                                                             >  1:5531973/1‑93 (MQ=255)
           gCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGtgctg                                                                             >  1:4810313/1‑93 (MQ=255)
             gttgttCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTgcg                                                                           <  1:991737/93‑1 (MQ=255)
              ttgttCGCTAAAGCCGCGCAGCCGCGTTTGCGCTTCAGTGATCCGTAGCAACTCACCGAAGAACAGTACCAGGCGGCCCTGAAAGTGCTgcgc                                                                          <  1:7660157/93‑1 (MQ=255)
                 ttCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCt                                                                       >  2:607981/1‑93 (MQ=255)
                    gTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCaa                                                                    >  1:3344583/1‑93 (MQ=255)
                     ttAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAac                                                                   >  2:2569993/1‑93 (MQ=255)
                     ttAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAac                                                                   >  1:6385210/1‑93 (MQ=255)
                      tAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAaca                                                                  <  1:1890969/93‑1 (MQ=255)
                      tAAAGCCACTCAGCCCCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAaca                                                                  <  1:4476910/93‑1 (MQ=255)
                       aaaGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAacac                                                                 >  1:5824981/1‑93 (MQ=255)
                       aaaGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAacac                                                                 >  1:3812560/1‑93 (MQ=255)
                         aGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAg                                                               <  2:1465725/93‑1 (MQ=255)
                          gCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGt                                                              <  1:1966807/93‑1 (MQ=255)
                             aCTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTg                                                           <  2:6599878/93‑1 (MQ=255)
                             aCTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTg                                                           >  2:3873500/1‑93 (MQ=255)
                              cTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGa                                                          <  2:3774582/93‑1 (MQ=255)
                               tCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGac                                                         >  1:6011405/1‑92 (MQ=255)
                               tCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGAt                                                         <  1:3704034/93‑1 (MQ=255)
                               tCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGAt                                                         >  1:7905403/1‑93 (MQ=255)
                               tCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGAt                                                         >  2:6745079/1‑93 (MQ=255)
                                  gCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCa                                                      >  2:4804505/1‑93 (MQ=255)
                                   ccGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCAt                                                     <  1:575523/93‑1 (MQ=255)
                                    cGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATc                                                    <  2:1280617/93‑1 (MQ=255)
                                    cGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATc                                                    <  1:2918645/93‑1 (MQ=255)
                                    cGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATc                                                    <  2:3509591/93‑1 (MQ=255)
                                    cGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATc                                                    <  2:6292086/93‑1 (MQ=255)
                                    cGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATc                                                    <  1:3543611/93‑1 (MQ=255)
                                        gTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTa                                                >  1:6196418/1‑93 (MQ=255)
                                        gTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTa                                                >  2:7893327/1‑93 (MQ=255)
                                         ttGGGGATGAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTAc                                               <  1:2533933/93‑1 (MQ=255)
                                         ttGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTAc                                               <  2:5862018/93‑1 (MQ=255)
                                         ttGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGGACCGGGCGGCGCTGAAAGTACTGCGCGCTCAACACAGTTTGATCCATCGCTAc                                               <  2:5862020/93‑1 (MQ=255)
                                            ggCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTgg                                            <  2:3473238/93‑1 (MQ=255)
                                               aTCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCAt                                         <  2:5424884/93‑1 (MQ=255)
                                               aTCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCAt                                         <  1:2492279/93‑1 (MQ=255)
                                                tCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCTGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGTCAta                                        >  2:3975453/1‑92 (MQ=255)
                                                tCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATg                                        >  1:4648914/1‑93 (MQ=255)
                                                tCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATg                                        >  1:1711025/1‑93 (MQ=255)
                                                 cAGCGTTCCGTATCACCACACCGAAGAACAGTACCAGGCGGCGATGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGa                                       <  2:7047985/93‑1 (MQ=255)
                                                 cAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGa                                       >  1:5151840/1‑93 (MQ=255)
                                                 cAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGa                                       <  1:5798074/93‑1 (MQ=255)
                                                   gCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGaca                                     >  2:725856/1‑93 (MQ=255)
                                                   gCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGaca                                     <  1:5684391/93‑1 (MQ=255)
                                                     gATCCGTATCAACTCGCCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACt                                   <  2:4554876/93‑1 (MQ=255)
                                                     gATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACt                                   >  2:7149637/1‑93 (MQ=255)
                                                      aTCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTa                                  <  2:4885503/93‑1 (MQ=255)
                                                       tCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTAc                                 <  2:2801040/93‑1 (MQ=255)
                                                        ccGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTAcc                                <  1:3180600/93‑1 (MQ=255)
                                                         cGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCg                               <  1:2141347/93‑1 (MQ=255)
                                                         cGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCg                               >  1:6982571/1‑93 (MQ=255)
                                                         cGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCg                               <  1:6300208/93‑1 (MQ=255)
                                                          gTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGt                              <  2:1707350/93‑1 (MQ=255)
                                                          gTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCAACGCTATTGGCATGACACTACCGt                              >  1:7046505/1‑93 (MQ=255)
                                                           aaTCAACACACCAAAGAACAGTACCAGGCGGAACTGAAAAAGAAGCGCACTCAAAACAGATAGATCCATAAATACTGACATAACACTACCGTg                             <  1:7022560/92‑1 (MQ=255)
                                                            aTCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGc                            >  2:262578/1‑93 (MQ=255)
                                                            aTCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGc                            >  2:3232525/1‑93 (MQ=255)
                                                            aTCAACTCACCGAAGAACAGTACCAGGCGGAGATGAAAGTGCTGCGCGCTCAACACCGTTTGATCCATCGCTACTGGCATGACACTACCGTGc                            >  2:6450677/1‑93 (MQ=255)
                                                             tCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCa                           <  1:4719017/93‑1 (MQ=255)
                                                             tCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCa                           >  2:2933460/1‑93 (MQ=255)
                                                             tCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCa                           <  2:2548480/93‑1 (MQ=255)
                                                             tCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCa                           >  1:7925953/1‑93 (MQ=255)
                                                              cAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCaa                          <  1:5465650/93‑1 (MQ=255)
                                                               aaCTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCaaa                         >  2:571599/1‑93 (MQ=255)
                                                                aCTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAAt                        >  2:1926238/1‑93 (MQ=255)
                                                                 cTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATg                       >  1:5678846/1‑93 (MQ=255)
                                                                  tCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGa                      <  2:431619/93‑1 (MQ=255)
                                                                   cACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAg                     <  1:7725224/93‑1 (MQ=255)
                                                                   cACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAg                     <  1:3078231/93‑1 (MQ=255)
                                                                    aCCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGc                    <  2:843591/93‑1 (MQ=255)
                                                                       gaagaaCAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGAt                 >  2:7545091/1‑93 (MQ=255)
                                                                       gaagaaCAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGAt                 <  1:1936806/93‑1 (MQ=255)
                                                                       gaagaaCAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGAt                 <  2:2096246/93‑1 (MQ=255)
                                                                         agaaCAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGAttt               <  1:1111053/93‑1 (MQ=255)
                                                                          gaaCAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTc              <  1:5053407/93‑1 (MQ=255)
                                                                          gaaCAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTc              >  1:2989695/1‑93 (MQ=255)
                                                                            aCAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCaa            >  2:5440589/1‑93 (MQ=255)
                                                                            aCAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCaa            >  1:904780/1‑93 (MQ=255)
                                                                             cAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCaaa           <  1:3023962/93‑1 (MQ=255)
                                                                              aGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCaaaa          >  1:2113296/1‑93 (MQ=255)
                                                                              aGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCaaaa          <  1:5120595/93‑1 (MQ=255)
                                                                               gTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCaaaaa         >  1:6503169/1‑93 (MQ=255)
                                                                               gTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCaaaaa         >  1:6075834/1‑93 (MQ=255)
                                                                                tACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAAc        <  2:4313115/93‑1 (MQ=255)
                                                                                  ccAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTAGTGGCATGACACTACCTTGCAAATGAGCGATTTCAAAAACGa      <  2:5069839/93‑1 (MQ=255)
                                                                                  ccAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGa      >  2:3094034/1‑93 (MQ=255)
                                                                                   cAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGAg     <  2:6908096/93‑1 (MQ=255)
                                                                                   cAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGAg     >  1:1070507/1‑93 (MQ=255)
                                                                                    aGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGAgg    <  2:6656256/93‑1 (MQ=255)
                                                                                    aGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGAgg    >  1:4177778/1‑93 (MQ=255)
                                                                                    aGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGAgg    >  2:5574123/1‑93 (MQ=255)
                                                                                      gcggcgCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGAGGgt  <  2:2527927/93‑1 (MQ=255)
                                                                                          |                                                                                        
TTGCGGACGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGTGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGAGGGT  >  NC_000913/1510249‑1510427

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: