Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,808,078 A→G intergenic (+6/‑157) yniC → / → b1728 putative phosphatase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,808,0780AG100.0% 875.2 / NA 268intergenic (+6/‑157)yniC/b1728putative phosphatase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (119/149);  total (119/149)

Only 100 of 283 total aligned reads displayed.

AGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACAC  >  NC_000913/1807988‑1808166
                                                                                          |                                                                                        
aGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCggg                                                                                        >  1:6647874/1‑93 (MQ=255)
 gcgcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCgggg                                                                                       <  2:5347313/93‑1 (MQ=255)
  cgcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGgggg                                                                                      >  1:149665/1‑93 (MQ=255)
    cAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCa                                                                                    >  1:6529343/1‑93 (MQ=255)
     aaaaTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAg                                                                                   >  1:2125927/1‑93 (MQ=255)
      aaaTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGt                                                                                  >  1:2693267/1‑93 (MQ=255)
      aaaTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCAGTTAATGACCGGGGGCAGt                                                                                  >  2:529604/1‑93 (MQ=255)
        aTGATCCACGTTTTGTAGTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGa                                                                                <  1:2901846/93‑1 (MQ=255)
          gATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGAtc                                                                              >  1:4724241/1‑93 (MQ=255)
           aTCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGAtct                                                                             >  2:2997771/1‑93 (MQ=255)
            tCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGAtctc                                                                            >  2:2876418/1‑93 (MQ=255)
              cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGc                                                                          <  2:2433586/93‑1 (MQ=255)
              cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGc                                                                          <  1:1334725/93‑1 (MQ=255)
              cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGc                                                                          <  1:3975397/93‑1 (MQ=255)
              cACGTTGTGTATTAGCAGAAGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGc                                                                          <  1:6241132/93‑1 (MQ=255)
               aCGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCt                                                                         >  2:3159398/1‑93 (MQ=255)
                cGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTg                                                                        >  2:5809244/1‑93 (MQ=255)
                 gTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAACAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGc                                                                       >  1:2272734/1‑93 (MQ=255)
                 gTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGc                                                                       >  1:1511637/1‑93 (MQ=255)
                   tttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGccc                                                                     <  1:5226555/93‑1 (MQ=255)
                    ttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGcccc                                                                    <  2:2046667/93‑1 (MQ=255)
                      gTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTg                                                                  <  2:3745354/93‑1 (MQ=255)
                      gTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTg                                                                  <  2:4495074/93‑1 (MQ=255)
                       tATTAGCAGACGTCAAACTTTCATCGGTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTgg                                                                 >  1:3463051/1‑93 (MQ=255)
                        aTTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGt                                                                <  1:1681893/93‑1 (MQ=255)
                        aTTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGt                                                                >  2:6047155/1‑93 (MQ=255)
                        aTTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGt                                                                <  1:4842873/93‑1 (MQ=255)
                         ttAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGtt                                                               <  2:3407407/93‑1 (MQ=255)
                         ttAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGtt                                                               <  1:5450875/93‑1 (MQ=255)
                         ttAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGtt                                                               >  1:563408/1‑93 (MQ=255)
                          tAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTc                                                              <  1:4301866/93‑1 (MQ=255)
                           aGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCt                                                             >  1:749556/1‑93 (MQ=255)
                              aGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTa                                                          <  2:1883915/93‑1 (MQ=255)
                              aGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTa                                                          <  2:6204856/93‑1 (MQ=255)
                              aGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTa                                                          <  1:258668/93‑1 (MQ=255)
                              aGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTa                                                          <  2:5652372/93‑1 (MQ=255)
                               gACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTAt                                                         >  2:4982253/1‑93 (MQ=255)
                                aCGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATc                                                        >  1:2354543/1‑93 (MQ=255)
                                 cGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCt                                                       >  2:1556565/1‑93 (MQ=255)
                                    cAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGaa                                                    <  1:3407544/93‑1 (MQ=255)
                                     aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGaaa                                                   >  2:2236727/1‑93 (MQ=255)
                                     aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGaaa                                                   <  1:6106611/93‑1 (MQ=255)
                                      aaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAt                                                  <  1:5774650/93‑1 (MQ=255)
                                      aaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAt                                                  >  1:6040750/1‑93 (MQ=255)
                                        cTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTg                                                >  2:1676364/1‑93 (MQ=255)
                                         tttCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGc                                               <  1:5711940/93‑1 (MQ=255)
                                             aTCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTc                                           >  1:3707672/1‑93 (MQ=255)
                                             aTCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTc                                           <  2:4798919/93‑1 (MQ=255)
                                                gCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAAc                                        <  1:108896/93‑1 (MQ=255)
                                                gCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAAc                                        >  1:1850198/1‑93 (MQ=255)
                                                   gACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGa                                     <  2:4461128/93‑1 (MQ=255)
                                                   gACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGa                                     >  1:4547012/1‑93 (MQ=255)
                                                   gACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGa                                     >  2:1726144/1‑93 (MQ=255)
                                                    aCAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGAc                                    <  1:3348394/93‑1 (MQ=255)
                                                     cAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACg                                   >  2:5241911/1‑93 (MQ=255)
                                                      aGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACgg                                  >  2:1809514/1‑93 (MQ=255)
                                                        aaCTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGAt                                <  2:2926777/93‑1 (MQ=255)
                                                         aCTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGAtt                               >  1:4176656/1‑93 (MQ=255)
                                                           tCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTaa                             <  1:5848875/93‑1 (MQ=255)
                                                           tCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTaa                             >  2:549877/1‑93 (MQ=255)
                                                            cACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAAt                            <  2:4625843/93‑1 (MQ=255)
                                                             aCCGCAAAAGACCTTCTCGTTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATc                           <  1:4426500/93‑1 (MQ=255)
                                                             aCCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCGTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATc                           <  2:2217961/93‑1 (MQ=255)
                                                             aCCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATc                           <  1:1179506/93‑1 (MQ=255)
                                                             aCCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATc                           <  1:3237478/93‑1 (MQ=255)
                                                             aCCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATc                           <  1:5058070/93‑1 (MQ=255)
                                                             aCCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATc                           <  1:1954126/93‑1 (MQ=255)
                                                              ccGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCg                          >  2:3653457/1‑93 (MQ=255)
                                                              ccGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCg                          >  2:3889011/1‑93 (MQ=255)
                                                                gCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTc                        >  2:6568067/1‑93 (MQ=255)
                                                                gCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTc                        >  1:6421128/1‑93 (MQ=255)
                                                                gCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTc                        <  1:5179360/93‑1 (MQ=255)
                                                                gCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTc                        <  2:5698715/93‑1 (MQ=255)
                                                                 cAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCa                       <  1:3282044/93‑1 (MQ=255)
                                                                  aaaaGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGCTTAATCGTCat                      >  1:1800289/1‑92 (MQ=255)
                                                                  aaaaGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCaa                      <  1:791399/93‑1 (MQ=255)
                                                                  aaaaGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCaa                      >  2:5165416/1‑93 (MQ=255)
                                                                   aaaGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAt                     <  1:6607442/93‑1 (MQ=255)
                                                                     aGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAttt                   >  1:991211/1‑93 (MQ=255)
                                                                      gACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTa                  <  2:1169045/93‑1 (MQ=255)
                                                                       aCCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTaa                 <  2:3259908/93‑1 (MQ=255)
                                                                        ccTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAag                >  2:2611194/1‑93 (MQ=255)
                                                                          ttCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAagag              <  2:792482/93‑1 (MQ=255)
                                                                          ttCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAagag              >  1:3897105/1‑93 (MQ=255)
                                                                           tctcGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAgaga             <  2:2280332/93‑1 (MQ=255)
                                                                           tctcGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAgaga             >  1:6298286/1‑93 (MQ=255)
                                                                            ctcGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGaa            <  2:4181991/93‑1 (MQ=255)
                                                                            ctcGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGaa            <  1:3513557/93‑1 (MQ=255)
                                                                             tcGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGaaa           <  1:4943324/93‑1 (MQ=255)
                                                                               ggTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAaga         >  1:1006890/1‑93 (MQ=255)
                                                                               ggTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAaga         <  1:3606985/93‑1 (MQ=255)
                                                                                gTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAagag        >  2:4064112/1‑93 (MQ=255)
                                                                                gTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAagag        <  2:1397522/93‑1 (MQ=255)
                                                                                gTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAagag        <  2:1976811/93‑1 (MQ=255)
                                                                                 ttAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGt       <  1:1223782/93‑1 (MQ=255)
                                                                                 ttAATGACCGGGGGCAGCGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGt       >  2:1000422/1‑93 (MQ=255)
                                                                                  tAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGtt      <  2:4928948/93‑1 (MQ=255)
                                                                                  tAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGtt      <  2:1336633/93‑1 (MQ=255)
                                                                                      gACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTacac  >  2:3708809/1‑93 (MQ=255)
                                                                                      gACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTacac  <  2:6498746/93‑1 (MQ=255)
                                                                                          |                                                                                        
AGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACAC  >  NC_000913/1807988‑1808166

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