Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,215,178 Δ1 bp coding (899/930 nt) metA → homoserine O‑transsuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,215,1780A.97.5% 706.8 / 9.8 161coding (899/930 nt)metAhomoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base . (80/77);  total (81/80)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 161 total aligned reads displayed.

TGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTT  >  NC_000913/4215080‑4215272
                                                                                                  |                                                                                              
tGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAt                                                                                               >  2:120999/1‑100 (MQ=255)
   tCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTAc                                                                                           >  1:2754291/1‑100 (MQ=255)
    ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACg                                                                                          >  1:4015392/1‑100 (MQ=255)
     cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACgg                                                                                         >  2:289305/1‑100 (MQ=255)
     cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACgg                                                                                         >  1:3110712/1‑100 (MQ=255)
       cAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGca                                                                                       <  1:3930744/100‑1 (MQ=255)
        aaaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGcac                                                                                      <  2:1892943/100‑1 (MQ=255)
        aaaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGcac                                                                                      <  2:350591/100‑1 (MQ=255)
           aTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATg                                                                                   <  2:665842/100‑1 (MQ=255)
            tACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACACGa                                                                                  <  2:2145523/100‑1 (MQ=255)
             acacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGaa                                                                                 >  1:1693356/1‑100 (MQ=255)
              cacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAAt                                                                                <  1:561013/100‑1 (MQ=255)
               acCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAAt                                                                                <  1:2374619/100‑1 (MQ=255)
               acCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATc                                                                               <  1:2611656/100‑1 (MQ=255)
               acCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATc                                                                               <  2:2255970/100‑1 (MQ=255)
                ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATcc                                                                              <  2:1930922/100‑1 (MQ=255)
                  gcgagcgaGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCaa                                                                            <  2:614496/100‑1 (MQ=255)
                   cgagcgagCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAAc                                                                           <  1:2847519/100‑1 (MQ=255)
                     agcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGc                                                                         <  2:1633164/100‑1 (MQ=255)
                     agcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGc                                                                         <  1:2703009/100‑1 (MQ=255)
                       cgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTg                                                                       >  2:3475422/1‑100 (MQ=255)
                       cgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTg                                                                       >  2:932892/1‑100 (MQ=255)
                        gagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTgg                                                                      >  1:928761/1‑100 (MQ=255)
                        gagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTgg                                                                      >  1:2348573/1‑100 (MQ=255)
                          gcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGAt                                                                    <  2:727334/100‑1 (MQ=255)
                           cTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGAtt                                                                   <  2:1324321/100‑1 (MQ=255)
                           cTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGAtt                                                                   <  2:3721609/100‑1 (MQ=255)
                             ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTaa                                                                 >  1:1765967/1‑100 (MQ=255)
                              gCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAt                                                                <  2:214240/100‑1 (MQ=255)
                               cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtc                                                               <  1:3466135/100‑1 (MQ=255)
                               cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtc                                                               >  2:360815/1‑100 (MQ=255)
                                  aGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtcttc                                                            >  2:73327/1‑100 (MQ=255)
                                   gTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtcttct                                                           >  1:2139728/1‑100 (MQ=255)
                                   gTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtcttct                                                           <  1:4052538/100‑1 (MQ=255)
                                   gTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtcttct                                                           <  2:2359922/100‑1 (MQ=255)
                                     cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                         <  2:3070307/100‑1 (MQ=255)
                                     cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                         >  2:3378365/1‑100 (MQ=255)
                                     cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                         <  1:2424628/100‑1 (MQ=255)
                                     cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                         <  2:1538184/100‑1 (MQ=255)
                                     cACGGTAATTCACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCCCATGAATCCAACGCTGGATTAATCTTCtgt                                                         <  2:2753065/100‑1 (MQ=255)
                                      aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgtg                                                        <  2:622103/100‑1 (MQ=255)
                                       cGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGa                                                       >  1:2422525/1‑100 (MQ=255)
                                            aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTc                                                  >  1:2623880/1‑100 (MQ=255)
                                            aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTc                                                  <  2:1302232/100‑1 (MQ=255)
                                            aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTc                                                  >  2:2234982/1‑100 (MQ=255)
                                              ttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCg                                                 <  2:3926963/99‑1 (MQ=255)
                                                 cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                                             <  1:541429/100‑1 (MQ=255)
                                                    tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGTTTAATCTTCTGTGATAGTCGATCGtt                                           >  1:1101834/1‑99 (MQ=255)
                                                    tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTCATCTTCTGTGATAGTCGATCGTTa                                          >  1:1548629/1‑100 (MQ=255)
                                                    tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTa                                          <  2:1101827/100‑1 (MQ=255)
                                                    tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTa                                          <  2:1794693/100‑1 (MQ=255)
                                                      tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                        >  2:3723096/1‑100 (MQ=255)
                                                       aCCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc                                       >  2:1038749/1‑100 (MQ=255)
                                                        ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCg                                      >  1:1642928/1‑100 (MQ=255)
                                                        ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCg                                      <  1:2233933/100‑1 (MQ=255)
                                                        ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCg                                      >  2:3002282/1‑100 (MQ=255)
                                                        ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCg                                      <  1:2922673/100‑1 (MQ=255)
                                                         cAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGa                                     >  1:1510735/1‑100 (MQ=255)
                                                         cAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGa                                     >  2:2445954/1‑100 (MQ=255)
                                                          aaCTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACCTGAATCAAACGCTGGATTAATCTTCTGTGCTAGTCGATCGTTACGCGAt                                    >  1:1024251/1‑100 (MQ=255)
                                                            atctcTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTc                                  <  1:80692/96‑1 (MQ=255)
                                                             tGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCa                                 <  1:3077068/100‑1 (MQ=255)
                                                              ggCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg                                >  1:3512563/1‑100 (MQ=255)
                                                               gCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGc                               <  2:3477721/100‑1 (MQ=255)
                                                                cTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAgat                              >  2:80692/1‑98 (MQ=255)
                                                                cTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCa                              >  1:89106/1‑100 (MQ=255)
                                                                  cAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAcc                            >  1:3313632/1‑100 (MQ=255)
                                                                  cAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAcc                            >  1:2341636/1‑100 (MQ=255)
                                                                  cAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGAAAGTCGATCGTTAAGCGATTCAGCAcc                            >  2:1996550/1‑100 (MQ=255)
                                                                   aaCTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCt                           >  1:1676801/1‑100 (MQ=255)
                                                                   aaCTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCt                           <  1:1851886/100‑1 (MQ=255)
                                                                     cTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTa                         <  1:998959/100‑1 (MQ=255)
                                                                      tATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTAc                        >  1:1354403/1‑100 (MQ=255)
                                                                       aTTACGTCTATCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTTATTGTCGATCGTTAAGCGATACAGCACCTTAcc                       <  2:2742754/100‑1 (MQ=255)
                                                                         tACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTc                     >  1:2095377/1‑100 (MQ=255)
                                                                         tACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTc                     >  2:2641537/1‑100 (MQ=255)
                                                                          aCGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCa                    <  2:2870621/100‑1 (MQ=255)
                                                                           cGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCa                    <  1:3760993/100‑1 (MQ=255)
                                                                           cGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCa                    >  2:117587/1‑99 (MQ=255)
                                                                             tCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGc                 <  2:294646/100‑1 (MQ=255)
                                                                             tCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGc                 >  1:70244/1‑100 (MQ=255)
                                                                              cTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCa                >  1:507524/1‑100 (MQ=255)
                                                                               tACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAc               >  2:3849295/1‑100 (MQ=255)
                                                                                aCCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAcc              >  1:65946/1‑100 (MQ=255)
                                                                                aCCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAcc              >  2:3590428/1‑100 (MQ=255)
                                                                                aCCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAcc              >  2:1187600/1‑100 (MQ=255)
                                                                                 ccAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCt             >  2:3196086/1‑100 (MQ=255)
                                                                                   aGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTc           >  1:2689049/1‑100 (MQ=255)
                                                                                    gATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCg          >  1:499819/1‑100 (MQ=255)
                                                                                    gATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGAATACACTTCTGTGATAGTCGATCGATAAGAGATTCAGCACCTTACCTCAGGCACCTTCg          >  2:3354518/1‑100 (MQ=255)
                                                                                     aTCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCgg         >  2:2447906/1‑100 (MQ=255)
                                                                                     aTCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCgg         >  2:2379614/1‑100 (MQ=255)
                                                                                      tCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCggg        <  1:3795564/100‑1 (MQ=255)
                                                                                      tCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCggg        >  1:3117647/1‑100 (MQ=255)
                                                                                       cACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGt       <  2:2337799/100‑1 (MQ=255)
                                                                                       cACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCATCTGTGATCGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCggg        >  1:2170357/1‑99 (MQ=255)
                                                                                         cGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGc     >  1:171115/1‑100 (MQ=255)
                                                                                         cGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGc     <  1:2172622/100‑1 (MQ=255)
                                                                                         cGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGc     >  1:2446764/1‑100 (MQ=255)
                                                                                             aTACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCt   <  1:120999/99‑1 (MQ=255)
                                                                                                  |                                                                                              
TGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTT  >  NC_000913/4215080‑4215272

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: