Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,356,994 C→T A206T (GCA→ACA)  cadB ← transport of lysine/cadaverine

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,356,9940CT99.8% 1731.0 / ‑6.2 479A206T (GCA→ACA) cadBtransport of lysine/cadaverine
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (231/246);  total (231/248)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 487 total aligned reads displayed.

TGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAG  >  NC_000913/4356904‑4357083
                                                                                           |                                                                                         
tGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTgg                                                                                         >  2:6906872/1‑93 (MQ=255)
 gTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGa                                                                                        <  1:5985525/93‑1 (MQ=255)
 gTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGa                                                                                        <  2:7131827/93‑1 (MQ=255)
   aaaCAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGAtt                                                                                      <  2:456705/93‑1 (MQ=255)
   aaaCAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGAtt                                                                                      >  1:5492019/1‑93 (MQ=255)
   aaaCAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGAtt                                                                                      >  1:382614/1‑93 (MQ=255)
    aaCAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTc                                                                                     >  1:716953/1‑93 (MQ=255)
     aCAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCa                                                                                    >  1:1484945/1‑93 (MQ=255)
     aCAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCa                                                                                    >  2:9193029/1‑93 (MQ=255)
      cAATCCCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCaa                                                                                   <  1:9891812/93‑1 (MQ=255)
        aTACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAaca                                                                                 >  2:1015532/1‑93 (MQ=255)
        aTACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAaca                                                                                 <  2:5309675/93‑1 (MQ=255)
        aTAACTGCTAAAACAGTAACCAGCATGGTTGCCAGCGGGACGGGACGGTTCGGGGTTTTTAACATAACAGTACTTACAACTGTGGGTTTaaaa                                                                                 >  2:6286792/1‑91 (MQ=255)
           ccTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACAccc                                                                              <  1:6355244/93‑1 (MQ=255)
            cTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCa                                                                             <  2:5608331/93‑1 (MQ=255)
             tGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCAc                                                                            >  1:1885791/1‑93 (MQ=255)
             tGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCAc                                                                            <  1:5058037/93‑1 (MQ=255)
             tGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCAc                                                                            >  1:2505114/1‑93 (MQ=255)
             tGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCAc                                                                            >  1:2589909/1‑93 (MQ=255)
              gCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACg                                                                           >  2:551971/1‑93 (MQ=255)
               cTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGa                                                                          >  1:1202073/1‑93 (MQ=255)
               cTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGa                                                                          >  1:7060994/1‑93 (MQ=255)
                tAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGaa                                                                         <  1:9572485/93‑1 (MQ=255)
                tAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGaa                                                                         <  2:5184691/93‑1 (MQ=255)
                tAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGaa                                                                         <  2:7452637/93‑1 (MQ=255)
                 aaaCCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAg                                                                        >  1:3445294/1‑93 (MQ=255)
                 aaaCCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAg                                                                        <  1:1623947/93‑1 (MQ=255)
                  aaCCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCCCGAAgg                                                                       <  1:3288230/93‑1 (MQ=255)
                  aaCCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAgg                                                                       >  1:9064409/1‑93 (MQ=255)
                   aCCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGc                                                                      >  1:4930592/1‑93 (MQ=255)
                    ccAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGcc                                                                     <  2:9061094/93‑1 (MQ=255)
                    ccAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGcc                                                                     >  2:7898881/1‑93 (MQ=255)
                     cAGTGCCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGccc                                                                    <  1:214306/93‑1 (MQ=255)
                       gTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCcac                                                                  <  1:9231545/93‑1 (MQ=255)
                        tACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCcaca                                                                 >  1:7399372/1‑93 (MQ=255)
                        tACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCcaca                                                                 <  1:3973770/93‑1 (MQ=255)
                        tACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCcaca                                                                 <  2:2154057/93‑1 (MQ=255)
                        tACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCcaca                                                                 <  2:3179332/93‑1 (MQ=255)
                         aCCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAg                                                                <  1:421024/93‑1 (MQ=255)
                          cccAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAgg                                                               <  2:6497638/93‑1 (MQ=255)
                          cccAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAgg                                                               <  2:1757516/93‑1 (MQ=255)
                          cccAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAgg                                                               <  2:503575/93‑1 (MQ=255)
                            cAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGgca                                                             <  2:9608690/93‑1 (MQ=255)
                            cAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGgca                                                             <  1:8937128/93‑1 (MQ=255)
                            cAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGgca                                                             <  2:8499660/93‑1 (MQ=255)
                             aGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGgcag                                                            >  2:3797235/1‑93 (MQ=255)
                               cATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGgcagag                                                          <  1:3409745/93‑1 (MQ=255)
                                aTGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGgcagagc                                                         >  2:7632087/1‑93 (MQ=255)
                                 tGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGgcagagca                                                        >  2:8850920/1‑93 (MQ=255)
                                     tGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAAt                                                    <  1:5167397/93‑1 (MQ=255)
                                      gCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATa                                                   <  1:1134274/93‑1 (MQ=255)
                                           cGGAACGGTATGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCTCAGTCAGAGCAGAATTCtttt                                              <  2:8215038/93‑1 (MQ=255)
                                            ggAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTa                                             >  1:3075011/1‑93 (MQ=255)
                                            ggAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTa                                             <  1:8480597/93‑1 (MQ=255)
                                             gAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTaa                                            <  2:2008821/93‑1 (MQ=255)
                                              aaCGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAAt                                           <  1:10019608/93‑1 (MQ=255)
                                              aaCGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAAt                                           <  2:2269166/93‑1 (MQ=255)
                                              aaCGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAAt                                           <  2:3661901/93‑1 (MQ=255)
                                               aCGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATg                                          >  2:2761586/1‑93 (MQ=255)
                                                cGGTACGTTTCGGGTTTTTAACCATACCAGTTCTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGa                                         <  1:3792326/93‑1 (MQ=255)
                                                cGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAATACCCACAAAGGCACACAGGCACAGCGGAATCTTTTTAATgt                                         >  1:7570432/1‑92 (MQ=255)
                                                cGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGa                                         <  2:1489641/93‑1 (MQ=255)
                                                cGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGa                                         <  1:9436366/93‑1 (MQ=255)
                                                cGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGa                                         <  1:275299/93‑1 (MQ=255)
                                                 ggTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGAt                                        >  2:8641853/1‑93 (MQ=255)
                                                 cgTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGAt                                        <  2:5404364/92‑1 (MQ=255)
                                                    aCGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGc                                     <  2:3964644/93‑1 (MQ=255)
                                                    aCGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGc                                     >  2:4750483/1‑93 (MQ=255)
                                                    aCGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGc                                     <  2:1129430/93‑1 (MQ=255)
                                                     cttttCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCa                                    <  2:2552548/91‑1 (MQ=255)
                                                     cGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCa                                    <  1:6868955/93‑1 (MQ=255)
                                                      gTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCAt                                   >  1:3656154/1‑93 (MQ=255)
                                                       tttCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGTATACTTTTAATGATCGCATg                                  >  1:7692214/1‑93 (MQ=255)
                                                       tttCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATg                                  >  1:6644805/1‑93 (MQ=255)
                                                        ttCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGa                                 >  2:4806637/1‑93 (MQ=255)
                                                         tCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGAc                                >  1:8820578/1‑93 (MQ=255)
                                                           gggTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCa                              >  2:5910889/1‑93 (MQ=255)
                                                           gggTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCa                              <  1:450120/93‑1 (MQ=255)
                                                            ggTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCAt                             >  1:4668646/1‑93 (MQ=255)
                                                             gTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATc                            >  1:6461258/1‑93 (MQ=255)
                                                             gTTTTAACCCATCCCAGACTTAACAGCTGTGAATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATc                            <  2:4079064/93‑1 (MQ=255)
                                                               ttttAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAg                          >  1:819546/1‑93 (MQ=255)
                                                                   aaCCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGt                      >  1:4350710/1‑93 (MQ=255)
                                                                   aaCCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGt                      >  2:3460117/1‑93 (MQ=255)
                                                                      cATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATc                   <  1:9669371/93‑1 (MQ=255)
                                                                       aTACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATcc                  >  2:7557515/1‑93 (MQ=255)
                                                                         aCCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGc                <  2:3910207/93‑1 (MQ=255)
                                                                          ccAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCa               >  2:9434578/1‑93 (MQ=255)
                                                                           cAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAg              >  1:8013463/1‑93 (MQ=255)
                                                                           cAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAg              <  2:9820852/93‑1 (MQ=255)
                                                                             gTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTa            >  1:7788977/1‑93 (MQ=255)
                                                                              tACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTAt           <  2:2467564/93‑1 (MQ=255)
                                                                               aCTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTAtt          <  1:6705552/93‑1 (MQ=255)
                                                                               aCTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTAtt          <  2:1380575/93‑1 (MQ=255)
                                                                                cTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTc         >  2:6985267/1‑93 (MQ=255)
                                                                                 ttACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTcc        >  2:4267522/1‑93 (MQ=255)
                                                                                     cAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGt     >  2:1570506/1‑93 (MQ=255)
                                                                                      aGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGtt    <  1:5658841/93‑1 (MQ=255)
                                                                                       gCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTa   >  2:6696759/1‑93 (MQ=255)
                                                                                       gCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTa   >  2:2075166/1‑93 (MQ=255)
                                                                                           |                                                                                         
TGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAG  >  NC_000913/4356904‑4357083

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