Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,129,641 C→T V261V (GTG→GTA b2983 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,129,6410CT99.7% 1207.5 / ‑5.8 333V261V (GTG→GTAb2983hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (187/145);  total (187/146)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 352 total aligned reads displayed.

CAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGC  >  NC_000913/3129551‑3129730
                                                                                          |                                                                                         
cAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAatac                                                                                         <  2:5622733/93‑1 (MQ=255)
cAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAatac                                                                                         >  1:2127308/1‑93 (MQ=255)
 aaaCCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGTAGCGGGTCCTAACACGATCCCCACCAATACt                                                                                        >  2:4093088/1‑93 (MQ=255)
 aaaCCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACt                                                                                        <  2:3556222/93‑1 (MQ=255)
  aaCCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACtg                                                                                       <  2:2951848/93‑1 (MQ=255)
  aaCCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACtg                                                                                       >  2:4280693/1‑93 (MQ=255)
  aaCCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACtg                                                                                       <  1:2902244/93‑1 (MQ=255)
     ccGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGtgc                                                                                    <  1:2447394/93‑1 (MQ=255)
     ccGCTGACGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTTtgc                                                                                    >  2:9882537/1‑93 (MQ=255)
      cGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGtgct                                                                                   >  1:2243660/1‑93 (MQ=255)
       gCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGtgctg                                                                                  >  1:5599613/1‑93 (MQ=255)
       gCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGtgctg                                                                                  >  2:1545520/1‑93 (MQ=255)
          ggCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCAc                                                                               >  2:9301039/1‑93 (MQ=255)
            cGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGa                                                                             <  1:6468617/93‑1 (MQ=255)
             gTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAg                                                                            <  2:6641216/93‑1 (MQ=255)
              ttCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGt                                                                           >  2:7273766/1‑93 (MQ=255)
                ccGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTa                                                                         >  2:1371349/1‑93 (MQ=255)
                ccGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTa                                                                         <  2:9089509/93‑1 (MQ=255)
                 cggcggCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTAc                                                                        >  1:6738374/1‑93 (MQ=255)
                  ggcggcGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTAcg                                                                       >  1:4504568/1‑93 (MQ=255)
                  ggcggcGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTAcg                                                                       >  1:9887278/1‑93 (MQ=255)
                     ggcgTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACgcgc                                                                    >  1:6923430/1‑93 (MQ=255)
                     ggcgTCGAAGAATGTCATGGCAATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACgcgc                                                                    <  2:5480690/93‑1 (MQ=255)
                      gcgTCAAAGAATGTCATGGCGATTTTAAATACCCCGGCGGCAGCGGGTCCTAACACGACCCCCACCAATACTGTGCTGCACGAGTTACGCGcc                                                                   <  2:3677165/93‑1 (MQ=255)
                        gTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGa                                                                 <  1:3753048/93‑1 (MQ=255)
                         tCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGAc                                                                <  1:4621161/93‑1 (MQ=255)
                         tCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGAc                                                                >  2:7098340/1‑93 (MQ=255)
                           gaagaaTGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCa                                                              >  1:9525583/1‑93 (MQ=255)
                           gaagaaTGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCa                                                              >  1:7756280/1‑93 (MQ=255)
                           gaagaaTGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCa                                                              >  1:3105085/1‑93 (MQ=255)
                             agaaTGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGa                                                            >  1:3007114/1‑93 (MQ=255)
                             agaaTGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGa                                                            >  2:4926012/1‑93 (MQ=255)
                               aaTGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATg                                                          >  2:3343941/1‑93 (MQ=255)
                                aTGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATgg                                                         <  1:2727937/93‑1 (MQ=255)
                                 tGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGa                                                        <  1:7159554/93‑1 (MQ=255)
                                 tGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTACACGAGTTATGCGACGACCAGATGGa                                                        >  1:4796528/1‑93 (MQ=255)
                                  gTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAg                                                       <  1:9036775/93‑1 (MQ=255)
                                  gTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAg                                                       >  1:7404936/1‑93 (MQ=255)
                                  gTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAg                                                       >  1:7274342/1‑93 (MQ=255)
                                  gTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAg                                                       <  1:3996911/93‑1 (MQ=255)
                                   tCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGt                                                      <  2:2823944/93‑1 (MQ=255)
                                   tCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGt                                                      <  1:8592841/93‑1 (MQ=255)
                                   tCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGt                                                      <  1:9433816/93‑1 (MQ=255)
                                    cATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTg                                                     >  2:4640195/1‑93 (MQ=255)
                                      tGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTggg                                                   >  2:7392788/1‑93 (MQ=255)
                                      tGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTggg                                                   >  1:1406170/1‑93 (MQ=255)
                                       ggagATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGc                                                  <  1:6020507/90‑1 (MQ=255)
                                          gATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAAt                                               >  1:5511880/1‑93 (MQ=255)
                                           aTTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATg                                              >  2:4813584/1‑93 (MQ=255)
                                                aaaTAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATgttggt                                         <  2:7512270/93‑1 (MQ=255)
                                                aaaTAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCCGCACGAGTTACGCGACGACCAGATGGAGTGGGCAATgttggt                                         >  1:8475395/1‑93 (MQ=255)
                                                 aaTAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATgttggtt                                        >  2:3431750/1‑93 (MQ=255)
                                                 aaTAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATgttggtt                                        >  1:6232296/1‑93 (MQ=255)
                                                 aaTAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGCGGGCAATgttggtt                                        >  2:5324148/1‑93 (MQ=255)
                                                  aTAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATgttggttg                                       <  2:789136/93‑1 (MQ=255)
                                                  aTAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATgttggttg                                       >  2:3698576/1‑93 (MQ=255)
                                                    aaCCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGAc                                     <  2:248985/93‑1 (MQ=255)
                                                     aCCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGAcc                                    <  1:4086082/93‑1 (MQ=255)
                                                      cccGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACGGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCa                                   <  1:1814512/93‑1 (MQ=255)
                                                        cggcggCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGa                                 >  2:2626422/1‑93 (MQ=255)
                                                         ggcggcAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGAc                                >  1:7565251/1‑93 (MQ=255)
                                                         ggcggcAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGAc                                <  2:2395177/93‑1 (MQ=255)
                                                           gggcAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAGTGTTGGTTGACCAGACaa                              <  1:9283801/92‑1 (MQ=255)
                                                             gcAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACaaac                            >  1:8881733/1‑92 (MQ=255)
                                                             gcAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACaaaa                            >  1:7912172/1‑93 (MQ=255)
                                                             gcAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACaaaa                            <  1:1348811/93‑1 (MQ=255)
                                                              cAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAAc                           <  2:9003562/93‑1 (MQ=255)
                                                               aGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACt                          >  1:3806005/1‑93 (MQ=255)
                                                                gCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTc                         >  2:1769322/1‑93 (MQ=255)
                                                                gCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTc                         >  2:6329529/1‑93 (MQ=255)
                                                                  gggTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCa                       <  1:1976125/93‑1 (MQ=255)
                                                                   ggcccTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCAc                      <  2:1863093/90‑1 (MQ=255)
                                                                   ggTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCAc                      >  2:3896833/1‑93 (MQ=255)
                                                                   ggTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCAc                      >  1:3399012/1‑93 (MQ=255)
                                                                    gTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCAcg                     >  1:1078750/1‑93 (MQ=255)
                                                                    gTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCAcg                     <  1:6995746/93‑1 (MQ=255)
                                                                     tCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCAcgc                    >  1:655062/1‑93 (MQ=255)
                                                                       cTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACgcgc                  <  2:6989226/93‑1 (MQ=255)
                                                                        tAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGcc                 >  2:1169375/1‑93 (MQ=255)
                                                                         aaCACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCt                <  2:8288831/93‑1 (MQ=255)
                                                                         aaCACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCt                <  2:10035801/93‑1 (MQ=255)
                                                                          acacGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCtt               >  2:8444269/1‑93 (MQ=255)
                                                                          acacGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCtt               >  2:6970073/1‑93 (MQ=255)
                                                                           cacGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCttt              <  1:8171720/93‑1 (MQ=255)
                                                                              gATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAAt           <  2:8049130/93‑1 (MQ=255)
                                                                                tCCCCACCAATACTGTGCGTCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGt         <  1:6150597/93‑1 (MQ=255)
                                                                                 ccccACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTa        <  1:2360734/93‑1 (MQ=255)
                                                                                 ccccACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTa        <  2:9515946/93‑1 (MQ=255)
                                                                                 ccccACCAATACTGTGCTGCACGAGTTACGAGCCGACCAGATGGAGTGGGAAATGTTGGTTGACCAGACAAATCTCCACGCGCCTTTAATGTa        <  1:4314530/93‑1 (MQ=255)
                                                                                  cccACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTAt       >  2:5932611/1‑93 (MQ=255)
                                                                                  cccACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTAt       <  2:783876/93‑1 (MQ=255)
                                                                                   ccaccaATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATc      >  1:1001232/1‑93 (MQ=255)
                                                                                   ccaccaATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATc      <  1:5767051/93‑1 (MQ=255)
                                                                                   ccaccaATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATc      <  1:9784314/93‑1 (MQ=255)
                                                                                     cccaATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCgg    <  2:2203852/92‑1 (MQ=255)
                                                                                     accaATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCgg    >  2:81043/1‑93 (MQ=255)
                                                                                     accaATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGAATCgg    >  2:3917765/1‑93 (MQ=255)
                                                                                      ccaATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCggg   >  1:813496/1‑93 (MQ=255)
                                                                                       caATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGc  <  1:505849/93‑1 (MQ=255)
                                                                                       caATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGc  >  2:2013511/1‑93 (MQ=255)
                                                                                          |                                                                                         
CAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGC  >  NC_000913/3129551‑3129730

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