Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,095,111 G→A intergenic (‑82/‑206) rhaB ← / → rhaS rhamnulokinase/RhaS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,095,1110GA100.0% 946.9 / NA 252intergenic (‑82/‑206)rhaB/rhaSrhamnulokinase/RhaS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (139/113);  total (139/113)

Only 100 of 255 total aligned reads displayed.

AAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAAT  >  NC_000913/4095021‑4095200
                                                                                          |                                                                                         
agtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGc                                                                                         <  2:470887‑M2/81‑1 (MQ=255)
 gtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGcc                                                                                        >  1:892889‑M2/12‑93 (MQ=255)
  tcgacatatgGTGATCCTGCTTAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCa                                                                                       <  1:4911323‑M2/83‑1 (MQ=255)
   cgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAt                                                                                      <  2:5334415‑M2/84‑1 (MQ=255)
   cgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAt                                                                                      <  2:2878496‑M2/84‑1 (MQ=255)
    gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAtt                                                                                     <  1:2065837‑M2/85‑1 (MQ=255)
    gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAtt                                                                                     >  2:13369‑M2/9‑93 (MQ=255)
     acatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTg                                                                                    >  1:5530097‑M2/8‑93 (MQ=255)
     acatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTg                                                                                    <  1:441237‑M2/86‑1 (MQ=255)
      catatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTgg                                                                                   <  2:227823‑M2/87‑1 (MQ=255)
       atatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGc                                                                                  >  2:1474704‑M2/6‑93 (MQ=255)
        tatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCCTTTGCa                                                                                 >  2:6071338‑M2/5‑93 (MQ=255)
        tatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCa                                                                                 >  1:5844338‑M2/5‑93 (MQ=255)
        tatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCa                                                                                 >  2:6001189‑M2/5‑93 (MQ=255)
        tatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCa                                                                                 >  1:2002752‑M2/5‑93 (MQ=255)
           gGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAAcc                                                                              <  2:1963178‑M2/92‑1 (MQ=255)
            gtgATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCa                                                                             >  2:5801494/1‑93 (MQ=255)
            gtgATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCa                                                                             >  1:4555178/1‑93 (MQ=255)
            gtgATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCa                                                                             <  2:85700/93‑1 (MQ=255)
             tgATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAg                                                                            <  2:2276353/93‑1 (MQ=255)
              gATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAgg                                                                           <  1:4601372/93‑1 (MQ=255)
               aTCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAggg                                                                          <  1:5321585/93‑1 (MQ=255)
                tCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGa                                                                         <  1:5553555/93‑1 (MQ=255)
                tCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGa                                                                         >  2:3724604/1‑93 (MQ=255)
                 ccTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGaa                                                                        >  1:276142/1‑93 (MQ=255)
                  ctgctgAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGaaa                                                                       >  2:4740209/1‑93 (MQ=255)
                  ctgctgAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGaaa                                                                       >  1:4636534/1‑93 (MQ=255)
                   tgctgAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCCTTGGCAACCAGGGaaaa                                                                      >  2:2392570/1‑92 (MQ=255)
                        aaTTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGaa                                                                 <  1:5003413/93‑1 (MQ=255)
                         aTTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAAc                                                                <  2:1925924/93‑1 (MQ=255)
                          tttCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACg                                                               <  2:1894129/93‑1 (MQ=255)
                             cATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGtga                                                            >  1:329205/1‑93 (MQ=255)
                             cATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGtga                                                            >  1:1963646/1‑93 (MQ=255)
                              aTTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGtgat                                                           >  1:3790878/1‑93 (MQ=255)
                                   gACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGt                                                      >  2:3405448/1‑93 (MQ=255)
                                   gACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGt                                                      >  2:4431981/1‑93 (MQ=255)
                                   gACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGt                                                      >  1:2948391/1‑93 (MQ=255)
                                    aCCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGtt                                                     <  1:1055936/93‑1 (MQ=255)
                                     ccAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTc                                                    >  2:3272435/1‑93 (MQ=255)
                                     ccAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTc                                                    >  2:4017981/1‑93 (MQ=255)
                                     ccAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTc                                                    >  2:5578562/1‑93 (MQ=255)
                                     ccAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTc                                                    <  1:636548/93‑1 (MQ=255)
                                      cAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACAAGGGAAAGATGAACGTGATGATGTTca                                                   >  1:3974199/1‑93 (MQ=255)
                                        gTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTcaca                                                 >  2:4850610/1‑93 (MQ=255)
                                         tCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACaa                                                <  2:5380869/93‑1 (MQ=255)
                                         tCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACaa                                                >  1:1090859/1‑93 (MQ=255)
                                         tCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACaa                                                >  1:3494176/1‑93 (MQ=255)
                                         tCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACaa                                                >  1:5159126/1‑93 (MQ=255)
                                          cTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAAt                                               <  2:5249314/93‑1 (MQ=255)
                                            aaaaaGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAAttt                                             <  2:2134058/93‑1 (MQ=255)
                                             aaaaGCGCCTGAATTCGCGACCTTCTCGTTACCGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTg                                            <  2:3416302/93‑1 (MQ=255)
                                              aaaGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGc                                           <  1:728660/93‑1 (MQ=255)
                                                aGCGCCTGAATTCGCGAGCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTg                                         >  2:3371243/1‑93 (MQ=255)
                                                aGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTg                                         >  1:2660828/1‑93 (MQ=255)
                                                   gcCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAAt                                      >  1:2323557/1‑93 (MQ=255)
                                                   gcCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAAt                                      <  2:2758925/93‑1 (MQ=255)
                                                    ccTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGCGATGATGTTCACAATTTGCTGaaat                                     >  2:1129623/1‑91 (MQ=255)
                                                    ccTGAATTAGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGTAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAAtt                                     >  2:1043582/1‑93 (MQ=255)
                                                     cTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATtg                                    >  1:3569114/1‑93 (MQ=255)
                                                     cTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATtg                                    >  1:3890620/1‑93 (MQ=255)
                                                      tGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATtgt                                   >  2:4310008/1‑93 (MQ=255)
                                                      tGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATtgt                                   <  2:4390381/93‑1 (MQ=255)
                                                      tGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATtgt                                   <  1:370648/93‑1 (MQ=255)
                                                      tGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATtgt                                   >  2:4616227/1‑93 (MQ=255)
                                                      tGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATtgt                                   >  1:2402650/1‑93 (MQ=255)
                                                      tGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGGTCACAATTTGCTGAATtgt                                   >  1:2467252/1‑93 (MQ=255)
                                                       gAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATtgtg                                  <  1:613183/93‑1 (MQ=255)
                                                       gAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATtgtg                                  >  2:1506763/1‑93 (MQ=255)
                                                        aaTTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTgtgg                                 <  2:2449458/93‑1 (MQ=255)
                                                         aTTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTgtggt                                <  1:4214352/93‑1 (MQ=255)
                                                          ttCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTgtggtg                               >  1:5911520/1‑93 (MQ=255)
                                                          ttCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTgtggtg                               >  1:4026813/1‑93 (MQ=255)
                                                           tCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtga                              <  2:6238303/93‑1 (MQ=255)
                                                            cgcgACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtgat                             <  1:47543/93‑1 (MQ=255)
                                                            cgcgACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtgat                             <  2:165554/93‑1 (MQ=255)
                                                              cgACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtgatgt                           >  2:1815207/1‑93 (MQ=255)
                                                              cgACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtgatgt                           <  1:3553794/93‑1 (MQ=255)
                                                              cgACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtgatgt                           <  2:3817235/93‑1 (MQ=255)
                                                               gACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtgatgtg                          >  2:710185/1‑93 (MQ=255)
                                                                aCCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtgatgtga                         <  1:3427663/93‑1 (MQ=255)
                                                                  cTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGtgatgtgatg                       <  2:4136671/93‑1 (MQ=255)
                                                                  cTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGtgatgtgatg                       >  1:3347604/1‑93 (MQ=255)
                                                                   ttCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGc                      >  1:6125778/1‑93 (MQ=255)
                                                                      tcGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCa                   >  1:2183099/1‑93 (MQ=255)
                                                                      tcGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCa                   >  1:319023/1‑93 (MQ=255)
                                                                      tcGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCa                   <  2:2520786/93‑1 (MQ=255)
                                                                          tACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGc               >  2:5901181/1‑93 (MQ=255)
                                                                          tACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGc               >  2:5432327/1‑93 (MQ=255)
                                                                           aCTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCa              >  2:5151339/1‑93 (MQ=255)
                                                                            cTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCAt             >  2:2608202/1‑93 (MQ=255)
                                                                              gACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCAttt           >  1:1387026/1‑93 (MQ=255)
                                                                              gACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCAttt           >  1:2720411/1‑93 (MQ=255)
                                                                              gACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCAttt           <  2:307679/93‑1 (MQ=255)
                                                                                  ggAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTg       <  2:2033813/93‑1 (MQ=255)
                                                                                     aaaTGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGaaa    >  1:4721119/1‑93 (MQ=255)
                                                                                      aaTGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGaaaa   <  2:5540113/93‑1 (MQ=255)
                                                                                      aaTGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGaaaa   <  2:1562527/93‑1 (MQ=255)
                                                                                      aaTGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGaaaa   >  1:4286868/1‑93 (MQ=255)
                                                                                       aTGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAAt  >  1:5700187/1‑93 (MQ=255)
                                                                                       aTGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAAt  <  1:1474704/93‑1 (MQ=255)
                                                                                          |                                                                                         
AAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAAT  >  NC_000913/4095021‑4095200

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