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.1% 772.2 / ‑4.6 221A206T (GCA→ACA) cadBtransport of lysine/cadaverine
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (108/111);  total (108/113)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 225 total aligned reads displayed.

TGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAG  >  NC_000913/4356904‑4357083
                                                                                          |                                                                                         
tGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTgg                                                                                         >  2:768255/1‑93 (MQ=255)
tGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTgg                                                                                         >  1:2823521/1‑93 (MQ=255)
tGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTgg                                                                                         >  2:173543/1‑93 (MQ=255)
tGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTgg                                                                                         >  2:2611923/1‑93 (MQ=255)
tGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTgg                                                                                         >  1:2107154/1‑93 (MQ=255)
  tAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGAt                                                                                       >  1:2292176/1‑93 (MQ=255)
  tAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGAt                                                                                       >  1:3264357/1‑93 (MQ=255)
  tAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGAt                                                                                       <  2:2168588/93‑1 (MQ=255)
  tAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGAt                                                                                       >  2:4586167/1‑93 (MQ=255)
   aaaCAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGAtt                                                                                      >  1:3564916/1‑93 (MQ=255)
   aaaCAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGAtt                                                                                      >  2:1144087/1‑93 (MQ=255)
    aaCAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTc                                                                                     <  2:1260316/93‑1 (MQ=255)
    aaCAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTc                                                                                     >  1:1241179/1‑93 (MQ=255)
      cAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTcta                                                                                   >  2:4726610/1‑91 (MQ=255)
       aaTACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAac                                                                                  <  1:2951740/93‑1 (MQ=255)
        aTACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAaca                                                                                 >  2:3036462/1‑93 (MQ=255)
         aaCCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAacac                                                                                <  2:3771953/92‑1 (MQ=255)
          aCCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACAcc                                                                               <  2:333629/93‑1 (MQ=255)
           ccTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATGCCAGTACTTACAGCTGTGGATTCAACAccc                                                                              <  2:331541/93‑1 (MQ=255)
           ccTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACAccc                                                                              >  2:2149983/1‑93 (MQ=255)
            cTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCa                                                                             >  2:2896368/1‑93 (MQ=255)
            cTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCa                                                                             <  1:843518/93‑1 (MQ=255)
            cTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCa                                                                             <  1:4031528/93‑1 (MQ=255)
            cTGCTAAACCAGTACCCAGCATGGTGGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATCCCAGTACTTACAGCTGTGGATTCAACACCCa                                                                             <  2:4490935/93‑1 (MQ=255)
             tGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCAc                                                                            >  1:1167700/1‑93 (MQ=255)
             tGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCAc                                                                            <  2:1251373/93‑1 (MQ=255)
             tGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCAc                                                                            >  2:649173/1‑93 (MQ=255)
             tGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCAc                                                                            >  1:369874/1‑93 (MQ=255)
             tGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCAc                                                                            >  1:2889782/1‑93 (MQ=255)
                 aaaCCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAg                                                                        >  1:15644/1‑93 (MQ=255)
                 aaaCCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAg                                                                        <  1:1654703/93‑1 (MQ=255)
                 aaaCCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTAGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACTCACGAAg                                                                        >  2:2221659/1‑93 (MQ=255)
                  aaCCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAgg                                                                       >  1:657981/1‑93 (MQ=255)
                   aCCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGc                                                                      >  1:2530229/1‑93 (MQ=255)
                       gTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCcac                                                                  >  1:2657911/1‑93 (MQ=255)
                        tACCCAGCCTGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCcaca                                                                 <  2:3957793/93‑1 (MQ=255)
                        tACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCcaca                                                                 >  1:718530/1‑93 (MQ=255)
                           ccAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGgc                                                              <  1:782117/93‑1 (MQ=255)
                           ccAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAACTGTGGATTCAACACCCACGAAGGCCCACAGgc                                                              >  2:3573717/1‑93 (MQ=255)
                           ccAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGGTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGgc                                                              >  1:620819/1‑93 (MQ=255)
                             aGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGgcag                                                            >  1:48429/1‑93 (MQ=255)
                              gCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGgcagg                                                           >  2:340213/1‑92 (MQ=255)
                               cATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGgcagag                                                          <  1:2176476/93‑1 (MQ=255)
                               cATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGgcagag                                                          <  1:1893701/93‑1 (MQ=255)
                                 tGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGgcagagca                                                        >  2:3704486/1‑93 (MQ=255)
                                 tGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGgcagagca                                                        >  1:1951412/1‑93 (MQ=255)
                                 tGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGgcagagca                                                        >  1:3905757/1‑93 (MQ=255)
                                   gTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGgcagagcaga                                                      <  2:2060769/93‑1 (MQ=255)
                                    ttGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGaa                                                     <  1:205393/93‑1 (MQ=255)
                                    ttGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGaa                                                     <  1:1824152/93‑1 (MQ=255)
                                    ttGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGaa                                                     <  2:834780/93‑1 (MQ=255)
                                     tGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAAt                                                    >  1:4562392/1‑93 (MQ=255)
                                     tGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAAt                                                    <  2:1705531/93‑1 (MQ=255)
                                      gCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATa                                                   <  1:322335/93‑1 (MQ=255)
                                         aGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACtt                                                <  1:518411/93‑1 (MQ=255)
                                          gCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACttt                                               >  2:3907911/1‑93 (MQ=255)
                                          gCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACttt                                               <  2:3014354/93‑1 (MQ=255)
                                          gCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACttt                                               <  1:379008/93‑1 (MQ=255)
                                            ggAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTa                                             <  1:1324427/93‑1 (MQ=255)
                                             gAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTaa                                            <  1:3044177/93‑1 (MQ=255)
                                             gAACGGGACGTTTCGGGTTTTTAACCATACCAGTCCTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTaa                                            >  2:344967/1‑93 (MQ=255)
                                              aaCGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAAt                                           <  1:1002043/93‑1 (MQ=255)
                                               aCGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATg                                          >  2:2281558/1‑93 (MQ=255)
                                                 ggTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGAt                                        >  2:3535857/1‑93 (MQ=255)
                                                 ggTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGAt                                        >  1:3615559/1‑93 (MQ=255)
                                                  gTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATc                                       >  2:1463122/1‑93 (MQ=255)
                                                   tACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCg                                      <  2:1754910/93‑1 (MQ=255)
                                                    tcGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGc                                     <  1:4594944/92‑1 (MQ=255)
                                                    aCGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGc                                     >  1:451242/1‑93 (MQ=255)
                                                      gTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCAt                                   <  2:2016126/93‑1 (MQ=255)
                                                      gTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCAt                                   >  2:669653/1‑93 (MQ=255)
                                                       tttCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATg                                  >  2:2763090/1‑93 (MQ=255)
                                                         tCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGAc                                >  1:4644756/1‑93 (MQ=255)
                                                           gggTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCa                              >  2:854980/1‑93 (MQ=255)
                                                            ggTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCAt                             <  1:4647622/93‑1 (MQ=255)
                                                            ggTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCAt                             <  1:946461/93‑1 (MQ=255)
                                                             gTTTTTAACCATACCAGTACTTACAGCTGTGGGTTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATc                            >  1:2053556/1‑93 (MQ=255)
                                                             gTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATc                            <  2:3136309/93‑1 (MQ=255)
                                                               ttttAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAg                          <  2:4677832/93‑1 (MQ=255)
                                                                tttAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGt                         <  1:1518787/93‑1 (MQ=255)
                                                                 ttAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTg                        <  1:2633664/93‑1 (MQ=255)
                                                                  tAACAATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGGATACTTTTAATGATCGCATGACCATCAGTgg                       >  2:2873748/1‑93 (MQ=255)
                                                                   aaCCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGt                      <  2:1557363/93‑1 (MQ=255)
                                                                    aCCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTa                     <  2:1083411/93‑1 (MQ=255)
                                                                     ccATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTAt                    <  2:4430544/93‑1 (MQ=255)
                                                                      cATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATc                   >  1:3403252/1‑93 (MQ=255)
                                                                          ccAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCa               <  1:904575/93‑1 (MQ=255)
                                                                           cAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAg              >  1:2641961/1‑93 (MQ=255)
                                                                           cAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAg              >  2:538075/1‑93 (MQ=255)
                                                                            aGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGt             >  2:2412944/1‑93 (MQ=255)
                                                                              tACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTAt           >  2:362930/1‑93 (MQ=255)
                                                                               aCTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTAtt          >  1:3164191/1‑93 (MQ=255)
                                                                                    cAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGt     >  1:2345809/1‑93 (MQ=255)
                                                                                    aaGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGt     <  2:4340039/92‑1 (MQ=255)
                                                                                     aGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGtt    <  2:2105131/93‑1 (MQ=255)
                                                                                     aGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGtt    >  2:898423/1‑93 (MQ=255)
                                                                                      gCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCATTTa   >  1:4075480/1‑93 (MQ=255)
                                                                                      gCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTa   >  2:2491139/1‑93 (MQ=255)
                                                                                       cTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAg  >  2:1741832/1‑93 (MQ=255)
                                                                                       cTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAg  >  2:1178817/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAG  >  NC_000913/4356904‑4357083

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