Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,215,178 Δ1 bp 29.9% 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.29.9% 239.1 / 109.0 97coding (899/930 nt)metAhomoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand):  ref base A (41/26);  new base . (18/11);  total (59/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

GATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTA  >  NC_000913/4215081‑4215277
                                                                                                 |                                                                                                   
gATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTCCCAGATCACGCCATACGATc                                                                                                   >  1:1276138/1‑100 (MQ=255)
   ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTAc                                                                                                >  1:1405845/1‑100 (MQ=255)
    cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACg                                                                                               <  2:487830/100‑1 (MQ=255)
     gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACgg                                                                                              >  1:260683/1‑100 (MQ=255)
     gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGc                                                                                             >  2:2030296/1‑100 (MQ=255)
      cAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGc                                                                                             <  1:368477/100‑1 (MQ=255)
       aaaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGcac                                                                                           <  2:1384961/100‑1 (MQ=255)
        aaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGcac                                                                                           <  2:2319383/100‑1 (MQ=255)
        aaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGcaca                                                                                          >  1:103/1‑100 (MQ=255)
         aaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACAt                                                                                         <  2:822245/100‑1 (MQ=255)
          aTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATg                                                                                        >  1:2134173/1‑100 (MQ=255)
           tACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTAGTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGa                                                                                       >  2:1252576/1‑100 (MQ=255)
           tACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATg                                                                                        <  1:534333/100‑1 (MQ=255)
           tACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGa                                                                                       >  1:1307178/1‑100 (MQ=255)
            acacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGa                                                                                       >  2:1185317/1‑100 (MQ=255)
             cacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACAt                                                                                         >  1:266479/1‑97 (MQ=255)
              acCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAAt                                                                                     <  1:1276989/100‑1 (MQ=255)
               ccGCGATCGAGCTGGCGTAGTCACGGTAATTTACTGCTTACCAACTGGCTCAATTATTACGTCTACCATATCACGCCATACGATCTACGGCACATGAATc                                                                                    <  2:531235/100‑1 (MQ=255)
               ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATc                                                                                    >  2:346869/1‑100 (MQ=255)
               ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATc                                                                                    >  1:1862830/1‑100 (MQ=255)
                cgcgAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATcc                                                                                   <  2:579803/100‑1 (MQ=255)
                       gagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTgg                                                                           >  1:769018/1‑100 (MQ=255)
                       gagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACGGGCGCAACTATTACGTCTACCAGATCACACCATCCCATCTACGTCACATGAATCCa                                                                                  >  1:2043146/1‑94 (MQ=255)
                           tGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGAtt                                                                        <  2:867970/100‑1 (MQ=255)
                            ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTa                                                                       >  1:1004664/1‑100 (MQ=255)
                            ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTa                                                                       >  1:2091484/1‑100 (MQ=255)
                              cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAt                                                                     >  1:2332315/1‑100 (MQ=255)
                              cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtc                                                                    >  2:1678713/1‑100 (MQ=255)
                               gTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtc                                                                    <  2:341701/100‑1 (MQ=255)
                                 aGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtctt                                                                  <  1:2349744/100‑1 (MQ=255)
                                  gTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtcttc                                                                 >  2:858478/1‑100 (MQ=255)
                                  gTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtcttc                                                                 <  2:1030734/100‑1 (MQ=255)
                                   tCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACAATCTACGGAACATGAATCCAACGCGGGATTAAttttc                                                                 >  1:2340583/1‑99 (MQ=255)
                                   tCACGGTAATTTACTGTTTACCAACTGGCAAAACTATTACGTCTACCAGACCACGCCATACGTCTACGGCACATGAATCCAACACTGGATTAAtcttc                                                                 <  2:1430239/98‑1 (MQ=255)
                                    naCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                                                               >  2:150723/2‑100 (MQ=255)
                                    cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGCATCCAACGCTGGATTAAtcttct                                                                >  1:114663/1‑99 (MQ=255)
                                    cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtcttct                                                                >  1:1390278/1‑99 (MQ=255)
                                    cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgn                                                              <  2:103/100‑2 (MQ=255)
                                    cACGGTAATTTACTCTTTACCAACTGGCTCAACTATTCCGGCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                                                               <  2:710120/100‑1 (MQ=255)
                                      agttattttACTTTTTCCCAATTTGGTCAACTATTAGTTCTACGATACAACGCCATACGTCTATGTCACATGAATCCAACGCTGGATTAATGTTAtg                                                               <  2:483364/91‑1 (MQ=255)
                                        gTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGAt                                                           <  2:1336568/100‑1 (MQ=255)
                                        gTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGAt                                                           <  1:1461789/100‑1 (MQ=255)
                                        gTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGAt                                                           >  1:2036258/1‑100 (MQ=255)
                                        gTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGAt                                                           <  1:250851/100‑1 (MQ=255)
                                         tAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGAt                                                           >  1:766654/1‑99 (MQ=255)
                                         tAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAg                                                         <  1:904168/100‑1 (MQ=255)
                                          aaTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAg                                                         >  2:815193/1‑100 (MQ=255)
                                           aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGAAAATGAATCCAACGCTGGATTAATCTTCTGTGATAGt                                                        >  1:1697778/1‑100 (MQ=255)
                                            tttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGTTatt                                                        >  1:631959/1‑97 (MQ=255)
                                             ttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCg                                                      >  1:2199694/1‑100 (MQ=255)
                                             ttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCg                                                      <  2:2036243/100‑1 (MQ=255)
                                             ttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGa                                                     >  1:1390922/1‑100 (MQ=255)
                                                cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                                                  >  2:1642358/1‑100 (MQ=255)
                                                cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                                                  >  2:2263691/1‑100 (MQ=255)
                                                cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                                                  >  1:20495/1‑100 (MQ=255)
                                                  gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGCTAGTcttttgt                                                 >  1:1516391/1‑94 (MQ=255)
                                                  gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATc                                                   >  2:562019/1‑98 (MQ=255)
                                                  gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGt                                                 >  1:1890600/1‑100 (MQ=255)
                                                  gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGt                                                 <  2:487570/100‑1 (MQ=255)
                                                   tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGtt                                                >  2:2207622/1‑100 (MQ=255)
                                                    ttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTaa                                              <  1:891654/100‑1 (MQ=255)
                                                     tACCAACTGGCTCAACTATTACGTCTACTAGATTACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTTTGATAGTTGATCGTTaa                                              >  2:1291407/1‑100 (MQ=255)
                                                     tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATCCGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTa                                               <  1:1404071/99‑1 (MQ=255)
                                                     tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTaa                                              <  2:1180478/100‑1 (MQ=255)
                                                     naCCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTaa                                              >  2:2178001/2‑100 (MQ=255)
                                                       ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTaatc                                            >  1:1360996/1‑98 (MQ=255)
                                                       ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc                                            >  2:889276/1‑100 (MQ=255)
                                                       ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc                                            <  1:2055492/100‑1 (MQ=255)
                                                        cAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCg                                           <  1:2000561/100‑1 (MQ=255)
                                                          aCTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGAt                                         <  1:2126469/100‑1 (MQ=255)
                                                           cTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTc                                       <  2:1737839/100‑1 (MQ=255)
                                                            tGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCa                                      >  1:933569/1‑100 (MQ=255)
                                                            tGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCa                                      >  1:2010362/1‑100 (MQ=255)
                                                            tGGCTCAACTATTACGTCTACCAGATCACACCATACGATCTACGGCACATGAATCCAACGCTGGATGAATCTTCTGTGATAGTCGATCGTTAAGCGATTc                                       <  2:1415156/100‑1 (MQ=255)
                                                             ggCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTACTCTTCTGTGATAGTCGATCGTTAAGCGATTCa                                      >  2:339538/1‑100 (MQ=255)
                                                             ggCTCAACTATTAAGTCTACCAGATCACGCCATACGATATACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCa                                      >  2:1453106/1‑100 (MQ=255)
                                                                 cAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCCGGATTAATCTTCGGTGCTAGTCGATCGTTAAGCGCTTCAg                                     >  2:815331/1‑97 (MQ=255)
                                                                 cAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTGAGCAcc                                 >  2:1360946/1‑100 (MQ=255)
                                                                  aaCTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAcc                                 <  1:740292/100‑1 (MQ=255)
                                                                    cTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTa                              >  1:2084433/1‑100 (MQ=255)
                                                                     tATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTa                              >  1:1441095/1‑100 (MQ=255)
                                                                      aTTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCtt                               >  1:1522441/1‑98 (MQ=255)
                                                                         aCGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATAGTTAAGCGATTCAGCACCTTACCTCa                         >  1:2132013/1‑100 (MQ=255)
                                                                            tCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGc                      >  1:318205/1‑100 (MQ=255)
                                                                            tCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGc                      <  1:1740437/100‑1 (MQ=255)
                                                                              tACCAGATCACGCCATACGCTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCa                     <  1:2028806/100‑1 (MQ=255)
                                                                               aCCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAcc                   >  1:909824/1‑100 (MQ=255)
                                                                               aCCACATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGAGCGTTAAGCGATTCAGCACCTTCCCTCATGc                      <  2:420004/97‑1 (MQ=255)
                                                                                ccAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAcc                   >  1:119567/1‑100 (MQ=255)
                                                                                   gATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTc                >  1:2096683/1‑100 (MQ=255)
                                                                                    aTCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCg               <  1:448553/100‑1 (MQ=255)
                                                                                     tCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCgg              >  1:27405/1‑100 (MQ=255)
                                                                                     tCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGAGAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCggg             <  2:1307172/100‑1 (MQ=255)
                                                                                        cGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTg           >  2:1349323/1‑100 (MQ=255)
                                                                                           cATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGGTTCAGCACCTTACCTCAGGCACCTTAgg              >  2:1822449/1‑94 (MQ=255)
                                                                                           cATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGcc         <  2:1862817/99‑1 (MQ=255)
                                                                                           cATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCt        >  2:717939/1‑100 (MQ=255)
                                                                                                 aTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTa  <  1:645107/100‑1 (MQ=255)
                                                                                                 aTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTa  >  2:269490/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
GATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTA  >  NC_000913/4215081‑4215277

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