Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,215,147 A→C N290H (AAC→CAC)  metA → homoserine O‑transsuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,215,1470AC96.6% 412.3 / ‑0.8 119N290H (AAC→CAC) metAhomoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base C (61/54);  total (64/55)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 119 total aligned reads displayed.

GGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAG  >  NC_000913/4215050‑4215242
                                                                                                 |                                                                                               
ggATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCAc                                                                                               <  1:1949282/100‑1 (MQ=255)
 gATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACt                                                                                              <  1:1017107/100‑1 (MQ=255)
  aTGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTa                                                                                             <  1:54111/100‑1 (MQ=255)
     tACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTa                                                                                          <  1:328605/100‑1 (MQ=255)
     tACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTa                                                                                          <  1:2113431/100‑1 (MQ=255)
     tACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTa                                                                                          >  1:2000414/1‑100 (MQ=255)
     tACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTa                                                                                          >  2:2903443/1‑100 (MQ=255)
      aCCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTAc                                                                                         <  1:755211/100‑1 (MQ=255)
       ccGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACg                                                                                        <  1:2574718/100‑1 (MQ=255)
       ccGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACg                                                                                        <  1:1925889/100‑1 (MQ=255)
       ccGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACg                                                                                        >  2:2160418/1‑100 (MQ=255)
        cGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGt                                                                                       <  2:2903133/100‑1 (MQ=255)
         gTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTc                                                                                      >  2:427761/1‑100 (MQ=255)
         gTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTc                                                                                      >  1:1928668/1‑100 (MQ=255)
         gTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTc                                                                                      >  1:1935414/1‑100 (MQ=255)
           ataACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTa                                                                                    >  2:1344763/1‑100 (MQ=255)
           ataACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTa                                                                                    >  1:451211/1‑100 (MQ=255)
            taACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTAc                                                                                   >  2:399457/1‑100 (MQ=255)
             aaCTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTAcc                                                                                  <  1:1229724/100‑1 (MQ=255)
              aCTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCa                                                                                 <  1:770814/100‑1 (MQ=255)
                tATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAg                                                                                <  1:1433971/99‑1 (MQ=255)
                 aTTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGAt                                                                              >  1:784466/1‑100 (MQ=255)
                  tttCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATc                                                                             <  1:1330907/100‑1 (MQ=255)
                  tttCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATc                                                                             <  1:2814960/100‑1 (MQ=255)
                     cccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACg                                                                          <  2:2195793/100‑1 (MQ=255)
                     cccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACg                                                                          >  1:1054825/1‑100 (MQ=255)
                     cccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACg                                                                          >  1:2384629/1‑100 (MQ=255)
                      ccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGc                                                                         <  2:1917653/100‑1 (MQ=255)
                          acaATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATa                                                                     <  1:119313/100‑1 (MQ=255)
                          acaATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATa                                                                     <  2:2401198/100‑1 (MQ=255)
                          acaATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATa                                                                     >  2:1216729/1‑100 (MQ=255)
                           caATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATAc                                                                    <  2:2632669/100‑1 (MQ=255)
                            aaTGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACg                                                                   >  2:1739942/1‑100 (MQ=255)
                             aTGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGa                                                                  >  1:2694785/1‑100 (MQ=255)
                             aTGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGa                                                                  >  1:2093069/1‑100 (MQ=255)
                               gATCCGCAAAATACACCGCGAGCGTGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGAt                                                                 <  1:2284101/99‑1 (MQ=255)
                               gATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATc                                                                >  2:225825/1‑100 (MQ=255)
                                aTCCACAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCt                                                               >  1:2904966/1‑100 (MQ=255)
                                  ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTAc                                                             >  2:958846/1‑100 (MQ=255)
                                  ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTAc                                                             <  2:387518/100‑1 (MQ=255)
                                  ccGCAAAATAAACCGCGAGCGAGCTGCCGTAGTCACTGTAATTTACTGTTTACCAACTCGCTCCACTATTACGTCTACCAGATCACGCCATACGATc                                                                <  1:2508153/97‑1 (MQ=255)
                                   cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACg                                                            <  2:1290660/100‑1 (MQ=255)
                                   cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACg                                                            >  2:1333146/1‑100 (MQ=255)
                                     cAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGc                                                          <  1:539967/100‑1 (MQ=255)
                                      aaaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGca                                                         <  1:1385677/100‑1 (MQ=255)
                                          tACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATg                                                     <  2:320061/100‑1 (MQ=255)
                                              ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATc                                                 >  2:1751779/1‑100 (MQ=255)
                                              ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACCATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATc                                                 >  2:2885594/1‑100 (MQ=255)
                                              ccGCGAGCGAACTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATc                                                 <  2:3333/100‑1 (MQ=255)
                                                gcgagcgaGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCa                                               <  1:2725309/100‑1 (MQ=255)
                                                  gagcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAAc                                             >  1:1513229/1‑100 (MQ=255)
                                                  gagcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAAc                                             <  2:20851/100‑1 (MQ=255)
                                                    gcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGc                                           <  2:1740000/100‑1 (MQ=255)
                                                     ngtgcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGACATACGATCTACGGCACATGAATCCAACGCt                                          >  1:1959936/4‑100 (MQ=255)
                                                      gagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTg                                         >  1:1790946/1‑100 (MQ=255)
                                                       agcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTgg                                        <  1:966237/100‑1 (MQ=255)
                                                       agcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTgg                                        >  1:763828/1‑100 (MQ=255)
                                                         cTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGAt                                      >  1:371247/1‑100 (MQ=255)
                                                          tGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCttgttt                                     >  2:914595/1‑95 (MQ=255)
                                                           ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTa                                    <  1:399452/100‑1 (MQ=255)
                                                           ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTa                                    >  1:1216524/1‑100 (MQ=255)
                                                              gTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtc                                 <  2:880636/100‑1 (MQ=255)
                                                                 gTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtcttc                              >  2:1445987/1‑100 (MQ=255)
                                                                 gTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtcttc                              >  2:840550/1‑100 (MQ=255)
                                                                 gTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtcttc                              <  2:1481901/100‑1 (MQ=255)
                                                                 gTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtcttc                              >  1:1566571/1‑100 (MQ=255)
                                                                  tCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtcttct                             >  2:120732/1‑100 (MQ=255)
                                                                   cACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                            <  2:142530/100‑1 (MQ=255)
                                                                   cACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                            >  2:789810/1‑100 (MQ=255)
                                                                    aCGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                           <  2:13915/100‑1 (MQ=255)
                                                                     cGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgtg                          >  1:2200431/1‑100 (MQ=255)
                                                                      ggTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGa                         <  1:1333128/100‑1 (MQ=255)
                                                                       gTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGAt                        >  1:1016238/1‑100 (MQ=255)
                                                                       gTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGAt                        >  1:1036153/1‑100 (MQ=255)
                                                                        tAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATa                       <  1:2449039/100‑1 (MQ=255)
                                                                         aaTTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAg                      <  2:451216/100‑1 (MQ=255)
                                                                          aTTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGt                     <  2:786147/100‑1 (MQ=255)
                                                                          aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATatt                     >  2:960145/1‑98 (MQ=255)
                                                                           tttACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTc                    >  2:1825740/1‑100 (MQ=255)
                                                                            ttACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCg                   <  2:1458631/100‑1 (MQ=255)
                                                                              aCTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGAt                 <  2:784472/100‑1 (MQ=255)
                                                                              aCTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATATTCg                   >  2:202045/1‑98 (MQ=255)
                                                                               cTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATc                <  2:2905004/100‑1 (MQ=255)
                                                                               cTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGCGATAGTCGAt                 >  2:2441597/1‑99 (MQ=255)
                                                                                tGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg               <  1:156605/100‑1 (MQ=255)
                                                                                   ttACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTa            <  1:1758632/100‑1 (MQ=255)
                                                                                   ttACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTa            <  2:1355134/100‑1 (MQ=255)
                                                                                    tACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGtt             >  2:1206697/1‑98 (MQ=255)
                                                                                    tACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTaa           >  1:253829/1‑100 (MQ=255)
                                                                                      ccAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc         >  2:1634214/1‑100 (MQ=255)
                                                                                      ccAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc         >  2:1757937/1‑100 (MQ=255)
                                                                                      ccAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc         >  1:788533/1‑100 (MQ=255)
                                                                                      ccAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc         >  1:2064306/1‑100 (MQ=255)
                                                                                        aaCTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGa       >  2:1280758/1‑100 (MQ=255)
                                                                                        aaCTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGa       <  1:2861218/100‑1 (MQ=255)
                                                                                        aaCTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGa       >  2:2310043/1‑100 (MQ=255)
                                                                                         aCTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGTTAGTCGCTCGTTATGCTAt      >  2:1104636/1‑100 (MQ=255)
                                                                                            ggCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCa   <  2:782689/100‑1 (MQ=255)
                                                                                            ggCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCa   <  2:898179/100‑1 (MQ=255)
                                                                                             gCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg  >  2:2404585/1‑100 (MQ=255)
                                                                                                 |                                                                                               
GGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAG  >  NC_000913/4215050‑4215242

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: