Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,889,795 G→A A184A (GCG→GCA yidY → putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,889,7950GA91.9% 809.5 / 39.4 246A184A (GCG→GCAyidYputative transport protein
Reads supporting (aligned to +/- strand):  ref base G (10/10);  new base A (122/104);  total (132/114)
Fisher's exact test for biased strand distribution p-value = 8.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

Only 100 of 265 total aligned reads displayed.

TTCCCGTGGCAGAGTCTGTTCTGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAG  >  NC_000913/3889703‑3889884
                                                                                            |                                                                                         
  cccGTGGCAGAGTCTGTTCTGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCagcagc                                                                                         <  2:907426/93‑1 (MQ=255)
  cccGTGGCAGAGTCTGTTCTGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCagcagc                                                                                         <  2:1434906/93‑1 (MQ=255)
  cccGTGGCAGAGTCTGTTCTGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCagcagc                                                                                         <  1:987692/93‑1 (MQ=255)
   ccGTGGCAGAGTCTGTTCTGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGcc                                                                                        >  1:272139/1‑93 (MQ=255)
     gTGGCAGAGTCTGTTCTGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGcccc                                                                                      <  1:900322/93‑1 (MQ=255)
     gTGGCAGAGTCTGTTCTGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGcccc                                                                                      >  2:547380/1‑93 (MQ=255)
     gTGGCAGAGTCTGTTCTGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGcccc                                                                                      <  1:5711925/93‑1 (MQ=255)
      tGGCAGAGTCTGTTCTGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGccccc                                                                                     >  2:6062235/1‑93 (MQ=255)
      tGGCAGAGTCTGTTCTGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGccccc                                                                                     <  2:5617870/93‑1 (MQ=255)
        gCAGAGTCTGTTCTGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGc                                                                                   >  2:5270528/1‑93 (MQ=255)
        gCAGAGTCTGTTCTGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGc                                                                                   <  1:6410870/93‑1 (MQ=255)
        gCAGAGTCTGTTCTGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGc                                                                                   >  1:4795929/1‑93 (MQ=255)
         cAGAGTCTGTTCTGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCa                                                                                  <  1:3351714/93‑1 (MQ=255)
          agagTCTGTTCTGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAg                                                                                 >  1:3159688/1‑93 (MQ=255)
           gagTCTGTTCTGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGc                                                                                >  1:6507525/1‑93 (MQ=255)
           gagTCTGTTCTGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGc                                                                                >  1:2468005/1‑93 (MQ=255)
             gTCTGTTCTGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCtt                                                                              <  2:5314034/93‑1 (MQ=255)
                 tttCTGTGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCGCGCAGCTTCGGa                                                                          <  2:1602216/92‑1 (MQ=255)
                 tttCTGGGTGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCGGCCCCCGCAGCTTCGGa                                                                          <  1:1267086/92‑1 (MQ=255)
                 gTTCTGGGTGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCGGCCCCCGCAGCTTCGGa                                                                          <  2:316879/93‑1 (MQ=255)
                   tCTGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATa                                                                        >  2:6287091/1‑93 (MQ=255)
                   tCTGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATa                                                                        <  1:1373316/93‑1 (MQ=255)
                     tGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATaaa                                                                      >  2:6619336/1‑93 (MQ=255)
                     tGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATaaa                                                                      >  1:5120793/1‑93 (MQ=255)
                         ggATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCAc                                                                  <  2:5256359/92‑1 (MQ=255)
                          ggtGGCAATTATGGGCATTGGGGTACTTAATGTGTTTTTTTTTAATTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACg                                                                 <  2:3031799/91‑1 (MQ=255)
                           gtGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACga                                                                <  1:1897363/92‑1 (MQ=255)
                           aTGGCAATGATGGGCATCGCGGTACTTATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACga                                                                <  1:3834948/93‑1 (MQ=255)
                           aTGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACga                                                                >  2:3889401/1‑93 (MQ=255)
                           aTGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACga                                                                <  1:2022630/93‑1 (MQ=255)
                             ggCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACgaga                                                              <  1:5943386/93‑1 (MQ=255)
                              gCAATGATGGGCATCGCGGTACTGCTGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGaa                                                             >  1:6403505/1‑93 (MQ=255)
                              gCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGaa                                                             >  1:3523145/1‑93 (MQ=255)
                                aaTGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGaaaa                                                           <  2:5957533/93‑1 (MQ=255)
                                 atgatgGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAAt                                                          <  1:1058696/93‑1 (MQ=255)
                                  tgatgGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATa                                                         >  1:1449131/1‑93 (MQ=255)
                                   gatgGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAACCGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAg                                                        <  1:5756782/93‑1 (MQ=255)
                                   gatgGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAg                                                        <  1:5917967/93‑1 (MQ=255)
                                    atgGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGc                                                       >  1:6269622/1‑93 (MQ=255)
                                        gCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGt                                                   <  1:4936713/93‑1 (MQ=255)
                                         cATCGCGGTACTGATGTTGCCTTTGTTTATTTTAAAAGAAACGTGCCCAGCATCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTc                                                  <  1:944206/93‑1 (MQ=255)
                                           tCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCgc                                                >  1:4840840/1‑93 (MQ=255)
                                           tCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCgc                                                <  2:1029257/93‑1 (MQ=255)
                                             gcgGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCgctg                                              >  2:5373613/1‑93 (MQ=255)
                                              cgGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCgctgc                                             <  1:2522596/93‑1 (MQ=255)
                                              cgGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCgctgc                                             <  2:5736068/93‑1 (MQ=255)
                                              cgGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCgctgc                                             >  2:2614020/1‑93 (MQ=255)
                                              cgGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCgctgc                                             >  2:2028468/1‑93 (MQ=255)
                                              cgGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCgctgc                                             <  1:1030689/93‑1 (MQ=255)
                                               ggTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCgctgct                                            <  1:2317408/93‑1 (MQ=255)
                                                gTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCtt                                           >  1:5554426/1‑93 (MQ=255)
                                                gTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCtt                                           >  1:5245329/1‑93 (MQ=255)
                                                 tACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTa                                          >  1:3660834/1‑93 (MQ=255)
                                                 tACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTa                                          <  2:1832444/93‑1 (MQ=255)
                                                  aCTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTaa                                         >  2:1246047/1‑93 (MQ=255)
                                                   cTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAAc                                        <  1:3572781/93‑1 (MQ=255)
                                                     gATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCg                                      >  2:841302/1‑93 (MQ=255)
                                                     gATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCg                                      <  2:4788844/93‑1 (MQ=255)
                                                     gATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCg                                      <  1:5305382/93‑1 (MQ=255)
                                                      aTGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGt                                     >  2:46999/1‑93 (MQ=255)
                                                       tgttgtCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCgct                                    >  2:3943046/1‑91 (MQ=255)
                                                       tgttgtCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGtt                                    >  1:4589432/1‑93 (MQ=255)
                                                        gttgtCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGttt                                   >  2:2422396/1‑93 (MQ=255)
                                                        gttgtCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGttt                                   >  2:6110205/1‑93 (MQ=255)
                                                        gttgtCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGttt                                   <  2:3274001/93‑1 (MQ=255)
                                                            tCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTc                               <  1:4433828/93‑1 (MQ=255)
                                                            tCTTTGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTc                               >  1:4696011/1‑93 (MQ=255)
                                                               ttGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCta                            >  2:3155263/1‑92 (MQ=255)
                                                               ttGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTc                            <  1:4326854/93‑1 (MQ=255)
                                                               ttGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTc                            <  1:5426214/93‑1 (MQ=255)
                                                               ttGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTc                            <  2:1944745/93‑1 (MQ=255)
                                                                tGTTTATTTTAAAAGAAACGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCa                           >  2:6419983/1‑93 (MQ=255)
                                                                tGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCa                           >  1:5070243/1‑93 (MQ=255)
                                                                tGTTTATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCa                           <  1:4726874/93‑1 (MQ=255)
                                                                  tttATTTTAAAAGAAACGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGc                         <  1:673581/93‑1 (MQ=255)
                                                                  tttATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGc                         >  2:2258684/1‑93 (MQ=255)
                                                                   ttATTTTAAAAGAAACGCGCCCAGCAGCTCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGcc                        >  2:4245774/1‑93 (MQ=255)
                                                                   ttATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGcc                        <  2:2753441/93‑1 (MQ=255)
                                                                    tATTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCg                       >  1:1606976/1‑93 (MQ=255)
                                                                     aTTTTAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGCTTTTTCCTCAGCCgt                      >  2:2862345/1‑93 (MQ=255)
                                                                      ttttAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCgtg                     <  2:3996301/93‑1 (MQ=255)
                                                                      ttttAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCgtg                     >  1:3947790/1‑93 (MQ=255)
                                                                       tttAAAAGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCgtgt                    >  1:4347936/1‑93 (MQ=255)
                                                                         tAAAAGAAACGCGCCCAGCAGACCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGtgttg                  >  1:4226670/1‑93 (MQ=255)
                                                                          aaaaGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGtgttgt                 <  1:1778233/93‑1 (MQ=255)
                                                                          aaaaGAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGtgttgt                 >  2:6237492/1‑93 (MQ=255)
                                                                              gAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATc             <  2:4967527/93‑1 (MQ=255)
                                                                              gAAACGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATc             >  1:800825/1‑93 (MQ=255)
                                                                              aaaaCGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATc             <  2:4170492/92‑1 (MQ=255)
                                                                               aaaCGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATca            >  2:5046421/1‑93 (MQ=255)
                                                                               aaaCGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATca            <  2:4391563/93‑1 (MQ=255)
                                                                                aaCGCGCCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATcac           >  2:3599295/1‑93 (MQ=255)
                                                                                   gcgcCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATcaccac        <  1:2074115/93‑1 (MQ=255)
                                                                                   gcgcCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATcaccac        >  1:3465584/1‑93 (MQ=255)
                                                                                   gcgcCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATcaccac        >  2:3339500/1‑93 (MQ=255)
                                                                                     gcCCAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCAccc      <  2:1645459/93‑1 (MQ=255)
                                                                                      cccAGCAGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCt     <  1:4063940/93‑1 (MQ=255)
                                                                                        cagcagCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCa   >  1:3272567/1‑93 (MQ=255)
                                                                                        cAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCa   <  2:4081318/93‑1 (MQ=255)
                                                                                         agcagcCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAg  <  2:1317909/93‑1 (MQ=255)
                                                                                            |                                                                                         
TTCCCGTGGCAGAGTCTGTTCTGGGCGATGGCAATGATGGGCATCGCGGTACTGATGTTGTCTTTGTTTATTTTAAAAGAAACGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAG  >  NC_000913/3889703‑3889884

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