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,8090AG85.0% 721.2 / 111.8 307intergenic (+35/‑56)b3708/tnaBtryptophanase/low affinity tryptophan permease
Reads supporting (aligned to +/- strand):  ref base A (19/27);  new base G (128/132);  total (148/159)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 307 total aligned reads displayed.

TTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCACTCTGCATT  >  NC_000913/3887718‑3887899
                                                                                           |                                                                                          
 ttACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGct                                                                                          >  2:5888863/1‑93 (MQ=255)
    cgtacgAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGctctc                                                                                       <  2:7072841/93‑1 (MQ=255)
     gtacgAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGctctct                                                                                      <  2:521691/93‑1 (MQ=255)
     gtacgAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGctctct                                                                                      <  1:6824668/93‑1 (MQ=255)
     gtacgAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGctctct                                                                                      <  1:7033394/93‑1 (MQ=255)
     gtacgAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCActctct                                                                                      <  1:2541058/93‑1 (MQ=255)
       acgAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTa                                                                                    <  2:1360241/93‑1 (MQ=255)
       acgAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTa                                                                                    <  1:2725155/93‑1 (MQ=255)
         gAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTAcc                                                                                  <  1:115371/93‑1 (MQ=255)
           aCCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCt                                                                                >  2:678360/1‑93 (MQ=255)
            ccGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTa                                                                               <  1:4299976/93‑1 (MQ=255)
              gAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACa                                                                             >  1:4174271/1‑93 (MQ=255)
              gAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACa                                                                             >  1:5551294/1‑93 (MQ=255)
              gAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACa                                                                             >  2:4200244/1‑93 (MQ=255)
               aaaGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACAt                                                                            >  1:1944915/1‑93 (MQ=255)
                 aGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATcc                                                                          <  1:5767573/93‑1 (MQ=255)
                 aGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATcc                                                                          <  1:879831/93‑1 (MQ=255)
                  gTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCt                                                                         <  2:5051254/93‑1 (MQ=255)
                   tATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTc                                                                        >  1:1724534/1‑93 (MQ=255)
                    aTTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCa                                                                       >  1:5969963/1‑93 (MQ=255)
                     ttGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCaa                                                                      >  1:2090300/1‑93 (MQ=255)
                     ttGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCaa                                                                      >  1:727047/1‑93 (MQ=255)
                      tGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAAt                                                                     <  1:6525493/93‑1 (MQ=255)
                      tGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAAt                                                                     <  1:805126/93‑1 (MQ=255)
                       gCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATa                                                                    >  2:1647320/1‑93 (MQ=255)
                       gCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATa                                                                    <  2:5403884/93‑1 (MQ=255)
                         gTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAAc                                                                  <  2:4428732/93‑1 (MQ=255)
                          tCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACa                                                                 <  1:4591694/93‑1 (MQ=255)
                          tCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACa                                                                 >  2:3914187/1‑93 (MQ=255)
                            aCTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACaaa                                                               <  2:1876376/93‑1 (MQ=255)
                             cTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACaaaa                                                              >  2:7382566/1‑93 (MQ=255)
                              ttCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACaaaaa                                                             >  2:317/1‑93 (MQ=255)
                              ttCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACaaaaa                                                             >  1:1670169/1‑93 (MQ=255)
                              ttCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACaaaaa                                                             <  2:1800222/93‑1 (MQ=255)
                              ttCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACaaaaa                                                             <  2:6041489/93‑1 (MQ=255)
                                cACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATa                                                           <  1:2312780/93‑1 (MQ=255)
                                 aCCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAg                                                          >  1:3543579/1‑93 (MQ=255)
                                  ccGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGc                                                         <  1:2937063/93‑1 (MQ=255)
                                  ccGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGc                                                         >  1:3213422/1‑93 (MQ=255)
                                  ccGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGc                                                         <  2:3695631/93‑1 (MQ=255)
                                   cGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGcc                                                        <  1:4832189/93‑1 (MQ=255)
                                   cGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGcc                                                        >  1:1517767/1‑93 (MQ=255)
                                    gCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCt                                                       >  2:4795450/1‑93 (MQ=255)
                                    gCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCt                                                       <  1:7415616/93‑1 (MQ=255)
                                    gCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCt                                                       <  2:5511890/93‑1 (MQ=255)
                                     cAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCtt                                                      <  1:63081/93‑1 (MQ=255)
                                       aaaCTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTcc                                                    >  1:3677883/1‑93 (MQ=255)
                                       aaaCTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTcc                                                    <  1:4429580/93‑1 (MQ=255)
                                         aCTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCctc                                                  >  1:5807403/1‑93 (MQ=255)
                                         aCTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCctc                                                  <  1:4700394/93‑1 (MQ=255)
                                          cTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCctct                                                 <  2:704437/93‑1 (MQ=255)
                                            tAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTaa                                               >  2:4348105/1‑93 (MQ=255)
                                            tAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTaa                                               >  2:2681794/1‑93 (MQ=255)
                                            tAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTaa                                               <  2:6255548/93‑1 (MQ=255)
                                             aaaGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTaaa                                              <  1:5889853/93‑1 (MQ=255)
                                             aaaGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTaaa                                              <  2:5288723/93‑1 (MQ=255)
                                              aagaagTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAg                                             >  2:4758788/1‑93 (MQ=255)
                                              aagaagTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAg                                             >  1:7193157/1‑93 (MQ=255)
                                                gaagTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGt                                           <  1:2158926/93‑1 (MQ=255)
                                                gaagTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGt                                           <  1:6965890/93‑1 (MQ=255)
                                                gaagTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGt                                           <  2:5611904/93‑1 (MQ=255)
                                                gaagTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGt                                           >  2:2389570/1‑93 (MQ=255)
                                                  agTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTgg                                         <  2:6121026/93‑1 (MQ=255)
                                                    tttAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGca                                       >  1:3499653/1‑93 (MQ=255)
                                                      taattaatACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGcatc                                     >  2:6714561/1‑93 (MQ=255)
                                                      taattaatACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGcatc                                     <  1:3959115/93‑1 (MQ=255)
                                                       aattaatACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGcatca                                    <  1:1175141/93‑1 (MQ=255)
                                                       aattaatACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGcatca                                    >  2:3475693/1‑93 (MQ=255)
                                                       aattaatACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGcatca                                    <  1:663804/93‑1 (MQ=255)
                                                       aattaatACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGcatca                                    <  1:4904321/93‑1 (MQ=255)
                                                          taatACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGa                                 <  2:838289/93‑1 (MQ=255)
                                                             tactacAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTg                              >  2:2779395/1‑93 (MQ=255)
                                                             tactacAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTg                              >  2:2103365/1‑93 (MQ=255)
                                                             tactacAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTg                              >  1:22822/1‑93 (MQ=255)
                                                               ctacAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGAt                            <  1:5392889/93‑1 (MQ=255)
                                                                tacAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATc                           >  2:3518535/1‑93 (MQ=255)
                                                                tacAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATc                           >  1:5021683/1‑93 (MQ=255)
                                                                  cAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCaa                         >  2:6587258/1‑93 (MQ=255)
                                                                  cAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCaa                         <  1:6570274/93‑1 (MQ=255)
                                                                  cAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCaa                         <  2:2252986/93‑1 (MQ=255)
                                                                   agagTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAg                        <  1:182512/93‑1 (MQ=255)
                                                                   agagTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAg                        >  1:4134675/1‑93 (MQ=255)
                                                                    gagTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGc                       <  1:962731/93‑1 (MQ=255)
                                                                     agTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCt                      <  2:5122222/93‑1 (MQ=255)
                                                                       tGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGa                    <  1:138784/93‑1 (MQ=255)
                                                                       tGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGa                    <  1:5664202/93‑1 (MQ=255)
                                                                           tataAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGaaaaa                <  1:5192078/93‑1 (MQ=255)
                                                                           tataAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGaaaaa                >  2:3117502/1‑93 (MQ=255)
                                                                                ggATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCa           >  2:4622052/1‑93 (MQ=255)
                                                                                ggATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCa           <  2:1097057/93‑1 (MQ=255)
                                                                                ggATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCa           <  1:3625010/93‑1 (MQ=255)
                                                                                  aTGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAATAAGCAct         >  2:1943074/1‑93 (MQ=255)
                                                                                   tGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCActc        <  1:5283243/93‑1 (MQ=255)
                                                                                   tGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCActc        >  2:1423029/1‑93 (MQ=255)
                                                                                    gTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCActct       >  1:4069841/1‑93 (MQ=255)
                                                                                    gTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCActct       <  1:2984024/93‑1 (MQ=255)
                                                                                      tAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCACTCTGc     <  1:3034635/93‑1 (MQ=255)
                                                                                      tAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCACTCTGc     <  2:2959925/93‑1 (MQ=255)
                                                                                       aGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCACTCTGCa    <  2:2080679/93‑1 (MQ=255)
                                                                                        gCCGCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCACTCTGCAt   >  1:7340457/1‑93 (MQ=255)
                                                                                           |                                                                                          
TTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCACTCTGCATT  >  NC_000913/3887718‑3887899

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