Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,767,179 G→A *281* (TAG→TAA yibA → putative immunity protein for polymorphic toxin RhsA, HEAT‑domain protein, lethality reduction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,767,1790GA98.7% 248.2 / NA 78*281* (TAG→TAAyibAputative immunity protein for polymorphic toxin RhsA, HEAT‑domain protein, lethality reduction protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (44/33);  minor base T (0/1);  total (44/34)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCG  >  NC_000913/3767061‑3767311
                                                                                                                      |                                                                                                                                    
ttGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAg                                                                                                                  >  1:700619/1‑139 (MQ=255)
    aGCTGCGGGTGGATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCt                                                                                                                          <  2:1209099/127‑1 (MQ=255)
    aGCTGCGGGTGGATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCt                                                                                                                          >  1:1209099/1‑127 (MQ=255)
     gCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                               <  2:512114/121‑1 (MQ=255)
     gCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                               >  1:512114/1‑121 (MQ=255)
     gCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGc                                                                                                                            >  1:216669/1‑124 (MQ=255)
     gCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGc                                                                                                                            <  2:216669/124‑1 (MQ=255)
                  tAGGTGATAAAACGCTATTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGATAATGAAATTATAAATTCCGCTATTGGTAAGCGGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGTCGGc                                                                                                <  2:715356/139‑1 (MQ=255)
                   aGGTGATAAAACGGTACTTGCTGTTTTAGATACTATGTTGTACAAGTTTGATGACGATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                               <  1:1046368/139‑1 (MQ=255)
                        aTAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATg                                                                                                                       >  2:1241391/1‑110 (MQ=255)
                        aTAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATg                                                                                                                       <  1:1241391/110‑1 (MQ=255)
                         tAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGcaca                                                                                                                   >  2:140375/1‑113 (MQ=255)
                         tAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGcaca                                                                                                                   <  1:140375/113‑1 (MQ=255)
                              cGCTACTTCCTGTTTTAGATACTATGTCGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGTGATTCCCGCCGGCACGGTGCATATg                                                                                    >  1:924926/1‑139 (MQ=255)
                                    ttCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGAc                                                                              <  2:700619/139‑1 (MQ=255)
                                               gATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAg                                                                   >  1:443496/1‑139 (MQ=255)
                                               gATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAg                                                                   >  1:616507/1‑139 (MQ=255)
                                                    tATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGca                                                                                                                     <  2:1254703/84‑1 (MQ=255)
                                                    tATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGca                                                                                                                     <  2:641891/84‑1 (MQ=255)
                                                    tATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGca                                                                                                                     >  1:1254703/1‑84 (MQ=255)
                                                    tATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGca                                                                                                                     >  1:641891/1‑84 (MQ=255)
                                                    tATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGTCTGGAAAGTGaaa                                                              >  2:109634/1‑139 (MQ=255)
                                                    tATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTCAGTTGCCTATGCACAGAGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGaaa                                                              >  2:178275/1‑139 (MQ=255)
                                                    tATGTTGTACAAGTTTGAGGACAATGAAATTATAACTTCCGCTAATGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTTGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGaaa                                                              >  2:1012348/1‑139 (MQ=255)
                                                      tgttgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATATCGTGAGTTGCCTATGca                                                                                                                     <  1:896294/82‑1 (MQ=255)
                                                      tgttgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGca                                                                                                                     >  2:896294/1‑82 (MQ=255)
                                                      tgttgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGAtt                                                                                                         >  1:1023563/1‑94 (MQ=255)
                                                      tgttgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGAtt                                                                                                         <  2:1023563/94‑1 (MQ=255)
                                                         tgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTg                                                         >  1:674120/1‑139 (MQ=255)
                                                          gtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGTGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTgg                                                        >  2:1356737/1‑139 (MQ=255)
                                                           tACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGc                                                       >  1:450022/1‑139 (MQ=255)
                                                              aaGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGa                                                    >  1:983266/1‑139 (MQ=255)
                                                                gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGcc                                                                                                                           <  1:622990/66‑1 (MQ=255)
                                                                gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGcc                                                                                                                           >  2:622990/1‑66 (MQ=255)
                                                                gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGca                                                                                                                     >  2:252769/1‑72 (MQ=255)
                                                                gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGca                                                                                                                     <  1:252769/72‑1 (MQ=255)
                                                                gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCg                                                                                >  1:685179/1‑109 (MQ=255)
                                                                gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCg                                                                                <  2:685179/109‑1 (MQ=255)
                                                                gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTgg                                                        <  1:435991/133‑1 (MQ=255)
                                                                gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTgg                                                        >  2:435991/1‑133 (MQ=255)
                                                                gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGAcc                                                  >  1:1343996/1‑139 (MQ=255)
                                                                   tgatgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCggg                                               >  2:1069979/1‑139 (MQ=255)
                                                                     atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGcac                                                                                                                    >  1:707214/1‑68 (MQ=255)
                                                                     atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGcac                                                                                                                    <  2:707214/68‑1 (MQ=255)
                                                                          aaTGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgtgagtga                                        >  2:59155/1‑139 (MQ=255)
                                                                                   ataACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCgg                                                                           >  2:510277/1‑95 (MQ=255)
                                                                                   ataACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCgg                                                                           <  1:510277/95‑1 (MQ=255)
                                                                                          ccGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                               >  2:998281/1‑68 (MQ=255)
                                                                                          ccGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                               <  1:998281/68‑1 (MQ=255)
                                                                                          ccGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCgg                                                                           >  2:764163/1‑88 (MQ=255)
                                                                                          ccGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCgg                                                                           <  1:764163/88‑1 (MQ=255)
                                                                                                 tGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                               >  1:316413/1‑61 (MQ=255)
                                                                                                 tGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                               <  2:316413/61‑1 (MQ=255)
                                                                                                 tGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCAt                 <  1:178275/139‑1 (MQ=255)
                                                                                                     aaGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCgctggc                       >  2:626052/1‑129 (MQ=255)
                                                                                                     aaGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCgctggc                       <  1:626052/129‑1 (MQ=255)
                                                                                                     aaGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTg             <  2:443496/139‑1 (MQ=255)
                                                                                                          gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCt                                                      <  1:36526/93‑1 (MQ=255)
                                                                                                          gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCt                                                      >  2:36526/1‑93 (MQ=255)
                                                                                                          gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGa                                                    <  2:588564/95‑1 (MQ=255)
                                                                                                          gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGa                                                    >  1:588564/1‑95 (MQ=255)
                                                                                                          gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGccgccg                            >  1:559334/1‑119 (MQ=255)
                                                                                                          gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGccgccg                            <  2:559334/119‑1 (MQ=255)
                                                                                                          gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTg             >  2:169444/1‑134 (MQ=255)
                                                                                                          gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTg             <  1:169444/134‑1 (MQ=255)
                                                                                                          gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGcc        <  1:1162330/139‑1 (MQ=255)
                                                                                                          gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGcc        >  1:731909/1‑139 (MQ=255)
                                                                                                          gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGcc        >  2:1162330/1‑139 (MQ=255)
                                                                                                                ttCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGAcc                                                  <  2:197285/91‑1 (MQ=255)
                                                                                                                ttCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGAcc                                                  >  1:197285/1‑91 (MQ=255)
                                                                                                                ttCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGccgccg                            <  2:800062/113‑1 (MQ=255)
                                                                                                                ttCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGccgccg                            >  1:1024060/1‑113 (MQ=255)
                                                                                                                ttCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGccgccg                            <  2:1024060/113‑1 (MQ=255)
                                                                                                                ttCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGccgccg                            >  1:800062/1‑113 (MQ=255)
                                                                                                                 tCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCg  <  2:1020964/138‑1 (MQ=255)
                                                                                                                 tCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCg  >  1:1020964/1‑138 (MQ=255)
                                                                                                                  cATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCggg                                               >  1:1406950/1‑92 (MQ=255)
                                                                                                                  cATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCggg                                               <  2:1406950/92‑1 (MQ=255)
                                                                                                                      |                                                                                                                                    
TTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCG  >  NC_000913/3767061‑3767311

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