Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,736,916 A→G 7.7% intergenic (+33/‑238) bamD → / → raiA BamABCDE complex OM biogenesis lipoprotein/cold shock protein associated with 30S ribosomal subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,736,9160AG 7.7% 225.3 / 7.3 78intergenic (+33/‑238)bamD/raiABamABCDE complex OM biogenesis lipoprotein/cold shock protein associated with 30S ribosomal subunit
Reads supporting (aligned to +/- strand):  ref base A (36/36);  new base G (2/4);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

CCAGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTGA  >  NC_000913/2736769‑2737057
                                                                                                                                                   |                                                                                                                                             
ccAGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAg                                                                                                                                              >  1:26662/1‑149 (MQ=255)
 cAGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGc                                                                                                                                             >  6:197646/1‑149 (MQ=255)
      tACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCg                                                                                                                                        >  4:142137/1‑149 (MQ=255)
      tACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCg                                                                                                                                        >  5:322133/1‑149 (MQ=255)
          cgTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCTCAAGGGCTGCCGtttt                                                                                                                                    <  7:142135/149‑1 (MQ=255)
           gTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATTATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGttttt                                                                                                                                   <  6:83196/149‑1 (MQ=255)
              aTGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGttttttt                                                                                                                                 <  8:7603/148‑1 (MQ=255)
                 cgcTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTc                                                                                                                             <  5:132056/149‑1 (MQ=255)
                 cgcTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTc                                                                                                                             <  6:348823/149‑1 (MQ=255)
                        gCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGt                                                                                                                      >  2:271292/1‑149 (MQ=255)
                         cTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGtt                                                                                                                     <  5:361572/149‑1 (MQ=255)
                          tGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCTCAAGGGCTGCCGtttt                                                                                                                                    >  5:202073/1‑133 (MQ=255)
                          tGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCTCAAGGGCTGCCGtttt                                                                                                                                    >  5:207525/1‑133 (MQ=255)
                          tGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCTCAAGGGCTGCCGtttt                                                                                                                                    <  6:207525/133‑1 (MQ=255)
                          tGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCTCAAGGGCTGCCGtttt                                                                                                                                    <  6:202073/133‑1 (MQ=255)
                              ggAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTAtt                                                                                                                              >  2:18914/1‑135 (MQ=255)
                              ggAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTAtt                                                                                                                              <  1:18914/135‑1 (MQ=255)
                                    tGCATACCGTCAGATTCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGaa                                                                                                          <  1:215320/149‑1 (MQ=255)
                                    tGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGt                                                                                                                          <  7:27825/133‑1 (MQ=255)
                                    tGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGt                                                                                                                          >  8:27825/1‑133 (MQ=255)
                                    tGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGaa                                                                                                          <  6:295562/149‑1 (MQ=255)
                                       aTACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGc                                                                                                        >  7:193952/1‑148 (MQ=255)
                                          ccGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGAtt                                                                                                    >  6:262831/1‑149 (MQ=255)
                                          ccGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATAGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGAt                                                                                                     <  6:322133/148‑1 (MQ=255)
                                           cGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGAtt                                                                                                    >  6:257169/1‑148 (MQ=255)
                                            gTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTa                                                                                                  <  3:142137/149‑1 (MQ=255)
                                             tCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAg                                                                                                 <  8:9746/149‑1 (MQ=255)
                                                    cagatgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTa                                                                                           <  1:248879/148‑1 (MQ=255)
                                                      gatgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGaa                                                                                                          <  5:91061/131‑1 (MQ=255)
                                                      gatgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGaa                                                                                                          >  6:91061/1‑131 (MQ=255)
                                                            tGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCtca                                                                                  >  6:18576/1‑149 (MQ=255)
                                                              cgcAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCtcatc                                                                                <  8:124697/149‑1 (MQ=255)
                                                                  aGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAACCTGAAGCGATTTAGTCGCTATCGATCTCATCAAat                                                                            <  2:26662/149‑1 (MQ=255)
                                                                     tGAAAAAGTATCGAAAATCATCGCCGCAAACAGCATCAATACATAACAGAAACATGAAAAACAAAACGGCATCCCTTGAGATTACGTTTTGTTATTCTGTCAGTTGTGAAACTGAAGCGAGTTAGAAGCTATCGCGCTCATCAAAGATg                                                                         >  2:134509/1‑149 (MQ=255)
                                                                     tGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATg                                                                         >  8:323988/1‑149 (MQ=255)
                                                                     tGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATg                                                                         >  6:300987/1‑149 (MQ=255)
                                                                      gAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCa                                                                                                                        >  3:8021/1‑101 (MQ=255)
                                                                      gAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCa                                                                                                                        <  4:8021/101‑1 (MQ=255)
                                                                       aaaaaGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGc                                                                       >  3:47980/1‑149 (MQ=255)
                                                                       aaaaaGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGc                                                                       >  1:156177/1‑149 (MQ=255)
                                                                        aaaaGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGc                                                                       <  8:34414/148‑1 (MQ=255)
                                                                               gCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCttt                                                                >  4:286848/1‑148 (MQ=255)
                                                                                cGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTga                                                              <  5:257169/149‑1 (MQ=255)
                                                                                cGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTga                                                              <  5:18576/149‑1 (MQ=255)
                                                                                 gAAAATCATCGCCGAACACAGAAGCAATCCATAACAGAAACCTGAAACAAAAACCGGATTCAGTTGATCTTCCGTTTTTTTATTCTGTCAGTTGTGAAACTTAATCTATTTAGTCTCTATCGATCTCATCAAATATTGCTCGCTTTga                                                              <  1:134509/148‑1 (MQ=255)
                                                                                    aaTCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGata                                                          <  3:286848/149‑1 (MQ=255)
                                                                                           gccgcAAACAGCAGCAATACATAACCGAAACATGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTc                                                    <  8:167324/148‑1 (MQ=255)
                                                                                           gccgcAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTc                                                    <  5:300987/148‑1 (MQ=255)
                                                                                            ccgcAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTc                                                    >  4:296177/1‑147 (MQ=255)
                                                                                            ccgcAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCaa                                                  <  5:262831/149‑1 (MQ=255)
                                                                                                       agcaATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACa                                       >  7:298770/1‑149 (MQ=255)
                                                                                                        gcaATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACAt                                      <  5:197646/149‑1 (MQ=255)
                                                                                                          aaTACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCaa                                                  >  8:41639/1‑135 (MQ=255)
                                                                                                          aaTACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCaa                                                  <  7:41639/135‑1 (MQ=255)
                                                                                                               aTAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTcc                               >  2:364484/1‑149 (MQ=255)
                                                                                                               aTAACAGAAACCTGAAACACAAAACGGCAGCCCTTGACCTCCCGTGTTTTTATTATCTCAGTAGTGAAACTGAAGCCATTTAGACGCTATCGATCTCATCAAATATGGCTCGATTAGAGATATTCCTCAAGTAAAAAAACATCTCTTcc                               >  5:261657/1‑149 (MQ=255)
                                                                                                                tAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCt                              <  1:271292/149‑1 (MQ=255)
                                                                                                                 aaCAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTg                             >  7:265979/1‑149 (MQ=255)
                                                                                                                    aGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGa                          >  8:242991/1‑149 (MQ=255)
                                                                                                                       aaCCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGAttt                       >  1:108085/1‑149 (MQ=255)
                                                                                                                       aaCCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGAttt                       <  8:265979/149‑1 (MQ=255)
                                                                                                                        aCCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTg                             >  4:339156/1‑142 (MQ=255)
                                                                                                                        aCCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTg                             <  3:339156/142‑1 (MQ=255)
                                                                                                                           tGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTca                   >  7:105016/1‑149 (MQ=255)
                                                                                                                               acacAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACaa                >  2:305658/1‑148 (MQ=255)
                                                                                                                                cacaAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTc                                >  7:47380/1‑131 (MQ=255)
                                                                                                                                cacaAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTc                                <  8:47380/131‑1 (MQ=255)
                                                                                                                                 acaAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACaaaa              >  2:208943/1‑148 (MQ=255)
                                                                                                                                  caAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACaaaaaa            <  8:298770/149‑1 (MQ=255)
                                                                                                                                  aaaaaCGGCAGCTCAAGGGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACaaaaaa            <  1:315063/148‑1 (MQ=255)
                                                                                                                                   aaaaCGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGCTATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACaaaa              >  3:173712/1‑146 (MQ=255)
                                                                                                                                   aaaaCGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACaaaa              >  7:315615/1‑146 (MQ=255)
                                                                                                                                     aaCGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGAt         >  6:27980/1‑149 (MQ=255)
                                                                                                                                          cAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGtt    <  8:193952/149‑1 (MQ=255)
                                                                                                                                          cAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGtt    <  2:156177/149‑1 (MQ=255)
                                                                                                                                          cAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGtt    <  4:47980/149‑1 (MQ=255)
                                                                                                                                             atcTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTGa  >  6:102347/3‑148 (MQ=255)
                                                                                                                                                ttGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTaaaa                                            <  7:288069/103‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
CCAGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTGA  >  NC_000913/2736769‑2737057

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