Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,887,809 A→G intergenic (+35/‑56) b3708 → / → tnaB tryptophanase/low affinity tryptophan permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,887,8090AG99.7% 1290.0 / NA 375intergenic (+35/‑56)b3708/tnaBtryptophanase/low affinity tryptophan permease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (187/186);  minor base C (1/0);  total (189/186)
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 380 total aligned reads displayed.

TTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCACTCTGCAT  >  NC_000913/3887719‑3887898
                                                                                          |                                                                                         
 tacgtacgAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGctc                                                                                        >  2:2874423/1‑93 (MQ=255)
  acgtacgAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGctct                                                                                       >  1:1546432/1‑93 (MQ=255)
   cgtacgAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGctctc                                                                                      <  1:822445/93‑1 (MQ=255)
    gtacgAACCGAAAGTATTGCGTCACTTCACCGCAATGCTTAAAGAAGTTTAATTAATTCTACAGAGTGGCTATAAGGATGTTAGCCGctctct                                                                                     <  2:4056882/93‑1 (MQ=255)
     tacgAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCtt                                                                                    <  1:3345407/93‑1 (MQ=255)
     tacgAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCtt                                                                                    >  2:7637344/1‑93 (MQ=255)
     tacgAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCtt                                                                                    >  1:6819775/1‑93 (MQ=255)
     tacgAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCtt                                                                                    <  2:6772008/93‑1 (MQ=255)
         aaCCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCTCTCTCTTAccc                                                                                >  2:3038530/1‑93 (MQ=255)
         aaCCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTAccc                                                                                >  2:4761193/1‑93 (MQ=255)
           ccGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTa                                                                              <  1:7535446/93‑1 (MQ=255)
           ccGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTa                                                                              <  1:3441653/93‑1 (MQ=255)
            cGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTAc                                                                             >  1:7813313/1‑93 (MQ=255)
            cGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTAc                                                                             <  1:6433871/93‑1 (MQ=255)
            cGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACACTAc                                                                             <  2:6963312/93‑1 (MQ=255)
             gAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACa                                                                            <  2:1832025/93‑1 (MQ=255)
              aaaGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACAt                                                                           >  1:1043275/1‑93 (MQ=255)
              aaaGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACAt                                                                           >  1:1922866/1‑93 (MQ=255)
               aaGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATc                                                                          >  2:6118221/1‑93 (MQ=255)
                aGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATcc                                                                         >  1:1442562/1‑93 (MQ=255)
                    ttGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCaa                                                                     >  1:5591391/1‑93 (MQ=255)
                     tGCGTCCCTTCACCGCAAAACTTAAATAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAAt                                                                    <  1:3044086/93‑1 (MQ=255)
                      gCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATa                                                                   <  1:2289261/93‑1 (MQ=255)
                        gTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAAc                                                                 >  2:1943177/1‑93 (MQ=255)
                        gTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAAc                                                                 >  2:2177984/1‑93 (MQ=255)
                         tCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACa                                                                >  1:3185534/1‑93 (MQ=255)
                          cACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACaa                                                               >  1:8461027/1‑93 (MQ=255)
                          cACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACaa                                                               >  1:7220106/1‑93 (MQ=255)
                           aCTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACaaa                                                              >  2:3504391/1‑93 (MQ=255)
                           aCTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACaaa                                                              <  2:722484/93‑1 (MQ=255)
                            cTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACaaaa                                                             >  2:2075392/1‑93 (MQ=255)
                            cTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACaaaa                                                             <  2:4472967/93‑1 (MQ=255)
                             ttCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACacaaa                                                            >  2:360925/1‑93 (MQ=255)
                              tCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAAt                                                           >  1:3307463/1‑93 (MQ=255)
                               cACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATa                                                          >  2:3855217/1‑93 (MQ=255)
                               cACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATa                                                          >  1:2612932/1‑93 (MQ=255)
                               cACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATa                                                          <  1:6995585/93‑1 (MQ=255)
                               cACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATa                                                          <  1:2464697/93‑1 (MQ=255)
                               cACCGCAAAACTTAAAGAAGTTTAATGAATACTACAGAGTGGCTATAATGATGTTAGCCGCTCGCTTACCCTACATCCTCAATAACAAAAATa                                                          <  1:1167637/93‑1 (MQ=255)
                                 ccGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGc                                                        >  1:6558129/1‑93 (MQ=255)
                                 ccGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGc                                                        <  1:4087847/93‑1 (MQ=255)
                                   gCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCt                                                      <  2:4932270/93‑1 (MQ=255)
                                   gCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCt                                                      <  2:4182297/93‑1 (MQ=255)
                                     aaaaCTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTc                                                    >  2:1307611/1‑93 (MQ=255)
                                     aaaaCTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTc                                                    >  1:5348516/1‑93 (MQ=255)
                                      aaaCTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTcc                                                   >  1:3551176/1‑93 (MQ=255)
                                       aaCTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCct                                                  <  1:7648952/93‑1 (MQ=255)
                                       aaCTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCct                                                  >  1:4953333/1‑93 (MQ=255)
                                        aCTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCctc                                                 >  2:1223374/1‑93 (MQ=255)
                                         cTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCctct                                                <  2:6386883/93‑1 (MQ=255)
                                          ttAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTa                                               <  1:1280871/93‑1 (MQ=255)
                                          ttAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTa                                               >  1:4644352/1‑93 (MQ=255)
                                           tAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTaa                                              >  1:8154687/1‑93 (MQ=255)
                                           tAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTaa                                              <  2:2475177/93‑1 (MQ=255)
                                           tAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTaa                                              <  2:5482985/93‑1 (MQ=255)
                                              agaagTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAgg                                           >  2:6229496/1‑93 (MQ=255)
                                              agaagTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAgg                                           <  1:4445455/93‑1 (MQ=255)
                                               gaagTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGt                                          <  1:5490451/93‑1 (MQ=255)
                                               gaagTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGt                                          <  1:8754270/93‑1 (MQ=255)
                                               gaagTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGt                                          <  2:5799704/93‑1 (MQ=255)
                                                aagTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTg                                         <  1:4263995/93‑1 (MQ=255)
                                                 agTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTATAGCTgg                                        >  1:5921519/1‑93 (MQ=255)
                                                  gTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGc                                       >  1:1883204/1‑93 (MQ=255)
                                                    ttaattaaTACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGcat                                     <  2:5070238/93‑1 (MQ=255)
                                                     taattaatACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGcatc                                    <  2:4614012/93‑1 (MQ=255)
                                                     taattaatACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGcatc                                    >  2:4002704/1‑93 (MQ=255)
                                                      aattaatACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGcatca                                   <  1:6778740/93‑1 (MQ=255)
                                                      aattaatACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGcatca                                   <  2:3709023/93‑1 (MQ=255)
                                                      aattaatACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGcatca                                   <  1:2694026/93‑1 (MQ=255)
                                                      aattaatACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGcatca                                   <  1:3907935/93‑1 (MQ=255)
                                                       attaatACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGcatcat                                  >  2:3366177/1‑93 (MQ=255)
                                                       attaatACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGcatcat                                  >  2:4364031/1‑93 (MQ=255)
                                                          aatACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGAc                               <  1:1597838/93‑1 (MQ=255)
                                                            tactacAGTGTTGCTATAAGGATGTTAGCCGCTCTCTTATCCTACATTCTCAATAACAAAAATAGCCTTCCTCTTTAGGTGTCATCTTTACTg                             >  2:8053951/1‑93 (MQ=255)
                                                              ctacAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGAt                           >  1:1746439/1‑93 (MQ=255)
                                                              ctacAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGAt                           >  2:6541734/1‑93 (MQ=255)
                                                              ctacAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGAt                           <  1:1249366/93‑1 (MQ=255)
                                                               tacAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATc                          <  1:2050444/93‑1 (MQ=255)
                                                                acAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCa                         <  2:146903/93‑1 (MQ=255)
                                                                acAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCa                         <  1:5186646/93‑1 (MQ=255)
                                                                   gagTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGc                      >  2:8802345/1‑93 (MQ=255)
                                                                   gagTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGc                      >  1:8631822/1‑93 (MQ=255)
                                                                    agTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCt                     <  2:3376510/93‑1 (MQ=255)
                                                                     gTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTg                    >  2:7808492/1‑93 (MQ=255)
                                                                     gTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTg                    <  1:7885412/93‑1 (MQ=255)
                                                                     gTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTg                    <  1:417925/93‑1 (MQ=255)
                                                                      tGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGa                   <  1:5777213/93‑1 (MQ=255)
                                                                         cTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGaaaa                >  2:7486482/1‑93 (MQ=255)
                                                                          tataAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGaaaaa               >  2:7023988/1‑93 (MQ=255)
                                                                          tataAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTAGCATCATGACTGATCAAGCTGaaaaa               >  1:976326/1‑93 (MQ=255)
                                                                           ataAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGaaaaaa              >  2:2666039/1‑93 (MQ=255)
                                                                           ataAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGaaaaaa              <  2:7734466/93‑1 (MQ=255)
                                                                             aaGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAg            <  2:5926633/93‑1 (MQ=255)
                                                                                gATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCAc         <  1:1083447/93‑1 (MQ=255)
                                                                                 aTGTTAGCCGCCCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCAct        <  1:7420373/93‑1 (MQ=255)
                                                                                  tGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATTAAGCTGAAAAAAAGCActc       <  1:2810529/93‑1 (MQ=255)
                                                                                  tGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCActc       <  1:6112676/93‑1 (MQ=255)
                                                                                  tGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCActc       <  2:8409840/93‑1 (MQ=255)
                                                                                  tGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCActc       <  2:6286904/93‑1 (MQ=255)
                                                                                   gTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCActct      <  1:7975029/93‑1 (MQ=255)
                                                                                          |                                                                                         
TTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCACTCTGCAT  >  NC_000913/3887719‑3887898

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