Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,781,027 C→G intergenic (‑151/+58) 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,0270CG99.7% 969.7 / ‑5.0 322intergenic (‑151/+58)ypjA/b2648putative ATP‑binding component of a transport system/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (143/178);  minor base A (0/1);  total (143/179)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

Only 100 of 327 total aligned reads displayed.

AACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGCCGTTGCCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTT  >  NC_000913/2780938‑2781116
                                                                                          |                                                                                         
aaCGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCa                                                                                         >  2:3619667/1‑93 (MQ=255)
  cGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaa                                                                                       <  2:69151/93‑1 (MQ=255)
   gCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaa                                                                                      <  2:3754935/93‑1 (MQ=25)
   gCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaa                                                                                      >  2:6278534/1‑93 (MQ=255)
    cAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaat                                                                                     <  2:5706666/93‑1 (MQ=25)
     aGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaata                                                                                    >  1:3269018/1‑93 (MQ=255)
        gcgcTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaataaaa                                                                                 <  1:1281499/93‑1 (MQ=25)
           cTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGa                                                                              >  1:6922595/1‑93 (MQ=255)
           cTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGa                                                                              >  2:4989161/1‑93 (MQ=255)
            tGAAGTGTGGTTGGATAAAAAAGTCAATCCATTTAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAt                                                                             <  2:2018226/93‑1 (MQ=25)
            tGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAt                                                                             <  1:7759436/93‑1 (MQ=255)
            tGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAt                                                                             >  2:4179156/1‑93 (MQ=255)
             gAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAtt                                                                            <  1:4620255/93‑1 (MQ=37)
             gAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAtt                                                                            <  2:7266368/93‑1 (MQ=25)
              aaGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAttt                                                                           <  1:7432124/93‑1 (MQ=37)
                ttgtgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTGTCGACAATGAGGCCGTGGGCAAAATAAAATGATTTAc                                                                         <  2:4674792/92‑1 (MQ=25)
                 tgtgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACa                                                                        >  1:7084582/1‑93 (MQ=255)
                 tgtgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACa                                                                        >  2:3830/1‑93 (MQ=255)
                  gtgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCCACAATGAGGCCGTTGGCAAAATAAAATGATTTACAt                                                                       >  1:4182755/1‑93 (MQ=37)
                   tgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATa                                                                      <  2:4020390/93‑1 (MQ=37)
                    ggTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATaa                                                                     <  1:1053762/93‑1 (MQ=255)
                     gTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAAt                                                                    >  1:3350611/1‑93 (MQ=255)
                     gTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAAt                                                                    <  2:3468437/93‑1 (MQ=37)
                      ttGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATc                                                                   <  2:6625686/93‑1 (MQ=37)
                      ttGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATc                                                                   >  1:2438030/1‑93 (MQ=255)
                      ttGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATc                                                                   <  2:7330617/93‑1 (MQ=255)
                        ggATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGt                                                                 >  1:3824781/1‑93 (MQ=255)
                        ggATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGt                                                                 >  1:7858584/1‑93 (MQ=255)
                         gATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGtt                                                                >  1:4068553/1‑93 (MQ=255)
                          aTAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGttt                                                               <  2:1886756/93‑1 (MQ=37)
                          aTAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGttt                                                               <  2:3293635/93‑1 (MQ=38)
                            aaaaaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCt                                                             >  2:5486591/1‑93 (MQ=255)
                             aaaaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtg                                                            <  1:7280491/93‑1 (MQ=37)
                             aaaaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtg                                                            >  1:3888597/1‑93 (MQ=255)
                              aaaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtga                                                           >  1:4813164/1‑93 (MQ=255)
                               aaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtgat                                                          >  2:1049639/1‑93 (MQ=255)
                               aaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtgat                                                          <  2:7199913/93‑1 (MQ=255)
                               aaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtgat                                                          >  1:5414265/1‑93 (MQ=255)
                                 aGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtgatga                                                        <  2:4453316/93‑1 (MQ=255)
                                   tCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAat                                                      <  2:1447402/93‑1 (MQ=255)
                                   tCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAat                                                      <  1:3041130/93‑1 (MQ=255)
                                     aaTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAatat                                                    <  1:1893623/93‑1 (MQ=255)
                                     aaTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAatat                                                    >  2:4943568/1‑93 (MQ=255)
                                     aaTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAatat                                                    >  1:7576274/1‑93 (MQ=255)
                                     aaTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAatat                                                    <  2:119304/93‑1 (MQ=25)
                                     aaTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAatat                                                    >  2:608295/1‑93 (MQ=255)
                                      aTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtc                                                   <  2:7836584/93‑1 (MQ=25)
                                      aTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtc                                                   <  2:6023839/93‑1 (MQ=255)
                                      aTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtc                                                   >  2:1322782/1‑93 (MQ=255)
                                      aTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtc                                                   >  1:1627321/1‑93 (MQ=255)
                                      aTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtc                                                   >  1:5127805/1‑93 (MQ=255)
                                        ccATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtctt                                                 <  2:510762/93‑1 (MQ=255)
                                        ccATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtctt                                                 >  1:5333016/1‑93 (MQ=255)
                                         cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                                >  2:6623004/1‑93 (MQ=255)
                                         cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                                <  1:2026110/93‑1 (MQ=255)
                                         cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                                <  1:6328534/93‑1 (MQ=25)
                                          aTTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttct                                               <  1:6626165/93‑1 (MQ=255)
                                          aTTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttct                                               >  2:7777670/1‑93 (MQ=255)
                                          aTTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttct                                               <  2:7668421/93‑1 (MQ=255)
                                           ttCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTg                                              >  1:3641061/1‑93 (MQ=255)
                                            tCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGc                                             >  1:3422711/1‑93 (MQ=255)
                                               ggAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTca                                          <  2:5592364/93‑1 (MQ=255)
                                               ggAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTca                                          >  2:2776868/1‑93 (MQ=255)
                                                  aaTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACAt                                       <  2:6122889/93‑1 (MQ=255)
                                                  aaTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACAt                                       <  2:5547326/93‑1 (MQ=255)
                                                    tACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaa                                     >  1:5854120/1‑93 (MQ=18)
                                                     aCGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaaa                                    >  1:25768/1‑93 (MQ=18)
                                                        ggCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAAt                                 >  2:584501/1‑93 (MQ=255)
                                                        ggCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAAt                                 <  1:5694026/93‑1 (MQ=255)
                                                         gCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATc                                <  1:4303917/93‑1 (MQ=255)
                                                          cGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATca                               <  2:2472806/93‑1 (MQ=255)
                                                           gTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcac                              <  2:4524104/93‑1 (MQ=255)
                                                            tATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcaca                             >  1:6006246/1‑93 (MQ=17)
                                                            tATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcaca                             >  2:6526410/1‑93 (MQ=18)
                                                            tATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcaca                             <  2:290256/93‑1 (MQ=17)
                                                            tATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcaca                             >  2:5387090/1‑93 (MQ=255)
                                                             aTTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcacac                            <  1:1714796/93‑1 (MQ=255)
                                                              ttcttttcttTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcacaca                           >  1:6432116/1‑93 (MQ=17)
                                                              ttcttttcttTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcacaca                           >  1:6793765/1‑93 (MQ=17)
                                                               tcttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACaa                          <  2:5788284/93‑1 (MQ=255)
                                                                 ttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATa                        <  1:6548252/93‑1 (MQ=255)
                                                                 ttttctttCGACAATGAGGCCGCTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATa                        >  1:2692923/1‑93 (MQ=14)
                                                                   ttctgtcGACAATGAGGCCCTTGGCAAAATAAAATGATGTACATAATCGTTTCTGAGGAATATCTTCTGCTCACATAAAAATCACACAATAAc                      <  1:950005/93‑1 (MQ=14)
                                                                    tctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACt                     >  1:880559/1‑93 (MQ=17)
                                                                    tctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACt                     >  1:4006010/1‑93 (MQ=17)
                                                                    tctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACt                     <  1:4922650/93‑1 (MQ=255)
                                                                     ctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACtt                    <  2:4818363/93‑1 (MQ=255)
                                                                      tttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACttt                   <  2:399121/93‑1 (MQ=255)
                                                                         cGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTgag                >  2:3659227/1‑93 (MQ=18)
                                                                          gACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTgaga               <  1:3498722/93‑1 (MQ=255)
                                                                           aCAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGAt              <  2:2540364/93‑1 (MQ=255)
                                                                           aCAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGAt              <  2:2177708/93‑1 (MQ=255)
                                                                             aaTGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCg            >  2:6655445/1‑93 (MQ=18)
                                                                               tgaggCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCa          >  1:4582771/5‑93 (MQ=255)
                                                                               tgaggCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCa          <  2:2577984/89‑1 (MQ=255)
                                                                                   gCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCCCATAAAAATCACACAATAACTTTGAGACCGCAGAtt      >  1:5994788/1‑93 (MQ=21)
                                                                                      gTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGtt   <  2:2393158/93‑1 (MQ=255)
                                                                                      gTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGtt   >  2:172152/1‑93 (MQ=255)
                                                                                       ttGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGttt  >  1:4232412/1‑93 (MQ=255)
                                                                                       ttGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGttt  >  1:2261915/1‑93 (MQ=255)
                                                                                          |                                                                                         
AACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGCCGTTGCCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTT  >  NC_000913/2780938‑2781116

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