Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,787,434 G→A pseudogene (1038/2253 nt) ygaQ → uncharacterized protein, putative enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,787,4340GA98.8% 262.6 / NA 81pseudogene (1038/2253 nt)ygaQuncharacterized protein, putative enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (40/40);  minor base C (1/0);  total (41/40)
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

CATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAATGGGAA  >  NC_000913/2787308‑2787569
                                                                                                                              |                                                                                                                                       
cATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGtt                                                                                                                             <  2:1223147/139‑1 (MQ=255)
      ctaTTATAACTTTTAATTTTTGGGAATTGGGGGGAAGTGTAAAAAAGACCTTGATTTACCCCTTGGTGGTATTCGGGCGAAGATCAATAAAGGAATAGGAGGTTGTAATCCCGGCTTATAACAAGGTCAGGTTTTGAtt                                                                                                                       <  2:407902/137‑1 (MQ=255)
       gATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTc                                                                                                                      >  2:1272375/1‑139 (MQ=255)
            taCCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTACCAGGTTAGGGTTTTGATTCaatta                                                                                                                 >  1:407902/1‑138 (MQ=255)
                   aaCTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAg                                                                                                          >  1:1087581/1‑139 (MQ=255)
                       ttATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                      <  1:1272375/139‑1 (MQ=255)
                          tGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGt                                                                                                                                   >  2:400569/1‑107 (MQ=255)
                          tGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGt                                                                                                                                   <  1:400569/107‑1 (MQ=255)
                               aTTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTa                                                                                                             >  1:776068/1‑124 (MQ=255)
                               aTTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTa                                                                                                             <  2:776068/124‑1 (MQ=255)
                                  gggggTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATATGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACATAAg                                                                                           <  2:1087581/139‑1 (MQ=255)
                                        aGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAg                                                                                                    <  2:407029/124‑1 (MQ=255)
                                        aGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAg                                                                                                    >  1:407029/1‑124 (MQ=255)
                                                 aaTGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcattcat                                                                                                               <  1:82066/104‑1 (MQ=255)
                                                 aaTGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcattcat                                                                                                               >  2:82066/1‑104 (MQ=255)
                                                      ccTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcat                                                                                                                   >  1:1653042/1‑95 (MQ=255)
                                                      ccTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcat                                                                                                                   <  2:1653042/95‑1 (MQ=255)
                                                      ccTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                        <  2:1660955/138‑1 (MQ=255)
                                                      ccTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                        >  1:1660955/1‑138 (MQ=255)
                                                      ccTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAAtccatcca                                                                       >  1:203394/1‑139 (MQ=255)
                                                      ccTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAAtccatcca                                                                       >  1:502277/1‑139 (MQ=255)
                                                          gATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCagaag                                                                   <  2:203394/139‑1 (MQ=255)
                                                                     tgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAat                                                                              >  1:638111/1‑117 (MQ=255)
                                                                     tgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAat                                                                              <  2:638111/117‑1 (MQ=255)
                                                                     tgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACAtctc                                                        >  1:524512/1‑139 (MQ=255)
                                                                          atatACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCagaag                                                                   <  2:1256510/123‑1 (MQ=255)
                                                                          atatACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCagaag                                                                   >  1:1256510/1‑123 (MQ=255)
                                                                            ataCGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCaa                                                 >  1:585455/1‑139 (MQ=255)
                                                                               cGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAg                                              >  1:1088391/1‑139 (MQ=255)
                                                                                 aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAAc                                                                                        >  2:264660/1‑95 (MQ=255)
                                                                                 aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAAc                                                                                        <  1:264660/95‑1 (MQ=255)
                                                                                  gCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACt                                           >  2:142275/1‑139 (MQ=255)
                                                                                     aaGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTaa                                                                                                            >  2:1639681/1‑71 (MQ=255)
                                                                                     aaGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTaa                                                                                                            <  1:1639681/71‑1 (MQ=255)
                                                                                       gATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGaa                                      >  1:258385/1‑139 (MQ=255)
                                                                                           aGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCaa                                  <  2:585455/139‑1 (MQ=255)
                                                                                              aaaGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTa                               <  2:258385/139‑1 (MQ=255)
                                                                                                aGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCca                                                                                                                    <  1:1311240/52‑1 (MQ=255)
                                                                                                aGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCca                                                                                                                    >  2:1311240/1‑52 (MQ=255)
                                                                                                aGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAAtccatcca                                                                       >  1:1520331/1‑97 (MQ=255)
                                                                                                aGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAAtccatcca                                                                       <  2:1520331/97‑1 (MQ=255)
                                                                                                  atatAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGAtt                                                                                                                       <  1:824758/47‑1 (MQ=255)
                                                                                                  atatAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGAtt                                                                                                                       >  2:824758/1‑47 (MQ=255)
                                                                                                    ataGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGttt                                                                                                                            >  2:645761/1‑40 (MQ=39)
                                                                                                    ataGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGttt                                                                                                                            <  1:645761/40‑1 (MQ=255)
                                                                                                     taGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTaaa                             >  2:622124/1‑134 (MQ=255)
                                                                                                     taGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTaaa                             <  1:622124/134‑1 (MQ=255)
                                                                                                      aGAAGGTGGTATCCCTGTCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                       <  2:524512/139‑1 (MQ=255)
                                                                                                      aGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTTACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGAAAAc                       <  2:1088391/139‑1 (MQ=255)
                                                                                                      aGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                       <  2:502277/139‑1 (MQ=255)
                                                                                                       gAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGAtt                                                                                                                       <  1:1678527/42‑1 (MQ=255)
                                                                                                       gAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGAtt                                                                                                                       >  2:1678527/1‑42 (MQ=255)
                                                                                                        aaGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTa                                                                                                  >  1:249/1‑62 (MQ=255)
                                                                                                        aaGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTa                                                                                                  <  2:249/62‑1 (MQ=255)
                                                                                                            tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGt                                                                                                   >  2:1329986/1‑57 (MQ=255)
                                                                                                            tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGt                                                                                                   <  1:1329986/57‑1 (MQ=255)
                                                                                                            tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCaa                                  <  1:707925/122‑1 (MQ=255)
                                                                                                            tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCaa                                  >  2:707925/1‑122 (MQ=255)
                                                                                                             ggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTa                               >  2:525283/1‑124 (MQ=255)
                                                                                                             ggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTa                               <  1:525283/124‑1 (MQ=255)
                                                                                                               tATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCca                                                                                                                    >  2:628506/1‑37 (MQ=255)
                                                                                                               tATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCca                                                                                                                    <  1:628506/37‑1 (MQ=255)
                                                                                                                 tCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                                     >  1:422149/1‑82 (MQ=255)
                                                                                                                 tCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                                     <  2:422149/82‑1 (MQ=255)
                                                                                                                 tCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                       >  1:1454286/1‑128 (MQ=255)
                                                                                                                 tCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                       <  2:1454286/128‑1 (MQ=255)
                                                                                                                  cccTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGAt           >  2:223578/1‑139 (MQ=255)
                                                                                                                   ccTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCagaag                                                                   <  1:1624292/82‑1 (MQ=255)
                                                                                                                   ccTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCagaag                                                                   >  2:1624292/1‑82 (MQ=255)
                                                                                                                      ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAAtccatcca                                                                       <  2:1026386/75‑1 (MQ=255)
                                                                                                                      ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAAtccatcca                                                                       >  1:1026386/1‑75 (MQ=255)
                                                                                                                      ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCagaag                                                                   <  1:1067480/79‑1 (MQ=255)
                                                                                                                      ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCagaag                                                                   >  2:1067480/1‑79 (MQ=255)
                                                                                                                      ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGc                                    <  2:688853/110‑1 (MQ=255)
                                                                                                                      ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGc                                    >  1:688853/1‑110 (MQ=255)
                                                                                                                      ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTGAGGATGAAt       >  1:1081879/1‑139 (MQ=255)
                                                                                                                      ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAAt       >  1:208560/1‑139 (MQ=255)
                                                                                                                       gCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGat                                                                 <  2:250989/80‑1 (MQ=255)
                                                                                                                       gCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGat                                                                 >  1:250989/1‑80 (MQ=255)
                                                                                                                           ttAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACaaa                                                                                             <  2:96275/48‑1 (MQ=255)
                                                                                                                           ttAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAATGGGaa  <  2:1081879/139‑1 (MQ=255)
                                                                                                                              |                                                                                                                                       
CATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAATGGGAA  >  NC_000913/2787308‑2787569

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