Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,085,695 T→C K581R (AAG→AGG)  frvR ← putative frv operon regulatory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,085,6950TC99.8% 1674.3 / ‑5.6 515K581R (AAG→AGG) frvRputative frv operon regulatory protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (278/236);  total (278/237)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 540 total aligned reads displayed.

CGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCGATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTA  >  NC_000913/4085604‑4085784
                                                                                           |                                                                                         
cGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCtg                                                                                         <  1:7793717/93‑1 (MQ=255)
  aaaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCtgtg                                                                                       >  1:8734017/1‑93 (MQ=255)
   aaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGa                                                                                      >  2:9322506/1‑93 (MQ=255)
   aaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGa                                                                                      >  2:6879817/1‑93 (MQ=255)
    aTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGaa                                                                                     <  1:6596186/93‑1 (MQ=255)
     tGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAg                                                                                    <  1:5612263/93‑1 (MQ=255)
     tGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAg                                                                                    >  1:7603099/1‑93 (MQ=255)
     tGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAg                                                                                    >  1:10031584/1‑93 (MQ=255)
       cttGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTa                                                                                  <  1:1754949/92‑1 (MQ=255)
        ttGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGACGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTaa                                                                                 >  2:1149547/1‑93 (MQ=255)
          gatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACt                                                                               >  2:3308718/1‑93 (MQ=255)
           atgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTc                                                                              <  2:2986722/93‑1 (MQ=255)
           atgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTc                                                                              >  1:107885/1‑93 (MQ=255)
             gaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCa                                                                            >  2:6396545/1‑93 (MQ=255)
              aCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCAt                                                                           >  1:4020148/1‑93 (MQ=255)
                tataATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATaa                                                                         >  1:8200373/1‑93 (MQ=255)
                tataATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATaa                                                                         <  1:8077699/93‑1 (MQ=255)
                 ataataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATaaa                                                                        >  2:1264466/1‑93 (MQ=255)
                 ataataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATaaa                                                                        <  2:5222372/93‑1 (MQ=255)
                 ataataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATaaa                                                                        <  1:5985540/93‑1 (MQ=255)
                 ataataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATaaa                                                                        >  1:4122352/1‑93 (MQ=255)
                 ataataGCCTGCGCTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATaaa                                                                        >  2:8683260/1‑93 (MQ=255)
                  taataGCCTGCTTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCACCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAAt                                                                       <  2:783798/93‑1 (MQ=255)
                  taataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAAt                                                                       <  2:6660457/93‑1 (MQ=255)
                    ataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACGGTTCACCCCCTGTGAAGTAACTCCATAAATGc                                                                     <  2:981697/93‑1 (MQ=255)
                    ataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGc                                                                     >  2:3749112/1‑93 (MQ=255)
                    ataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGc                                                                     >  1:8500668/1‑93 (MQ=255)
                     taGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGcc                                                                    >  1:4652956/1‑93 (MQ=255)
                     taGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGcc                                                                    <  2:9158630/93‑1 (MQ=255)
                      aGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCt                                                                   >  2:2277635/1‑93 (MQ=255)
                       gCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTc                                                                  >  1:9983325/1‑93 (MQ=255)
                       gCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTc                                                                  >  1:4238201/1‑93 (MQ=255)
                        ccTGCGTTAAACGTTTTTGATGATAAATGCAGTGTGGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCa                                                                 <  1:2204647/93‑1 (MQ=255)
                          tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCata                                                               <  1:8340996/93‑1 (MQ=255)
                          tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCata                                                               >  2:1496628/1‑93 (MQ=255)
                          tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCata                                                               <  2:3939042/93‑1 (MQ=255)
                          tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCata                                                               >  2:5855711/1‑93 (MQ=255)
                          tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCata                                                               >  2:5998311/1‑93 (MQ=255)
                          tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCata                                                               >  2:7312986/1‑93 (MQ=255)
                          tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCata                                                               >  2:8297361/1‑93 (MQ=255)
                          tGCGTTAAACGTTTTTGATAATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCata                                                               <  1:7910262/93‑1 (MQ=255)
                            acttAAAAGTTTTTGATGAAAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATc                                                             <  1:6961215/91‑1 (MQ=255)
                             gTTAAACGTTTTTGATGATATTTGCATTGTCGGATGCGGCGCATCCGACTTAGTTCACCCCCTTTGTAGTAACTCCATAAACGTTTCATATcc                                                            >  1:8953220/1‑93 (MQ=255)
                             gTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATcc                                                            <  2:4417001/93‑1 (MQ=255)
                             gTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATcc                                                            <  2:7725119/93‑1 (MQ=255)
                              ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCt                                                           >  2:5769299/1‑93 (MQ=255)
                                aaaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGt                                                         >  2:2600487/1‑93 (MQ=255)
                                 aaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGtt                                                        <  2:8556666/93‑1 (MQ=255)
                                  aCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTa                                                       >  2:352310/1‑93 (MQ=255)
                                   catttttGATGATAATTGCAGTGTCGGATGAAGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCACATATCCTGTTaa                                                      <  2:8924380/91‑1 (MQ=255)
                                   cGTTTTGGATGATAAATGCGGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGGGAAGTAACTCCATAAATGCCTCATATCCTGTTaa                                                      <  1:4472679/93‑1 (MQ=255)
                                          gatgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCa                                               >  2:3056646/1‑93 (MQ=255)
                                          gatgatAAATGCAGTGTCGGATGCGGCGGATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACTCTGCa                                               >  1:9428708/1‑93 (MQ=255)
                                          gatgatAAATGCAGTGTCGGATGCGGCCCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCa                                               <  1:3692219/93‑1 (MQ=255)
                                            tgatAAATGCGGTGTCGGATGCGGCGCACCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaat                                             <  1:2094954/93‑1 (MQ=255)
                                            tgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaat                                             >  1:6667250/1‑93 (MQ=255)
                                              atAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataa                                           >  1:3338180/1‑93 (MQ=255)
                                               tAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataat                                          <  1:2815458/93‑1 (MQ=255)
                                                 aaTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtct                                        >  2:4699481/1‑93 (MQ=255)
                                                   tGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTcc                                      <  2:5128077/93‑1 (MQ=255)
                                                   tGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTACATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTcc                                      >  1:2371874/1‑93 (MQ=255)
                                                    gCAGTGTCGGATGCGGCGCATCCGACACAGTTCCCCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccg                                     <  1:5143518/93‑1 (MQ=255)
                                                    gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccg                                     >  2:8414536/1‑93 (MQ=255)
                                                    gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccg                                     >  1:2616846/1‑93 (MQ=255)
                                                    gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccc                                     >  2:7502305/1‑92 (MQ=255)
                                                      aGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgcc                                   >  1:9284080/1‑93 (MQ=255)
                                                      aGTGTCGGAGGCGCCCCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgcc                                   <  2:5441859/93‑1 (MQ=255)
                                                       gtgtCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgccg                                  >  1:4776513/1‑93 (MQ=255)
                                                       gtgtCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgccg                                  <  2:452784/93‑1 (MQ=255)
                                                        tgtCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCgg                                 >  2:3628000/1‑93 (MQ=255)
                                                         gtCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGa                                <  2:4381961/93‑1 (MQ=255)
                                                         gtCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGa                                <  2:626848/93‑1 (MQ=255)
                                                          tCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGAt                               >  1:1372278/1‑93 (MQ=255)
                                                            ggATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGc                             <  1:3242958/93‑1 (MQ=255)
                                                            ggATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGc                             <  2:9618819/93‑1 (MQ=255)
                                                             gATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt                            >  1:7189624/1‑93 (MQ=255)
                                                             gATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt                            <  2:4224161/93‑1 (MQ=255)
                                                              aTGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTg                           >  2:6518477/1‑93 (MQ=255)
                                                               tGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGa                          <  1:4292832/93‑1 (MQ=255)
                                                               tGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGa                          <  2:8160541/93‑1 (MQ=255)
                                                                 cggcgCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGaca                        >  2:5599956/1‑93 (MQ=255)
                                                                       cacccGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAata                   <  2:1873186/90‑1 (MQ=255)
                                                                       caTCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAata                   <  2:4906937/92‑1 (MQ=255)
                                                                       caTCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAata                   <  1:9646182/92‑1 (MQ=255)
                                                                        aTCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAatat                  >  1:4921218/1‑93 (MQ=255)
                                                                         tCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATg                 <  2:6729612/93‑1 (MQ=255)
                                                                         tCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATg                 <  1:3084294/93‑1 (MQ=255)
                                                                         tCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATg                 <  1:2957460/93‑1 (MQ=255)
                                                                           cGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCt               <  1:6354347/93‑1 (MQ=255)
                                                                           cGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCt               <  1:1639642/93‑1 (MQ=255)
                                                                            gACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCtt              >  2:4036335/1‑93 (MQ=255)
                                                                            gACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCtt              >  2:8759674/1‑93 (MQ=255)
                                                                            gACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCtt              >  1:5539392/1‑93 (MQ=255)
                                                                               acaGTTCCCCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGcc           <  2:2721755/93‑1 (MQ=255)
                                                                               acaGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGcc           <  2:7141505/93‑1 (MQ=255)
                                                                               acaGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGcc           >  1:9822723/1‑93 (MQ=255)
                                                                               acaGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGcc           >  1:6249418/1‑93 (MQ=255)
                                                                                   ttCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGAt       <  1:6801592/93‑1 (MQ=255)
                                                                                      cccccCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGc    <  2:2120105/92‑1 (MQ=255)
                                                                                      cccccCTGCGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGc    <  1:3531055/92‑1 (MQ=255)
                                                                                           |                                                                                         
CGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCGATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTA  >  NC_000913/4085604‑4085784

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