Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,215,122 T→G G281G (GGT→GGG metA → homoserine O‑transsuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,215,1220TG98.8% 584.6 / ‑4.7 173G281G (GGT→GGGmetAhomoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base G (82/87);  total (86/87)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 173 total aligned reads displayed.

GCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTG  >  NC_000913/4215025‑4215218
                                                                                                 |                                                                                                
gcgATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAGAATACACCGCGAGCGAGCTGGCGTAGTCACGGGaa                                                                                                >  2:3949190/1‑100 (MQ=255)
 cgATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGaag                                                                                               >  2:1574964/1‑99 (MQ=255)
 cgATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAAt                                                                                               <  2:417055/100‑1 (MQ=255)
  gATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAAtn                                                                                              <  1:3261380/100‑2 (MQ=255)
    tgtgGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTa                                                                                            <  1:1137999/100‑1 (MQ=255)
    tgtgGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTa                                                                                            >  2:332203/1‑100 (MQ=255)
     gtgGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTAc                                                                                           <  1:3312173/100‑1 (MQ=255)
       ggAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTg                                                                                         >  1:2241141/1‑100 (MQ=255)
       ggAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTg                                                                                         >  1:2335096/1‑100 (MQ=255)
       ggAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTg                                                                                         >  1:876881/1‑100 (MQ=255)
        gAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGt                                                                                        <  1:2848869/100‑1 (MQ=255)
          aGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGttt                                                                                      <  1:71961/100‑1 (MQ=255)
          aGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGttt                                                                                      >  1:1074818/1‑100 (MQ=255)
           gCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTa                                                                                     <  2:1262391/100‑1 (MQ=255)
            ccGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTAc                                                                                    <  2:1139779/100‑1 (MQ=255)
             cGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTAcc                                                                                   <  2:269076/100‑1 (MQ=255)
             cGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTAcc                                                                                   <  1:3332003/100‑1 (MQ=255)
             cGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTAcc                                                                                   >  1:3732481/1‑100 (MQ=255)
              ggACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCa                                                                                  >  1:846061/1‑100 (MQ=255)
                aCTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAAc                                                                                >  1:2928513/1‑100 (MQ=255)
                 cTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACt                                                                               <  1:1898677/100‑1 (MQ=255)
                 cTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACt                                                                               <  2:2603623/100‑1 (MQ=255)
                 cTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCAAGCGAGCTGGCGTAGTCACGGGACTTTACTGTTTACCAACt                                                                               >  2:2936352/1‑100 (MQ=255)
                  tAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTg                                                                              >  1:1480596/1‑100 (MQ=255)
                  tAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTg                                                                              <  1:2910641/100‑1 (MQ=255)
                   aGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTgg                                                                             >  1:1307945/1‑100 (MQ=255)
                    gACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGc                                                                            <  2:2991271/100‑1 (MQ=255)
                    gACCCGGATGTACCGTATAACTATTTCCCGCAAAATGATCCGCAAAATTCACCGCGAGCGCGCTGGAGTAGTTACGGGAATTTACTGTTTACCAACTGGc                                                                            <  2:150845/100‑1 (MQ=255)
                     aCCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACGGGCt                                                                           <  2:3403065/100‑1 (MQ=255)
                      cccGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTc                                                                          <  2:1743753/100‑1 (MQ=255)
                      cccGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTc                                                                          >  2:1187863/1‑100 (MQ=255)
                      cccGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTc                                                                          >  2:625420/1‑100 (MQ=255)
                       ccGGATGTGCCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCa                                                                         <  1:1427617/100‑1 (MQ=255)
                       ccGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCa                                                                         <  1:638708/100‑1 (MQ=255)
                          gATGTAGCGTATAACTATTTCCCGAACGATGACCCGCACAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACt                                                                      <  2:3998020/100‑1 (MQ=255)
                          gATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACt                                                                      <  2:1300296/100‑1 (MQ=255)
                             gTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTAtt                                                                   >  1:511920/1‑100 (MQ=255)
                              tACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTa                                                                  >  1:1041145/1‑100 (MQ=255)
                              tACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTa                                                                  <  2:2503070/100‑1 (MQ=255)
                              tACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTa                                                                  <  1:1234939/100‑1 (MQ=255)
                               aCCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTAc                                                                 >  1:1609322/1‑100 (MQ=255)
                                ccGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACg                                                                >  1:2421774/1‑100 (MQ=255)
                                   tataAGTATTTGCCGCACAGTGATCCGCAATATACACCGCGAGTGTGCTGGCGTAGTCAGGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCt                                                             <  2:3226361/100‑1 (MQ=255)
                                     taACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTAc                                                           >  1:1153312/1‑100 (MQ=255)
                                     taACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTAc                                                           >  1:1038472/1‑100 (MQ=255)
                                      aaCTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTAcc                                                          >  2:2673588/1‑100 (MQ=255)
                                        cTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAg                                                        <  2:2954338/100‑1 (MQ=255)
                                        cTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAg                                                        >  2:3471934/1‑100 (MQ=255)
                                        cTAGTTCCCGCACAATGATCCGCAAAATACAACGCTAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAg                                                        <  2:2870331/100‑1 (MQ=255)
                                         tATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGa                                                       <  2:3325833/100‑1 (MQ=255)
                                          aTTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGAt                                                      >  1:936166/1‑100 (MQ=255)
                                           tttCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATc                                                     <  1:1285953/100‑1 (MQ=255)
                                           tttCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGACTCAACTATTACGTCTACCatatg                                                     >  1:3028923/1‑96 (MQ=255)
                                             tCCCGCACAATGAGCCGCAAAATACACCGCGCGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCAc                                                   <  2:3136569/100‑1 (MQ=255)
                                              cccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACg                                                  <  2:56803/100‑1 (MQ=255)
                                              cccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGGCACGGGAAATTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACg                                                  <  2:644767/100‑1 (MQ=255)
                                               ccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCActc                                                 >  1:830375/1‑98 (MQ=255)
                                               ccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGc                                                 >  2:2135126/1‑100 (MQ=255)
                                                cGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGcc                                                >  2:2301565/1‑100 (MQ=255)
                                                cGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGcc                                                >  2:3562199/1‑100 (MQ=255)
                                                 gCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCa                                               <  1:2771186/100‑1 (MQ=255)
                                                   acaATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCCGATCACGCCATa                                             >  2:1116177/1‑100 (MQ=255)
                                                   acaATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATa                                             <  1:3447451/100‑1 (MQ=255)
                                                    caATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCCCAACTATTACGTCTACCAGATCACTCCATAc                                            >  1:1322079/1‑100 (MQ=255)
                                                      aTGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTTTTACGTCTACCAGATCACGCCATACGa                                          >  2:3861363/1‑100 (MQ=255)
                                                       tGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAt                                         <  2:1474669/100‑1 (MQ=255)
                                                           ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTAc                                     >  2:2191719/1‑100 (MQ=255)
                                                            cGCACAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACg                                    <  2:1210249/100‑1 (MQ=255)
                                                            cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGAAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACg                                    >  2:714841/1‑100 (MQ=255)
                                                              cAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGc                                  <  2:2373153/100‑1 (MQ=255)
                                                              cAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGc                                  <  2:2233192/100‑1 (MQ=255)
                                                               aaaaTACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGca                                 >  1:634130/1‑100 (MQ=255)
                                                                aaaTACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGcac                                <  1:1789735/100‑1 (MQ=255)
                                                                aaaTACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGcac                                <  2:29043/100‑1 (MQ=255)
                                                                 aaTACACCGCGAGCGAGCTGGCGTCGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGACTACCAGATCAAGCCATACGATCTACGGcaca                               >  2:1829366/1‑100 (MQ=255)
                                                                  aTACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACAt                              <  1:1905907/100‑1 (MQ=255)
                                                                          cgagcgagCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCaa                      >  2:1339720/1‑100 (MQ=255)
                                                                           gagcgagcTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAAc                     <  1:1344488/100‑1 (MQ=255)
                                                                           gagcgagcTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAAc                     <  1:3056942/100‑1 (MQ=255)
                                                                           gagcgagcTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGAACACGCCATACGATCTACGGCACATGAATCCAAc                     >  1:4038730/1‑100 (MQ=255)
                                                                            agcgagcTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACg                    >  1:2538140/1‑100 (MQ=255)
                                                                               gagcTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTg                 >  2:2266345/1‑100 (MQ=255)
                                                                                agcTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCTAACGCttg                >  1:405879/1‑98 (MQ=255)
                                                                                agcTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTgg                >  1:2383592/1‑100 (MQ=255)
                                                                                 gcTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGa               <  1:2851899/100‑1 (MQ=255)
                                                                                  cTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGcgcgaa              >  1:3916894/1‑95 (MQ=255)
                                                                                  cTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGAt              <  1:3833529/100‑1 (MQ=255)
                                                                                  cTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGAt              <  2:2092003/100‑1 (MQ=255)
                                                                                   tGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGAtt             >  1:2819530/1‑100 (MQ=255)
                                                                                    ggCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTTAATTATTACGTCTAACAGATCACGCCATACGATCTACGGCACATGAATGCAACGTTGGATTa            >  1:583222/1‑100 (MQ=255)
                                                                                    ggCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCAGGCCATACGATCTACGGCACATGAATCCAACGCTGGAtca            >  1:2807700/1‑98 (MQ=255)
                                                                                       gTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtc         <  2:3504367/100‑1 (MQ=255)
                                                                                       gTAGTCACGAGAATTTACTGTTTGCCAACTGGCTCAACTATTACGTCTACCAGAACACGCCATACGATCGACGGCACATGAATCTAACGCTGGATTAAtc         <  1:3548716/100‑1 (MQ=255)
                                                                                            cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg    >  2:2799066/1‑100 (MQ=255)
                                                                                            cACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg    <  1:2641152/100‑1 (MQ=255)
                                                                                            cACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg    >  2:3001430/1‑100 (MQ=255)
                                                                                            cACGGCAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACACCATACGATCTACGGCACATGAATCCAACGCTGGCTTAATCTTCtg    >  2:673865/1‑100 (MQ=255)
                                                                                             aCGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt   <  1:382544/100‑1 (MQ=255)
                                                                                              cGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgtg  >  2:2414769/1‑100 (MQ=255)
                                                                                              cGGGAATCTACTGTTTACCAACTGGCCCCACTATTACGTCTACCAAATCACGCCATACGATCAACGGCACATCAATCCAACCCTGGCTCAAccttccgag  >  2:1417274/1‑96 (MQ=255)
                                                                                                 |                                                                                                
GCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTG  >  NC_000913/4215025‑4215218

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: