Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,439,589 G→C 7.6% intergenic (‑80/+27) yhdN ← / ← rplQ DUF1992 family protein/50S ribosomal subunit protein L17

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,439,5890GC 7.6% 216.8 / 8.9 66intergenic (‑80/+27)yhdN/rplQDUF1992 family protein/50S ribosomal subunit protein L17
Reads supporting (aligned to +/- strand):  ref base G (30/31);  new base C (3/2);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.69e-01

GTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCT  >  NC_000913/3439441‑3439735
                                                                                                                                                    |                                                                                                                                                   
gTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCgg                                                                                                                                                     >  3:166624/1‑149 (MQ=255)
gTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCgg                                                                                                                                                     >  7:101001/1‑149 (MQ=255)
      aaaCTCACCTTTCGCTTGCGCTTCTTCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGTGGATACTACTGGTATAAAAAAACCCGCCGGGGCggg                                                                                                                                               >  7:212144/1‑149 (MQ=255)
       aaCTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑t                                                                                                                                             >  2:140344/1‑149 (MQ=255)
        aCTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑tt                                                                                                                                            >  2:245436/1‑149 (MQ=255)
            aCCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑tttttt                                                                                                                                        <  7:47896/149‑1 (MQ=255)
            aCCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑tttttt                                                                                                                                        >  7:150844/1‑149 (MQ=255)
            aCCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑ttttt                                                                                                                                         <  1:140344/148‑1 (MQ=255)
            aCCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGG‑tttttt                                                                                                                                        <  3:83634/149‑1 (MQ=255)
              cTTTCGCTTGCGCTTCTGCTCTATGGCGCTCTGCCCACTGGTCCAGTAACCACATAAATACTCCAGGGGCTGAACAAAACGCGGACAGCTATTGTAGCTAAGTGGGGCTACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTa                                                                                                                                      <  3:130076/149‑1 (MQ=255)
              cTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACAGAAATACTCCAGGGGATGAACAATAAGAGTCGAGCTATTGTAGATACGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTGTTa                                                                                                                                      >  6:105945/1‑149 (MQ=255)
                    cTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGc                                                                                                                                >  3:84121/1‑149 (MQ=255)
                       gcgcTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCtt                                                                                                                              >  6:23505/1‑148 (MQ=255)
                        cgcTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCtt                                                                                                                              >  7:143545/1‑147 (MQ=255)
                           ttCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGAt                                                                                                                         <  4:60311/149‑1 (MQ=255)
                             cTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGAtt                                                                                                                        <  3:222592/148‑1 (MQ=255)
                                 tataTGGCGCTCTGCCCACTGATCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGTTAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTActct                                                                                                                   >  6:207406/1‑149 (MQ=255)
                                   tatGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTg                                                                                                                  <  4:200254/148‑1 (MQ=255)
                                   tatGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTg                                                                                                                  <  8:2928/148‑1 (MQ=255)
                                     tGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTgcag                                                                                                               <  3:79808/149‑1 (MQ=255)
                                       gcgcTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTgcagca                                                                                                             <  8:228783/149‑1 (MQ=255)
                                             tGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGTGATACTACGTGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTAAGTTGCTTCAGATTACTCTGCAGCAGCTTCt                                                                                                       >  2:202599/1‑149 (MQ=255)
                                             tGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCtt                                                                                                         <  2:85753/147‑1 (MQ=255)
                                               cccACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTc                                                                                                        <  2:145410/146‑1 (MQ=255)
                                                  aCTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCttt                                                                                                  <  1:245436/149‑1 (MQ=255)
                                                      gTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTtctc                                                                                               <  4:166624/148‑1 (MQ=255)
                                                         aaGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGTTAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGaa                                                                                           <  4:20420/149‑1 (MQ=255)
                                                            tAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGa                                                                                        <  5:192927/149‑1 (MQ=255)
                                                              aCCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGAt                                                                                       <  8:150844/148‑1 (MQ=255)
                                                               ccACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCa                                                                                     <  8:101001/149‑1 (MQ=255)
                                                                 acaTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAc                                                                                   <  4:84121/149‑1 (MQ=255)
                                                                   aTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAc                                                                                   >  5:76211/1‑147 (MQ=255)
                                                                   aTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCa                                                                                 <  8:179997/149‑1 (MQ=255)
                                                                     aaaTACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAACAg                                                                                >  5:195241/1‑148 (MQ=255)
                                                                      aaTACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCt                                                                              >  7:104802/1‑149 (MQ=255)
                                                                           tCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATg                                                                         >  2:10194/1‑149 (MQ=255)
                                                                            ccAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGt                                                                        <  5:23505/149‑1 (MQ=255)
                                                                                gggATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGc                                                                    >  1:77786/1‑149 (MQ=255)
                                                                                  gATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCACCTTCTGCGTTCTCTGAACTATAAACCATCTCGATGt                                                                        >  7:111866/1‑143 (MQ=255)
                                                                                    tGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATc                                                                >  1:82447/1‑149 (MQ=255)
                                                                                          aaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAACCCCGGCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGATCGATCAACCAGCTCGATGTAAGCCATCGGcgcg                                                          <  3:227864/149‑1 (MQ=255)
                                                                                           aaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGt                                                         <  2:77786/149‑1 (MQ=255)
                                                                                                 tACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGc                                                   >  7:60744/1‑149 (MQ=255)
                                                                                                 tACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGc                                                   <  7:87059/149‑1 (MQ=255)
                                                                                                 tACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGc                                                   >  8:89137/1‑149 (MQ=255)
                                                                                                    aGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCAGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATAACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCAACGGCGCGTTGTCGCATg                                                >  1:196319/1‑149 (MQ=255)
                                                                                                      cTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCa                                              >  4:103275/1‑149 (MQ=255)
                                                                                                        aTTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACg                                            >  2:141632/1‑149 (MQ=255)
                                                                                                          tGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACgg                                           >  6:170461/1‑148 (MQ=255)
                                                                                                            tAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGaa                                         <  7:89137/148‑1 (MQ=255)
                                                                                                               aTAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCca                                     <  1:10194/149‑1 (MQ=255)
                                                                                                                      gggATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCa                              <  5:170461/149‑1 (MQ=255)
                                                                                                                             tacGGGTATAAAAAAACCCGCCGGGGCGGGCTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGAt                                                                          <  6:198321/99‑1 (MQ=255)
                                                                                                                             tacGGGTATAAAAAAACCCGCCGGGGCGGGCTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGAt                                                                          >  5:198321/1‑99 (MQ=255)
                                                                                                                              acGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGGCGCCTGCACGGAAGCcac                                    >  1:161160/1‑135 (MQ=255)
                                                                                                                              acGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGCCGCCTGCACGGAAGCcac                                    <  2:161160/135‑1 (MQ=255)
                                                                                                                               cGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGTACGGAAGCCACACTTCAGAATACGa                      <  6:76211/148‑1 (MQ=255)
                                                                                                                                  gTATAAAAAAACCCGCCTGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGTAGCCACACTTCAGAATACGAgtgt                  <  1:202599/149‑1 (MQ=255)
                                                                                                                                    ataAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTag                >  1:166549/1‑148 (MQ=255)
                                                                                                                                      aaaaaaaCCCGCCCCGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTg                                                                                                                  <  7:110783/49‑1 (MQ=21)
                                                                                                                                      aaaaaaaCCCGCCCCGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTg                                                                                                                  >  8:110783/1‑49 (MQ=21)
                                                                                                                                        aaaaaCCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccac            >  6:2246/1‑149 (MQ=255)
                                                                                                                                        aaaaaCCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccac            >  1:161560/1‑149 (MQ=255)
                                                                                                                                        aaaaaCCCGCCCCGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGGACGGAAGCCACACTTCAGAATACGAGTGTcaccac            >  3:56105/1‑149 (MQ=255)
                                                                                                                                         aaaaCCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccacc           >  3:22918/1‑149 (MQ=255)
                                                                                                                                            aCCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCgg         >  2:125569/1‑148 (MQ=255)
                                                                                                                                             cccGCCCCGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCa       >  4:179521/1‑149 (MQ=255)
                                                                                                                                                  cgGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCt  <  3:103275/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGG‑TTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCT  >  NC_000913/3439441‑3439735

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: ‑