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
                                                                                               |                                                                                             
tCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTAc                                                                                          >  2:204673/1‑100 (MQ=255)
 ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACg                                                                                         <  2:3902042/100‑1 (MQ=255)
  cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACgg                                                                                        <  1:325604/100‑1 (MQ=255)
   gcgaaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGc                                                                                       <  1:1030638/97‑1 (MQ=255)
   gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGc                                                                                       >  1:3057876/1‑100 (MQ=255)
    cAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGca                                                                                      <  1:11016/100‑1 (MQ=255)
      aaaTACACCGCTAGCGAGCTGGCGTAGGCACTGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGcaca                                                                                    <  1:3273320/100‑1 (MQ=255)
      aaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGcaca                                                                                    >  1:3096624/1‑100 (MQ=255)
       aaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGcaca                                                                                    <  1:1427293/100‑1 (MQ=255)
        aTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATg                                                                                  >  1:3374347/1‑100 (MQ=255)
        aTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATg                                                                                  <  1:353024/100‑1 (MQ=255)
          acacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGaa                                                                                >  1:3417582/1‑100 (MQ=255)
          acacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTAACAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACAGCACATGaa                                                                                >  2:2579006/1‑100 (MQ=255)
           cacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACAGCACATTAAt                                                                               >  2:2916580/1‑100 (MQ=255)
             ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATcc                                                                             >  1:1959853/1‑100 (MQ=255)
             ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATcc                                                                             >  1:2033613/1‑100 (MQ=255)
                     gagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTgg                                                                     >  1:451732/1‑100 (MQ=255)
                      agcTGGCGTAGTCACGGTAATTTACTGTTTACAAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTCCGGCACATGAATCCAAAGCTGGa                                                                    >  2:239724/1‑100 (MQ=255)
                       gcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGAt                                                                   <  2:2011719/100‑1 (MQ=255)
                        cTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTAAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGAtt                                                                  >  1:3860159/1‑100 (MQ=255)
                          ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTaa                                                                >  1:2104165/1‑100 (MQ=255)
                          ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACACCATACGTCTACGGCACATGAATCCAACGCTGGATTaa                                                                >  1:36749/1‑100 (MQ=255)
                               aGTCACGGTAATTTACTGTTTACCGACTGGCTTAACTATGAAGTCTACCAGATCACTCCATACGTCTACGGCACATGATTCCAACGCTGTATTAAtcttc                                                           <  2:3862041/100‑1 (MQ=255)
                               aGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtcttc                                                           <  2:2178832/100‑1 (MQ=255)
                                 tCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                                                         <  2:1805670/100‑1 (MQ=255)
                                 tCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                                                         <  1:3563366/100‑1 (MQ=255)
                                 tCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                                                         <  1:1619368/100‑1 (MQ=255)
                                 tCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                                                         <  1:2955510/100‑1 (MQ=255)
                                 tCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                                                         <  1:2059806/100‑1 (MQ=255)
                                  cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                        >  2:1505751/1‑100 (MQ=255)
                                   aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgtg                                                       <  1:114816/100‑1 (MQ=255)
                                   aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                        <  2:526269/99‑1 (MQ=255)
                                      gTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATa                                                    <  2:763650/100‑1 (MQ=255)
                                       tAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAg                                                   <  2:282262/100‑1 (MQ=255)
                                        aaTTTACTGTTTTCCAACTGGCTCAACTATTACTTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCGGTGAt                                                     <  2:483783/97‑1 (MQ=255)
                                        aaTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGt                                                  >  1:2479565/1‑100 (MQ=255)
                                         aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTc                                                 <  2:1552280/100‑1 (MQ=255)
                                         aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATAAATACAACGCTGGATTAATATTCTGTGATAGTc                                                 >  1:1904298/1‑100 (MQ=255)
                                         aTTTACTGTTTACCAACTGGATCAACTATTACGTCTACCAGAAAAAGCCATACGTCTACGGCACATGAATCCAAAGCTGGATTAATCTTCTGTGATAGTc                                                 >  2:1148586/1‑100 (MQ=255)
                                          tttACTGTTTACCAACTGGCTCAACTATTACGACTACCAGATCACGCCATACGTCTACGGCACATGAATACAACGCTGGATTAATCTTCTGTGATAGt                                                  >  2:1279813/1‑98 (MQ=255)
                                          tttACTGTTTACCAACTGGCTCAAATATTACGTCTACCAGATCACGCAATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCg                                                >  2:1028960/1‑100 (MQ=255)
                                            tACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGAt                                              <  1:3010645/100‑1 (MQ=255)
                                            tACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGAt                                              <  2:1649404/100‑1 (MQ=255)
                                              cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                                            >  2:1446784/1‑100 (MQ=255)
                                              cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                                            >  1:2850621/1‑100 (MQ=255)
                                              cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAATAGATCg                                            >  1:1338753/1‑100 (MQ=255)
                                               tGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGt                                           <  2:356356/100‑1 (MQ=255)
                                                gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGtt                                          <  1:1505766/100‑1 (MQ=255)
                                                   tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTatg                                       <  1:2648477/100‑3 (MQ=255)
                                                    aCCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTaa                                        >  1:1265224/1‑98 (MQ=255)
                                                    aCCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc                                      <  2:1574358/100‑1 (MQ=255)
                                                    aCCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc                                      <  2:3361703/100‑1 (MQ=255)
                                                     ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCg                                     >  1:1753865/1‑100 (MQ=255)
                                                     ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCg                                     >  1:2450281/1‑100 (MQ=255)
                                                     ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCg                                     >  1:1282706/1‑100 (MQ=255)
                                                       aaCTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGAt                                   >  1:2391985/1‑100 (MQ=255)
                                                         cTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTc                                 <  1:758588/100‑1 (MQ=255)
                                                         cTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTc                                 <  1:3766938/100‑1 (MQ=255)
                                                           ggCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg                               >  1:1896822/1‑100 (MQ=255)
                                                           ggCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg                               >  1:2810451/1‑100 (MQ=255)
                                                           ggCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg                               >  2:591089/1‑100 (MQ=255)
                                                           ggCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGACACATGAATCCAACGCTGTATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTc                                 <  2:1888960/98‑1 (MQ=255)
                                                           ggCTCAACTATTACGTCTAACAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg                               >  2:1404987/1‑100 (MQ=255)
                                                            gCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGc                              <  2:2713095/100‑1 (MQ=255)
                                                            gCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGc                              <  2:2439149/100‑1 (MQ=255)
                                                             cTCAACTATTACGTCTACTATATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCa                             <  1:1937629/100‑1 (MQ=255)
                                                             cTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCa                             >  1:1171590/1‑100 (MQ=255)
                                                             cTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCa                             >  1:683999/1‑100 (MQ=255)
                                                             cTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCa                             >  1:847542/1‑100 (MQ=255)
                                                             cTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCa                             >  2:311998/1‑100 (MQ=255)
                                                             cTCAACTATTACGTCTACCAGATAACGCCATACGACTCCGGCACATGAATCCAACGCTGGATTAATCTTATGTAATAGTCGATAGTTAAACGATTCAGc                              >  2:876456/1‑99 (MQ=255)
                                                                aacgagtACGACGACCAGATCAGGCGAGAGGTGTACGGGACATGAATGCAACGCGGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAcc                           <  2:3004459/93‑1 (MQ=255)
                                                                aaCTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCt                          >  2:2187978/1‑100 (MQ=255)
                                                                    aTTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTAc                       >  1:1357984/1‑99 (MQ=255)
                                                                      tACGTCTACCAGATCACGCTATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTc                    <  1:3581507/100‑1 (MQ=255)
                                                                      tACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTc                    <  2:1053072/100‑1 (MQ=255)
                                                                       aCGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCa                   >  2:786716/1‑100 (MQ=255)
                                                                          tCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGc                <  2:543090/100‑1 (MQ=255)
                                                                           cTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCa               >  1:2497927/1‑100 (MQ=255)
                                                                           cTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCa               >  2:527499/1‑100 (MQ=255)
                                                                            tACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAc              >  1:3898852/1‑100 (MQ=255)
                                                                            tACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAc              >  1:392755/1‑100 (MQ=255)
                                                                            tACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAc              >  2:2356947/1‑100 (MQ=255)
                                                                              ccAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCt            >  2:3548894/1‑100 (MQ=255)
                                                                                aGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTc          >  2:222510/1‑100 (MQ=255)
                                                                                aGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTc          <  2:3709349/100‑1 (MQ=255)
                                                                                  aTCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCgg        >  1:1572637/1‑100 (MQ=255)
                                                                                   tCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCggg       <  2:1211136/100‑1 (MQ=255)
                                                                                    cACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGt      <  2:392740/100‑1 (MQ=255)
                                                                                    cACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGt      >  1:3289803/1‑100 (MQ=255)
                                                                                    cACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGt      >  2:3061055/1‑100 (MQ=255)
                                                                                     aCGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTg     <  2:385940/100‑1 (MQ=255)
                                                                                      cGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGCTACTCGATCGTTCAGCGATTCACCACCTTACTTCACGCACCTTCgg        >  2:3805378/1‑96 (MQ=255)
                                                                                      cGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGTTTCATCTCCTTACCTCAGGCACCTTCGGGTg     <  2:3414045/99‑1 (MQ=255)
                                                                                      cGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTg     <  1:2892380/99‑1 (MQ=255)
                                                                                      cGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGc    >  2:3248313/1‑100 (MQ=255)
                                                                                      cGCCATACGTCTACGGCACATAAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGCTAAGCGATTCCGAACCTTAACTCAGGCATCTTCggg       >  2:3622649/1‑97 (MQ=255)
                                                                                       gCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGcc   <  1:1953712/100‑1 (MQ=255)
                                                                                       gCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGcc   <  2:660700/100‑1 (MQ=255)
                                                                                        ccATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCt  <  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: