Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,586,436 A→G A17A (GCT→GCC ydeP ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,586,4360AG98.6% 258.3 / NA 73A17A (GCT→GCCydePputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (34/38);  minor base T (0/1);  total (34/39)
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

CGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGCACAGTACATACAAAGGGATTATCGCCAACATGG  >  NC_000913/1586304‑1586571
                                                                                                                                    |                                                                                                                                       
cGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCAccc                                                                                                                                   <  2:740324/139‑1 (MQ=255)
       tgCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAaccac                                                                                                                            >  2:577541/1‑139 (MQ=255)
          ttAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTg                                                                                                                         <  1:355926/139‑1 (MQ=255)
            aGGATCTGGCCATGCACAACCCGGACAGTCAAAGTCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTgca                                                                                                                       <  1:600430/139‑1 (MQ=255)
               aTCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTGCTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTgcagca                                                                                                                    <  2:544692/139‑1 (MQ=255)
                        tGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAg                                                                                                           >  1:603838/1‑139 (MQ=255)
                              aCCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCa                                                                                                     >  2:152713/1‑139 (MQ=255)
                                  ggACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAAttt                                                                                                 >  1:311217/1‑139 (MQ=255)
                                   gACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAAtttt                                                                                                >  2:63206/1‑139 (MQ=255)
                                       gTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTttc                                                                                            >  1:643868/1‑139 (MQ=255)
                                       gTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTttc                                                                                            >  1:216737/1‑139 (MQ=255)
                                           aaGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCa                                                                                        <  1:317700/139‑1 (MQ=255)
                                           aaGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCa                                                                                        <  2:603838/139‑1 (MQ=255)
                                           aaGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCa                                                                                        <  2:598665/139‑1 (MQ=255)
                                           aaGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACAGCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCa                                                                                        <  1:627455/139‑1 (MQ=255)
                                           aaGCCCTCTCGCTTATTCATGTCAAACATGCCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCCCTACGGATTTAACTGCCCCCCAACCTCCTGCAGCACTCTGGTAGGATTCAATTTTTTTCTTCa                                                                                        <  2:115016/139‑1 (MQ=255)
                                                ctctGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTAc                                                                                   >  1:5803/1‑139 (MQ=255)
                                                  ctGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCa                                                                                        <  1:401605/132‑1 (MQ=255)
                                                  ctGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCa                                                                                        >  2:401605/1‑132 (MQ=255)
                                                  ctGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACtt                                                                                 >  1:605896/1‑139 (MQ=255)
                                                     gCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATc                                                                              <  2:582184/139‑1 (MQ=255)
                                                         aTTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCtgt                                                                          >  1:722824/1‑139 (MQ=255)
                                                          ttCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCa                                                                                                                                 <  2:314393/83‑1 (MQ=255)
                                                          ttCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCa                                                                                                                                 >  1:314393/1‑83 (MQ=255)
                                                               gTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTgca                                                                                                                       <  1:566090/88‑1 (MQ=255)
                                                               gTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTgca                                                                                                                       >  2:566090/1‑88 (MQ=255)
                                                               gTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGaa                                                                    >  1:334032/1‑139 (MQ=255)
                                                               gTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGaa                                                                    >  2:50533/1‑139 (MQ=255)
                                                                         gCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAattatt                                                          >  2:304423/1‑139 (MQ=255)
                                                                         gCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAattatt                                                          >  1:727290/1‑139 (MQ=255)
                                                                          cAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAAttatta                                                         <  1:152713/139‑1 (MQ=255)
                                                                             tAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGa                                                      >  1:650896/1‑139 (MQ=255)
                                                                                  tCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGt                                                 >  2:109352/1‑139 (MQ=255)
                                                                                   ccTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTg                                                                                                                         <  1:508341/66‑1 (MQ=255)
                                                                                   ccTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTAAGGATTTAACGGCACTCCAACCACCTg                                                                                                                         >  2:508341/1‑66 (MQ=255)
                                                                                     tGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAAt                                              <  1:433595/139‑1 (MQ=255)
                                                                                       gCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAg                                            <  2:216737/139‑1 (MQ=255)
                                                                                        cGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGc                                           >  2:603894/1‑139 (MQ=255)
                                                                                           atatCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCaa                                                                                                                                >  1:119007/1‑51 (MQ=255)
                                                                                           atatCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCaa                                                                                                                                <  2:119007/51‑1 (MQ=255)
                                                                                            tatCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCgatg                                       >  1:207624/1‑139 (MQ=255)
                                                                                                cATCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGc                                   <  1:577541/139‑1 (MQ=255)
                                                                                                 aTCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAgg                                                       <  2:449522/118‑1 (MQ=255)
                                                                                                 aTCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAgg                                                       >  1:449522/1‑118 (MQ=255)
                                                                                                 aTCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGca                                  >  1:9536/1‑139 (MQ=255)
                                                                                                 aTCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGca                                  <  1:109352/139‑1 (MQ=255)
                                                                                                 aTCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGca                                  <  2:5803/139‑1 (MQ=255)
                                                                                                 aTCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGca                                  <  2:600988/139‑1 (MQ=255)
                                                                                                 aTCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGca                                  <  2:311217/139‑1 (MQ=255)
                                                                                                 aTCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGca                                  <  1:63206/139‑1 (MQ=255)
                                                                                                  tCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGcac                                 <  2:605896/139‑1 (MQ=255)
                                                                                                  tCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGcac                                 <  2:643868/139‑1 (MQ=255)
                                                                                                  tCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGcac                                 <  2:650896/139‑1 (MQ=255)
                                                                                                  tCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGcac                                 <  1:285152/139‑1 (MQ=255)
                                                                                                  tCTGCTTACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGcac                                 <  1:50533/139‑1 (MQ=255)
                                                                                                       ttACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGCACAGtac                            >  1:630995/1‑139 (MQ=255)
                                                                                                        tACGTACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGCACAGtaca                           >  2:745776/1‑139 (MQ=255)
                                                                                                            tACTGCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGCACAGtacataca                       <  1:12587/139‑1 (MQ=255)
                                                                                                                gCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCtgtg                                                                         >  2:406610/1‑85 (MQ=255)
                                                                                                                gCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCtgtg                                                                         <  1:406610/85‑1 (MQ=255)
                                                                                                                gCATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGCACAGTACATACAAAgg                   >  2:312427/1‑139 (MQ=255)
                                                                                                                 cATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGcac                                 <  2:375839/124‑1 (MQ=255)
                                                                                                                 cATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGcac                                 >  1:375839/1‑124 (MQ=255)
                                                                                                                 cATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGca                                  >  2:110999/1‑123 (MQ=255)
                                                                                                                 cATTCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGca                                  <  1:110999/123‑1 (MQ=255)
                                                                                                                   ttCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCGTGTGTGGAAATGAATTATTAGGACTTGTAATAGCAATGATGCACAGTACATACAAAGGGAt                >  2:491569/1‑139 (MQ=255)
                                                                                                                   ttCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAAttatta                                                         >  2:430550/1‑98 (MQ=255)
                                                                                                                   ttCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAAttatta                                                         <  1:430550/98‑1 (MQ=255)
                                                                                                                    tCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTgg                                                                                                              <  2:405321/44‑1 (MQ=255)
                                                                                                                    tCGCTACGGATTTAACGGCACCCCAACCACCTGCAGCACCCTgg                                                                                                              >  1:405321/1‑44 (MQ=255)
                                                                                                                                 aaCGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGCACAGTACATACAAAGGGATTATCGCCAACATgg  <  2:334032/139‑1 (MQ=255)
                                                                                                                                 aaCGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGCACAGTACATACAAAGGGATTATCGCCAACATgg  <  1:304423/139‑1 (MQ=255)
                                                                                                                                 aaCGGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACCTAACCTTTTTTGAAATGAATTAATAGGATTTGTAATAGCGATGATGGACAGTACATAAAAAGGGATTATCGCCAACATgg  <  2:722824/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                       
CGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGCACAGTACATACAAAGGGATTATCGCCAACATGG  >  NC_000913/1586304‑1586571

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