Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,100,447 (T)6→5 coding (79/1035 nt) rhaT ← L‑rhamnose:proton symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,100,4420T.93.8% 369.8 / 17.9 96coding (84/1035 nt)rhaTL‑rhamnose:proton symporter
Reads supporting (aligned to +/- strand):  ref base T (4/2);  new base . (40/50);  total (44/52)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

AAATTCGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTCTTATCCT  >  NC_000913/4100344‑4100536
                                                                                                  |                                                                                              
aaaTTCGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACtt                                                                                               <  2:1447407/100‑1 (MQ=255)
 aaTTCGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACttt                                                                                              >  1:1043339/1‑100 (MQ=255)
 aaTTCGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACttt                                                                                              <  2:1948473/100‑1 (MQ=255)
  aTTCGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACtttt                                                                                             >  1:1978807/1‑100 (MQ=255)
   ttCGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACttttt                                                                                            <  1:2137107/100‑1 (MQ=255)
   ttCGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACttttt                                                                                            <  1:1775120/100‑1 (MQ=255)
   ttCGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACttttt                                                                                            <  2:266918/100‑1 (MQ=255)
    tCGGTAGTAACAGGGCGCTGCTGGCCCACGGCAGAATATTCCACGAAACAATCCCAGCGACTGACCACATGGTTTCCCATGATCATTTTTTTACtttttg                                                                                           <  2:1725790/100‑2 (MQ=255)
     cGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACtttttga                                                                                          <  2:734151/100‑3 (MQ=255)
     cGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACtttttga                                                                                          <  1:1431815/100‑3 (MQ=255)
      ggTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTAtttttt                                                                                            >  2:1066789/1‑97 (MQ=255)
       gtagtaACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGaa                                                                                        >  1:613026/1‑100 (MQ=255)
        tagtaACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAACCAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACg                                                                                      <  2:1850608/100‑1 (MQ=255)
        tagtaACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACg                                                                                      <  2:1499578/100‑1 (MQ=255)
        tagtaACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACg                                                                                      <  1:2751508/100‑1 (MQ=255)
         agtaACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACgg                                                                                     <  1:1476760/100‑1 (MQ=255)
         agtaACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACgg                                                                                     >  2:2189587/1‑100 (MQ=255)
          gtaACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGCCCATTTTTTTACTTTTTGAACGGa                                                                                    >  1:1091391/1‑100 (MQ=255)
          gtaACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGa                                                                                    >  1:2051083/1‑100 (MQ=255)
           taACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAg                                                                                   >  2:1933823/1‑100 (MQ=255)
           taACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAg                                                                                   <  2:2818351/100‑1 (MQ=255)
           naaCAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTATCCCATGACCATTTTTTTACTTTTTGAACGGAg                                                                                   >  2:2270236/2‑100 (MQ=255)
            aaCAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATttttttttttttttttcggag                                                                                    >  2:1143943/1‑92 (MQ=255)
             aCAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGACCGGAGc                                                                                  >  1:800218/1‑99 (MQ=255)
             aCAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCg                                                                                 <  2:1180269/100‑1 (MQ=255)
              cAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGt                                                                                <  1:897164/100‑1 (MQ=255)
                   cgcTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTaaaaa                                                                           >  1:466821/1‑100 (MQ=255)
                      tGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAg                                                                        <  2:182119/100‑1 (MQ=255)
                      tGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAg                                                                        >  1:2235592/1‑100 (MQ=255)
                        aTGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGc                                                                      <  2:535884/100‑1 (MQ=255)
                          ggCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTg                                                                    >  1:1090434/1‑100 (MQ=255)
                           gCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATTACCATTTTTTGACTTTTTGAACGGAGCGTAAAAACAGGCttc                                                                   <  1:1950402/100‑3 (MQ=255)
                           gCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGc                                                                   <  2:1978793/100‑1 (MQ=255)
                           gCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGc                                                                   >  2:1389990/1‑100 (MQ=255)
                            cccACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTAACGGAGCGTAAAAACAGGCTGCa                                                                  >  2:336289/1‑100 (MQ=255)
                             ccACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCan                                                                 <  2:590729/100‑2 (MQ=255)
                              cACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCAc                                                                 <  1:1268939/100‑1 (MQ=255)
                              cACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTAACGGAGCGTAAAAACAGGCTGCACt                                                                >  2:1727185/1‑100 (MQ=255)
                               aCGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACt                                                                <  1:242636/99‑1 (MQ=255)
                               aCGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTg                                                               <  1:1782721/100‑1 (MQ=255)
                                 ggCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGc                                                             <  1:2674512/100‑1 (MQ=255)
                                  gCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGcc                                                            <  1:2532721/100‑1 (MQ=255)
                                  gCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGcc                                                            >  1:85558/1‑100 (MQ=255)
                                  gCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGcc                                                            >  2:1733586/1‑100 (MQ=255)
                                  gCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGcc                                                            <  1:2629753/100‑1 (MQ=255)
                                  gCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTACCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGcc                                                            <  2:196842/100‑1 (MQ=255)
                                     gAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCcgcg                                                         >  2:2785627/1‑100 (MQ=255)
                                      aataatCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCgcgc                                                        <  1:2688009/100‑1 (MQ=255)
                                       ataatCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGcc                                                       <  2:1053507/100‑1 (MQ=255)
                                        taatCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCg                                                      <  2:356964/100‑1 (MQ=255)
                                          atCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGAt                                                    <  1:490302/100‑1 (MQ=255)
                                           tCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATc                                                   <  2:1024275/100‑1 (MQ=255)
                                            ccACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATc                                                   <  2:2582848/100‑1 (MQ=255)
                                            ccACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCa                                                  >  2:2749295/1‑100 (MQ=255)
                                            ccACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCa                                                  <  2:1583019/100‑1 (MQ=255)
                                             cACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCaa                                                 >  1:2696182/1‑100 (MQ=255)
                                               cGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAAt                                               <  2:2201705/100‑1 (MQ=255)
                                                    cAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCaa                                          <  2:418655/100‑1 (MQ=255)
                                                      aTCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCaaaa                                        >  2:540498/1‑100 (MQ=255)
                                                      aTCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCaaaa                                        <  1:127152/100‑1 (MQ=255)
                                                          cACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATc                                    >  1:1417461/1‑100 (MQ=255)
                                                           aCCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATcc                                   >  1:1909738/1‑100 (MQ=255)
                                                              gACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCa                                >  1:416009/1‑100 (MQ=255)
                                                              gACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCa                                >  1:416008/1‑100 (MQ=255)
                                                              gACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCa                                >  2:1661961/1‑100 (MQ=255)
                                                                 tGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCATAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCg                             <  2:144251/100‑1 (MQ=255)
                                                                  gACCACATGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGt                            >  1:1083457/1‑100 (MQ=255)
                                                                     cacaTGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAAt                         >  2:1610662/1‑100 (MQ=255)
                                                                     cacaTGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAAt                         >  1:963768/1‑100 (MQ=255)
                                                                     cacaTGGTTTCCCATGACCATTTTTTTTTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAAt                         >  2:2842373/1‑100 (MQ=255)
                                                                      acaTGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATc                        >  2:2846386/1‑100 (MQ=255)
                                                                      acaTGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATc                        >  1:305462/1‑100 (MQ=255)
                                                                       caTGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATcg                       >  1:863670/1‑100 (MQ=255)
                                                                        aTGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATcgc                      <  2:2696162/100‑1 (MQ=255)
                                                                         tGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTCATcgcg                     <  2:2452349/100‑1 (MQ=255)
                                                                         tGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATcgcg                     <  2:601770/100‑1 (MQ=255)
                                                                         tGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATcgcg                     >  1:125155/1‑100 (MQ=255)
                                                                         tGGTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATcgcg                     <  1:1843600/100‑1 (MQ=255)
                                                                           gTTTCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGtt                   >  1:1329642/1‑100 (MQ=255)
                                                                            tttCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTa                  <  1:1727198/100‑1 (MQ=255)
                                                                             ttCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTAc                 <  2:2824706/100‑1 (MQ=255)
                                                                             ttCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTAc                 <  1:1126005/100‑1 (MQ=255)
                                                                             ttCCCATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTAc                 <  1:780042/100‑1 (MQ=255)
                                                                                ccATGACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCa              >  2:2272371/1‑100 (MQ=255)
                                                                                    gACCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAACGCCAAAATATCCCCATCGTAATCGCGTTACTCAtctt          >  2:1947358/1‑100 (MQ=255)
                                                                                     aCCATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCAtcttc         <  2:816593/100‑1 (MQ=255)
                                                                                      ccATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCAtcttct        <  2:2349257/100‑1 (MQ=255)
                                                                                      ccATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCAtcttc         <  2:800217/99‑1 (MQ=255)
                                                                                       cATTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATcttctt       >  1:976045/1‑100 (MQ=255)
                                                                                        aTTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTCTTa      <  1:2846407/100‑1 (MQ=255)
                                                                                        aTTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTCTTa      <  2:416007/100‑1 (MQ=255)
                                                                                        aTTTTTTTACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTCTTa      <  2:416008/100‑1 (MQ=255)
                                                                                           tttttACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTCTTATcc   <  2:322449/100‑1 (MQ=255)
                                                                                           tttttACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTCTTATcc   <  2:2430704/100‑1 (MQ=255)
                                                                                           tttttACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTCTTATcc   <  2:466819/100‑1 (MQ=255)
                                                                                           tttttACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTCTTATcc   <  1:1389995/100‑1 (MQ=255)
                                                                                           tttttACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTCTTATcc   <  2:1043337/100‑1 (MQ=255)
                                                                                           tttttACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTCTTATcc   <  2:1563157/100‑1 (MQ=255)
                                                                                            ttttACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTCTTATCCt  >  1:527861/1‑100 (MQ=255)
                                                                                            ttttACTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTCTTATCCt  >  1:527160/1‑100 (MQ=255)
                                                                                                  |                                                                                              
AAATTCGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTCTTATCCT  >  NC_000913/4100344‑4100536

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