Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,215,178 Δ1 bp 100% 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.0% 707.0 / 9.0 162coding (899/930 nt)metAhomoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base . (83/74);  total (84/78)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

Only 100 of 162 total aligned reads displayed.

TCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCT  >  NC_000913/4215083‑4215271
                                                                                               |                                                                                             
tCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTAc                                                                                          >  2:204673/1‑100 (MQ=255)
 ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACg                                                                                         <  2:3902042/100‑1 (MQ=255)
  cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACgg                                                                                        <  1:325604/100‑1 (MQ=255)
   gcgaaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGc                                                                                       <  1:1030638/97‑1 (MQ=255)
   gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGc                                                                                       >  1:3057876/1‑100 (MQ=255)
    cAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGca                                                                                      <  1:11016/100‑1 (MQ=255)
      aaaTACACCGCTAGCGAGCTGGCGTAGGCACTGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGcaca                                                                                    <  1:3273320/100‑1 (MQ=255)
      aaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGcaca                                                                                    >  1:3096624/1‑100 (MQ=255)
       aaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGcaca                                                                                    <  1:1427293/100‑1 (MQ=255)
        aTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATg                                                                                  >  1:3374347/1‑100 (MQ=255)
        aTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATg                                                                                  <  1:353024/100‑1 (MQ=255)
          acacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGaa                                                                                >  1:3417582/1‑100 (MQ=255)
          acacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTAACAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACAGCACATGaa                                                                                >  2:2579006/1‑100 (MQ=255)
           cacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACAGCACATTAAt                                                                               >  2:2916580/1‑100 (MQ=255)
             ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATcc                                                                             >  1:1959853/1‑100 (MQ=255)
             ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATcc                                                                             >  1:2033613/1‑100 (MQ=255)
                     gagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTgg                                                                     >  1:451732/1‑100 (MQ=255)
                      agcTGGCGTAGTCACGGTAATTTACTGTTTACAAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTCCGGCACATGAATCCAAAGCTGGa                                                                    >  2:239724/1‑100 (MQ=255)
                       gcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGAt                                                                   <  2:2011719/100‑1 (MQ=255)
                        cTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTAAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGAtt                                                                  >  1:3860159/1‑100 (MQ=255)
                          ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTaa                                                                >  1:2104165/1‑100 (MQ=255)
                          ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACACCATACG‑TCTACGGCACATGAATCCAACGCTGGATTaa                                                                >  1:36749/1‑100 (MQ=255)
                               aGTCACGGTAATTTACTGTTTACCGACTGGCTTAACTATGAAGTCTACCAGATCACTCCATACG‑TCTACGGCACATGATTCCAACGCTGTATTAAtcttc                                                           <  2:3862041/100‑1 (MQ=255)
                               aGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtcttc                                                           <  2:2178832/100‑1 (MQ=255)
                                 tCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                                                         <  2:1805670/100‑1 (MQ=255)
                                 tCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                                                         <  1:3563366/100‑1 (MQ=255)
                                 tCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                                                         <  1:1619368/100‑1 (MQ=255)
                                 tCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                                                         <  1:2955510/100‑1 (MQ=255)
                                 tCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                                                         <  1:2059806/100‑1 (MQ=255)
                                  cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                        >  2:1505751/1‑100 (MQ=255)
                                   aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgtg                                                       <  1:114816/100‑1 (MQ=255)
                                   aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                        <  2:526269/99‑1 (MQ=255)
                                      gTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATa                                                    <  2:763650/100‑1 (MQ=255)
                                       tAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAg                                                   <  2:282262/100‑1 (MQ=255)
                                        aaTTTACTGTTTTCCAACTGGCTCAACTATTACTTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCGGTGAt                                                     <  2:483783/97‑1 (MQ=255)
                                        aaTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGt                                                  >  1:2479565/1‑100 (MQ=255)
                                         aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTc                                                 <  2:1552280/100‑1 (MQ=255)
                                         aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATAAATACAACGCTGGATTAATATTCTGTGATAGTc                                                 >  1:1904298/1‑100 (MQ=255)
                                         aTTTACTGTTTACCAACTGGATCAACTATTACGTCTACCAGAAAAAGCCATACG‑TCTACGGCACATGAATCCAAAGCTGGATTAATCTTCTGTGATAGTc                                                 >  2:1148586/1‑100 (MQ=255)
                                          tttACTGTTTACCAACTGGCTCAACTATTACGACTACCAGATCACGCCATACG‑TCTACGGCACATGAATACAACGCTGGATTAATCTTCTGTGATAGt                                                  >  2:1279813/1‑98 (MQ=255)
                                          tttACTGTTTACCAACTGGCTCAAATATTACGTCTACCAGATCACGCAATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCg                                                >  2:1028960/1‑100 (MQ=255)
                                            tACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGAt                                              <  1:3010645/100‑1 (MQ=255)
                                            tACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGAt                                              <  2:1649404/100‑1 (MQ=255)
                                              cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                                            >  2:1446784/1‑100 (MQ=255)
                                              cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                                            >  1:2850621/1‑100 (MQ=255)
                                              cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAATAGATCg                                            >  1:1338753/1‑100 (MQ=255)
                                               tGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGt                                           <  2:356356/100‑1 (MQ=255)
                                                gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGtt                                          <  1:1505766/100‑1 (MQ=255)
                                                   tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTatg                                       <  1:2648477/100‑3 (MQ=255)
                                                    aCCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTaa                                        >  1:1265224/1‑98 (MQ=255)
                                                    aCCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc                                      <  2:1574358/100‑1 (MQ=255)
                                                    aCCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc                                      <  2:3361703/100‑1 (MQ=255)
                                                     ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCg                                     >  1:1753865/1‑100 (MQ=255)
                                                     ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCg                                     >  1:2450281/1‑100 (MQ=255)
                                                     ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCg                                     >  1:1282706/1‑100 (MQ=255)
                                                       aaCTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGAt                                   >  1:2391985/1‑100 (MQ=255)
                                                         cTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTc                                 <  1:758588/100‑1 (MQ=255)
                                                         cTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTc                                 <  1:3766938/100‑1 (MQ=255)
                                                           ggCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg                               >  1:1896822/1‑100 (MQ=255)
                                                           ggCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg                               >  1:2810451/1‑100 (MQ=255)
                                                           ggCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg                               >  2:591089/1‑100 (MQ=255)
                                                           ggCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGACACATGAATCCAACGCTGTATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTc                                 <  2:1888960/98‑1 (MQ=255)
                                                           ggCTCAACTATTACGTCTAACAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg                               >  2:1404987/1‑100 (MQ=255)
                                                            gCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGc                              <  2:2713095/100‑1 (MQ=255)
                                                            gCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGc                              <  2:2439149/100‑1 (MQ=255)
                                                             cTCAACTATTACGTCTACTATATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCa                             <  1:1937629/100‑1 (MQ=255)
                                                             cTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCa                             >  1:1171590/1‑100 (MQ=255)
                                                             cTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCa                             >  1:683999/1‑100 (MQ=255)
                                                             cTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCa                             >  1:847542/1‑100 (MQ=255)
                                                             cTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCa                             >  2:311998/1‑100 (MQ=255)
                                                             cTCAACTATTACGTCTACCAGATAACGCCATACGA‑CTCCGGCACATGAATCCAACGCTGGATTAATCTTATGTAATAGTCGATAGTTAAACGATTCAGc                              >  2:876456/1‑99 (MQ=255)
                                                                aacgagtACGACGACCAGATCAGGCGAGA‑GGTGTACGGGACATGAATGCAACGCGGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAcc                           <  2:3004459/93‑1 (MQ=255)
                                                                aaCTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCt                          >  2:2187978/1‑100 (MQ=255)
                                                                    aTTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTAc                       >  1:1357984/1‑99 (MQ=255)
                                                                      tACGTCTACCAGATCACGCTATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTc                    <  1:3581507/100‑1 (MQ=255)
                                                                      tACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTc                    <  2:1053072/100‑1 (MQ=255)
                                                                       aCGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCa                   >  2:786716/1‑100 (MQ=255)
                                                                          tCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGc                <  2:543090/100‑1 (MQ=255)
                                                                           cTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCa               >  1:2497927/1‑100 (MQ=255)
                                                                           cTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCa               >  2:527499/1‑100 (MQ=255)
                                                                            tACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAc              >  1:3898852/1‑100 (MQ=255)
                                                                            tACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAc              >  1:392755/1‑100 (MQ=255)
                                                                            tACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAc              >  2:2356947/1‑100 (MQ=255)
                                                                              ccAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCt            >  2:3548894/1‑100 (MQ=255)
                                                                                aGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTc          >  2:222510/1‑100 (MQ=255)
                                                                                aGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTc          <  2:3709349/100‑1 (MQ=255)
                                                                                  aTCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCgg        >  1:1572637/1‑100 (MQ=255)
                                                                                   tCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCggg       <  2:1211136/100‑1 (MQ=255)
                                                                                    cACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGt      <  2:392740/100‑1 (MQ=255)
                                                                                    cACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGt      >  1:3289803/1‑100 (MQ=255)
                                                                                    cACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGt      >  2:3061055/1‑100 (MQ=255)
                                                                                     aCGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTg     <  2:385940/100‑1 (MQ=255)
                                                                                      cGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGCTACTCGATCGTTCAGCGATTCACCACCTTACTTCACGCACCTTCgg        >  2:3805378/1‑96 (MQ=255)
                                                                                      cGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGTTTCATCTCCTTACCTCAGGCACCTTCGGGTg     <  2:3414045/99‑1 (MQ=255)
                                                                                      cGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTg     <  1:2892380/99‑1 (MQ=255)
                                                                                      cGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGc    >  2:3248313/1‑100 (MQ=255)
                                                                                      cGCCATACG‑TCTACGGCACATAAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGCTAAGCGATTCCGAACCTTAACTCAGGCATCTTCggg       >  2:3622649/1‑97 (MQ=255)
                                                                                       gCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGcc   <  1:1953712/100‑1 (MQ=255)
                                                                                       gCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGcc   <  2:660700/100‑1 (MQ=255)
                                                                                        ccATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCt  <  2:3912126/100‑1 (MQ=255)
                                                                                               |                                                                                             
TCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCT  >  NC_000913/4215083‑4215271

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