Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,951,101 A→G 5.7% L90L (TTA→TTG yecD → isochorismatase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,951,1010AG 5.7% 215.9 / 2.5 70L90L (TTA→TTGyecDisochorismatase family protein
Reads supporting (aligned to +/- strand):  ref base A (32/34);  new base G (2/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

GGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGC  >  NC_000913/1950954‑1951247
                                                                                                                                                   |                                                                                                                                                  
ggcgAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAg                                                                                                                                                   <  3:117119/149‑1 (MQ=255)
ggcgAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAg                                                                                                                                                   <  2:283000/149‑1 (MQ=255)
ggcgAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAg                                                                                                                                                   <  2:278306/149‑1 (MQ=255)
 gcgAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAgg                                                                                                                                                  <  8:264197/149‑1 (MQ=255)
  cgAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGt                                                                                                                                                 <  4:140151/149‑1 (MQ=255)
    aaaTTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTg                                                                                                                                                >  8:234498/1‑148 (MQ=255)
         tCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAAcc                                                                                                                                           <  4:115301/148‑1 (MQ=255)
                aGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCCACATCCTGCTGCATTAGGTGCAACCGACAGCg                                                                                                                                    >  6:64363/1‑148 (MQ=255)
                     tCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCg                                                                                                                              <  4:330367/149‑1 (MQ=255)
                         gcggtgtTTCTGGTGCGGGTTGGCTGGTGTGCCGTTTTCGCTGAAGCATAAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCGGTAAATTGGTGGCAAGATCCTGCTGCATTAGGTGCAACCGTCAGCGATATCGAAat                                                                                                                          <  2:323792/146‑1 (MQ=255)
                         cccGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAat                                                                                                                          >  2:296456/1‑149 (MQ=255)
                         cccGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAat                                                                                                                          <  4:42856/149‑1 (MQ=255)
                          ccGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAat                                                                                                                          >  8:284959/1‑148 (MQ=255)
                           cGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAatca                                                                                                                        <  6:317290/149‑1 (MQ=255)
                           cGTGTGTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCCGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAatca                                                                                                                        <  2:299991/149‑1 (MQ=255)
                            gtgtTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAatcat                                                                                                                       >  1:293576/1‑149 (MQ=255)
                             tgtTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAatcatc                                                                                                                      <  5:185913/149‑1 (MQ=255)
                                  cTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAAc                                                                                                                  <  4:77730/148‑1 (MQ=255)
                                     gTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTc                                                                                                               >  8:149704/1‑148 (MQ=255)
                                       gcgcGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGACATCGAAATCATCAAACGTCAAt                                                                                                            <  4:206580/149‑1 (MQ=255)
                                          cgTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATggg                                                                                                         <  5:73478/149‑1 (MQ=255)
                                          cgTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATggg                                                                                                         <  6:65089/149‑1 (MQ=255)
                                           gTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATgggg                                                                                                        <  1:70599/149‑1 (MQ=255)
                                                 tGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGtt                                                                                                  >  3:95175/1‑149 (MQ=255)
                                                   gTCTGCCGAGTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGGATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGttt                                                                                                 >  6:289756/1‑148 (MQ=255)
                                                      tGCCGAATACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTtacg                                                                                             >  4:330208/1‑149 (MQ=255)
                                                      tGCAGATTACGCTGAAGCATTAAAACACCCGTTCGATGCTCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGGAATCATCAAACTTCAATTGGGTGCGTTTtact                                                                                             >  2:180568/1‑148 (MQ=255)
                                                        ccGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTGGGTGCAACCGACAGCGATATCGAAATAATCAAACGTCAATGGGTTGCGTTTtacgg                                                                                            >  1:18422/1‑148 (MQ=255)
                                                        ccGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTGGGGGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGGGCGTGTtacggt                                                                                           >  1:12817/1‑149 (MQ=255)
                                                          gATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTtacggtac                                                                                         >  2:292614/1‑149 (MQ=255)
                                                                gCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCt                                                                                   >  7:80683/1‑149 (MQ=255)
                                                                gCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCt                                                                                   >  1:161506/1‑149 (MQ=255)
                                                                 cTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTTCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTg                                                                                  <  2:293576/149‑1 (MQ=255)
                                                                 cTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGGGGCAACATCCTGCTGCATTAGGGGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTg                                                                                  <  3:177017/149‑1 (MQ=255)
                                                                  tGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGtt                                                                                                  <  6:254252/132‑1 (MQ=255)
                                                                  tGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGtt                                                                                                  >  5:254252/1‑132 (MQ=255)
                                                                  tGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTg                                                                                  >  6:198738/1‑148 (MQ=255)
                                                                    aaGCATTAAAACATCCTGTCGATTCGCCCTCACCGGCTAAAGTGTTGACCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGTATCTgg                                                                                 <  3:173414/147‑1 (MQ=255)
                                                                    aaGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGa                                                                                <  7:234498/148‑1 (MQ=255)
                                                                    aaGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATAGGTGGCACCAGCGTGCTGCAGAAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGGGTTACGGTACGGATCTGGAg                                                                               >  4:219490/1‑149 (MQ=255)
                                                                       cATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTg                                                                            >  7:238950/1‑149 (MQ=255)
                                                                           aaaaCAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAGCCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAAt                                                                        >  2:21712/1‑149 (MQ=255)
                                                                            aaaCAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAAt                                                                        >  6:154271/1‑148 (MQ=255)
                                                                             aaCAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTa                                                                      <  1:296456/149‑1 (MQ=255)
                                                                                aGCCGGTCGATGCGCCCTCCCCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTAcgc                                                                   >  3:60742/1‑149 (MQ=255)
                                                                                        gATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAgg                                                           <  7:284959/149‑1 (MQ=255)
                                                                                        gATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGGCAATGGGGTGCGTTTTACGGTACGTATCTGGAGTTGCAATTACGCCGCCGAg                                                            <  1:121765/148‑1 (MQ=255)
                                                                                              cccTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGa                                                     <  4:127609/149‑1 (MQ=255)
                                                                                                   aCCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATcatgatta                                                >  2:283731/1‑142 (MQ=255)
                                                                                                       gCTAAAGTGTTGCCCGAAAATTGGTGGAAACAACCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGAGTTTTACGGTACGGATCTGGAGTCGCAATTACGCCGCCGAGGTATCGATACAATAgt                                            >  8:29412/1‑149 (MQ=255)
                                                                                                        cTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAgtg                                           >  6:134420/1‑149 (MQ=255)
                                                                                                         tAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAgtgt                                          >  3:216089/1‑149 (MQ=255)
                                                                                                           aaGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTa                                        <  7:149704/149‑1 (MQ=255)
                                                                                                           aaGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAAGAGTGTTa                                        >  6:241892/1‑149 (MQ=255)
                                                                                                                ttGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTAtgtg                                    <  8:80683/148‑1 (MQ=255)
                                                                                                                    ccGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGAtc                               >  3:164964/1‑149 (MQ=255)
                                                                                                                      gAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCTATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGAtctc                             <  7:306838/149‑1 (MQ=255)
                                                                                                                         aaTTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGAc                          <  1:292614/149‑1 (MQ=255)
                                                                                                                               tggCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAatat                    >  4:320829/1‑149 (MQ=255)
                                                                                                                                 gCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGTATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCg                  <  4:95175/149‑1 (MQ=255)
                                                                                                                                   aaCATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGgt                >  4:117318/1‑149 (MQ=255)
                                                                                                                                   aaCATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGgt                <  2:161506/149‑1 (MQ=255)
                                                                                                                                    aCATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGgtg               >  4:84792/1‑149 (MQ=255)
                                                                                                                                        ccTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTg            >  7:268569/1‑148 (MQ=255)
                                                                                                                                         ctgctgCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGaa          <  5:154271/149‑1 (MQ=255)
                                                                                                                                               cATTGGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGAGCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCa      <  2:12817/147‑1 (MQ=255)
                                                                                                                                               cATTGGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAAGGGGGTGCGTTTTACGGTACTGATCTGGAGTTTCAATTAAGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCa      <  2:18422/147‑1 (MQ=255)
                                                                                                                                               cATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCAc     <  6:58073/148‑1 (MQ=255)
                                                                                                                                                aTTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTAAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACTCACCCGAGGTATCGATACAATAGTGTTATGTGGTATCTCGACCAATATCGGTGTTTAATCAAcag   >  2:253064/1‑147 (MQ=255)
                                                                                                                                                 ttAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTATTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCTACCAATATCGGTGTTGAATCCAccgc  >  5:106095/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGC  >  NC_000913/1950954‑1951247

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