Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,781,020 +G intergenic (‑144/+65) ypjA ← / ← b2648 putative ATP‑binding component of a transport system/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,781,0191.G99.6% 756.5 / NA 226intergenic (‑143/+66)ypjA/b2648putative ATP‑binding component of a transport system/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base G (96/129);  minor base T (0/1);  total (96/130)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 250 total aligned reads displayed.

GGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGCCGTTGCCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCG  >  NC_000913/2780927‑2781106
                                                                                             |                                                                                       
ggCATTTTTTAAACGCAGTGCGCTGAGATGTGATAGGATAGAAACGTAAAAACAAAAAGAACATACGGGCGTATTCTTTTCTTTCGACAATGa                                                                                          <  2:5467396/93‑1 (MQ=255)
ggCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGa                                                                                          >  1:3463370/1‑93 (MQ=255)
ggCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGa                                                                                          <  1:1778872/93‑1 (MQ=255)
ggCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGa                                                                                          >  2:4949148/1‑93 (MQ=255)
ggCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGa                                                                                          <  2:3200895/93‑1 (MQ=255)
 gCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAg                                                                                        >  1:71682/1‑93 (MQ=255)
   aTTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAggc                                                                                      <  2:3309237/93‑3 (MQ=255)
    ttttttAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAggcc                                                                                     <  1:154918/93‑4 (MQ=255)
    ttttttAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAggcc                                                                                     >  2:5236950/1‑90 (MQ=255)
    ttttttAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAggcc                                                                                     >  2:3280648/1‑90 (MQ=255)
       tttAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGtt                                                                                   <  2:6060121/93‑1 (MQ=255)
        ttAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTg                                                                                  >  2:2897509/1‑93 (MQ=255)
         tAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTgg                                                                                 <  1:502854/93‑2 (MQ=255)
         tAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTgg                                                                                 <  1:4756582/93‑2 (MQ=255)
         tAAACGCAGTGCGCTGAAGTGTGATTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGGCAATGAGGCCGTTgg                                                                                 <  1:1494261/93‑2 (MQ=25)
          aaaCGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTggc                                                                                >  2:3867005/1‑91 (MQ=255)
           aaCGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCa                                                                               <  1:4453412/93‑1 (MQ=255)
            aCGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaa                                                                              <  2:4115255/93‑1 (MQ=255)
              gCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaa                                                                            >  1:668714/1‑93 (MQ=255)
                aGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaata                                                                          >  1:6103949/1‑93 (MQ=255)
                 gTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaataa                                                                         >  2:5987642/1‑93 (MQ=255)
                      cTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGa                                                                    <  2:361806/93‑1 (MQ=255)
                      cTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGa                                                                    <  2:6020269/93‑1 (MQ=25)
                       tGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAt                                                                   <  1:892280/93‑1 (MQ=255)
                        gAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGCAATACGGGCGTATTCTTTTCTTTCGACAATAAGGCCGTTGGCAAAATACAATGagt                                                                  >  2:4704759/1‑91 (MQ=255)
                        gAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAtt                                                                  <  1:2802324/93‑1 (MQ=25)
                        gAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAtt                                                                  >  1:4039387/1‑93 (MQ=255)
                        gAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAtt                                                                  >  2:1538137/1‑93 (MQ=255)
                        gAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAtt                                                                  >  2:180211/1‑93 (MQ=255)
                          aGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTa                                                                >  2:3528869/1‑93 (MQ=255)
                           aggtgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTTGCAAAATAAAATGATTTAc                                                               <  2:2182889/91‑1 (MQ=25)
                            tgtgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACa                                                              <  2:706849/93‑1 (MQ=255)
                              tgGTTGGATAACAAAGTCCATCCATTCAGGCAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCACAATCAAATGATTTACATa                                                            >  2:12820/1‑93 (MQ=25)
                              tgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATa                                                            >  1:1208765/1‑93 (MQ=255)
                              tgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATa                                                            >  1:2104542/1‑93 (MQ=255)
                                gTTGTATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAAt                                                          <  1:3241719/93‑1 (MQ=37)
                                gTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAAt                                                          <  1:5191718/93‑1 (MQ=37)
                                gTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAAt                                                          >  2:4553265/1‑93 (MQ=255)
                                gTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAAt                                                          <  2:4435044/93‑1 (MQ=255)
                                 ttGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATc                                                         <  2:4726526/93‑1 (MQ=255)
                                 ttGGATAAAAAAGGCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATc                                                         <  1:5272021/93‑1 (MQ=37)
                                   ggATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGt                                                       <  2:1457205/93‑1 (MQ=37)
                                    gATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGtt                                                      <  2:3935276/93‑1 (MQ=37)
                                    gATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGtt                                                      <  1:2212920/93‑1 (MQ=37)
                                      tAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTc                                                    >  2:1607512/1‑93 (MQ=255)
                                        aaaaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtg                                                  >  1:1900860/1‑93 (MQ=255)
                                          aaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtgat                                                <  1:5721860/93‑1 (MQ=37)
                                          aaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtgat                                                <  1:4906888/93‑1 (MQ=255)
                                          aaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtgat                                                <  1:1471301/93‑1 (MQ=37)
                                            aGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtgatga                                              <  1:5835389/93‑1 (MQ=255)
                                            aGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtgatga                                              <  2:1234956/93‑1 (MQ=37)
                                             gTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGaa                                             <  1:1085952/93‑1 (MQ=255)
                                             gTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGaa                                             >  2:5406552/1‑93 (MQ=255)
                                              tCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAat                                            <  1:5249976/93‑1 (MQ=37)
                                              tCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAat                                            <  1:4740006/93‑1 (MQ=37)
                                              tCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAat                                            <  1:4470870/93‑1 (MQ=37)
                                                aaTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAatat                                          >  2:4355561/1‑93 (MQ=255)
                                                 aTCCATTCAGGAAATACGGGCGTATTCTTTTCTTGCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtc                                         >  1:3975883/1‑93 (MQ=25)
                                                  tCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtct                                        >  2:2658864/1‑93 (MQ=255)
                                                   ccATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtctt                                       >  1:4613998/1‑93 (MQ=255)
                                                   ccATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtctt                                       <  2:5769366/93‑1 (MQ=25)
                                                    cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                      <  2:4200763/93‑1 (MQ=255)
                                                    cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                      <  2:4923993/93‑1 (MQ=255)
                                                    cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                      <  1:1709594/93‑1 (MQ=25)
                                                    cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                      <  1:3422250/93‑1 (MQ=255)
                                                    cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                      >  2:244320/1‑93 (MQ=255)
                                                     aTTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttct                                     >  1:3353707/1‑93 (MQ=255)
                                                     aTTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttct                                     <  1:1772424/93‑1 (MQ=255)
                                                     aTTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttct                                     >  2:3796052/1‑93 (MQ=255)
                                                      ttCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTg                                    >  1:3691999/1‑93 (MQ=255)
                                                      ttCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTg                                    >  2:2045401/1‑93 (MQ=255)
                                                          ggAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTca                                >  2:64491/1‑93 (MQ=255)
                                                          ggAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAGATACTTAATCGTTTCTGATGAATAGCTTCTGCTct                                >  1:5493282/1‑92 (MQ=18)
                                                            aaaTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTcaca                              <  1:379178/93‑1 (MQ=21)
                                                            aaaTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTGGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTcaca                              <  1:537472/93‑1 (MQ=18)
                                                             aaTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACAt                             <  2:2878388/93‑1 (MQ=18)
                                                                  gggCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaaaaa                        <  2:1770221/93‑1 (MQ=18)
                                                                  cggCGTATTCTTTTCTTTCGACAATGATGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaaaaa                        <  1:5340472/92‑1 (MQ=18)
                                                                   ggCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAAt                       >  2:412648/1‑93 (MQ=18)
                                                                     cGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGCTTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATca                     >  2:4283885/1‑93 (MQ=17)
                                                                       tATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcaca                   >  1:4495732/1‑93 (MQ=18)
                                                                        aTTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcacac                  <  2:856990/93‑1 (MQ=255)
                                                                        aTTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcacac                  >  2:4501290/1‑93 (MQ=17)
                                                                         ttcttttcttTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcacaca                 <  1:1394915/93‑1 (MQ=17)
                                                                          tcttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATACCACaa                >  1:1917698/1‑93 (MQ=14)
                                                                            ttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATa              <  2:4598862/93‑1 (MQ=255)
                                                                            ttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATa              <  1:1640365/93‑1 (MQ=255)
                                                                            ttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATa              <  1:2668862/93‑1 (MQ=255)
                                                                             tttctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATaa             <  1:5138188/93‑1 (MQ=255)
                                                                             tttctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATaa             >  1:978732/1‑93 (MQ=17)
                                                                             tttctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATaa             <  1:2478549/93‑1 (MQ=255)
                                                                                 tttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAAAttt         >  2:1087664/1‑93 (MQ=14)
                                                                                  ttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTg        <  2:2624865/93‑1 (MQ=255)
                                                                                   tcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTga       <  2:848897/93‑1 (MQ=255)
                                                                                   tcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTga       <  2:437104/93‑1 (MQ=255)
                                                                                       cAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATc   <  2:5032083/93‑1 (MQ=255)
                                                                                        aaTGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCg  >  2:3055117/1‑93 (MQ=18)
                                                                                        aaTGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCg  >  2:2688023/1‑93 (MQ=18)
                                                                                        aaTGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCg  <  1:2356272/93‑1 (MQ=255)
                                                                                        aaTGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCg  >  2:4067431/1‑93 (MQ=17)
                                                                                             |                                                                                       
GGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGCCGTTGCCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCG  >  NC_000913/2780927‑2781106

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