Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,664,599:1 +T 9.2% coding (94/615 nt) dmsD → twin‑argninine leader‑binding protein for DmsA and TorA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,664,5991.T 9.2% 257.6 / 21.6 65coding (94/615 nt)dmsDtwin‑argninine leader‑binding protein for DmsA and TorA
Reads supporting (aligned to +/- strand):  ref base . (33/26);  new base T (3/3);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GTTAACAAATTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGC  >  NC_000913/1664459‑1664747
                                                                                                                                             |                                                                                                                                                    
gTTAACAAATTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAG‑GTGCGGCTTCCGCActctct                                                                                                                                                  >  1:1524136/1‑140 (MQ=255)
         ttGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCgg                                                                                                                                          <  2:1643171/144‑1 (MQ=255)
         ttGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCgg                                                                                                                                          >  1:1643171/1‑144 (MQ=255)
               cGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGc                                                                                                                                            >  1:685642/1‑136 (MQ=255)
               cGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGc                                                                                                                                            <  2:685642/136‑1 (MQ=255)
                        ccGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAG‑GTGCGGCTTCCGCACTCTCTgga                                                                                                                                               <  3:184090/124‑2 (MQ=255)
                        ccGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAG‑GTGCGGCTTCCGCACTCTCTgga                                                                                                                                               >  4:184090/1‑123 (MQ=255)
                           cacaTTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAG‑GTGCGGCTTCCGCACTCTCTgga                                                                                                                                               <  3:352476/121‑2 (MQ=255)
                           cacaTTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAG‑GTGCGGCTTCCGCACTCTCTgga                                                                                                                                               >  4:352476/1‑120 (MQ=255)
                                           aGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTAc                                                                                                   >  4:548376/1‑149 (MQ=255)
                                             gAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCa                                                                                                 >  4:20780/1‑149 (MQ=255)
                                                   cccATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGt                                                                                                            <  1:152908/132‑1 (MQ=255)
                                                   cccATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGt                                                                                                            >  2:152908/1‑132 (MQ=255)
                                                    ccATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGt                                                                                          >  3:1554352/1‑149 (MQ=255)
                                                        ttttCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCAtt                                                                                      <  3:548376/149‑1 (MQ=255)
                                                          ttCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAg                                                                                    <  4:1554352/149‑1 (MQ=255)
                                                                   aGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCAt                                                                                       <  1:1089503/137‑1 (MQ=255)
                                                                   aGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCAt                                                                                       >  2:1089503/1‑137 (MQ=255)
                                                                    gATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGt                                                                          <  3:20780/149‑1 (MQ=255)
                                                                              cTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCtttt                                                                >  3:1507179/1‑149 (MQ=255)
                                                                                             tgtTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGc                                                                                            <  2:1387855/106‑1 (MQ=255)
                                                                                             tgtTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGc                                                                                            >  1:1387855/1‑106 (MQ=255)
                                                                                               ttGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGAc                                                                                                                            >  3:19363/1‑72 (MQ=255)
                                                                                               ttGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGAc                                                                                                                            <  4:19363/72‑1 (MQ=255)
                                                                                                 gggCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAg                                                                                                             >  2:1540282/1‑85 (MQ=255)
                                                                                                 gggCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAg                                                                                                             <  1:1540282/85‑1 (MQ=255)
                                                                                                   gcgcgcTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCa                                                                                                                                 <  3:175138/63‑1 (MQ=255)
                                                                                                   gcgcgcTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCa                                                                                                                                 >  4:175138/1‑63 (MQ=255)
                                                                                                     gcgcTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCa                                         >  3:1894094/1‑149 (MQ=255)
                                                                                                       gcTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACg                                       >  1:1549931/1‑149 (MQ=255)
                                                                                                       gcTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACg                                       >  1:76170/1‑149 (MQ=255)
                                                                                                          gTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGccc                                    >  2:665110/1‑149 (MQ=255)
                                                                                                          gTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGccc                                    >  4:252659/1‑149 (MQ=255)
                                                                                                          gTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGccc                                    >  4:564414/1‑149 (MQ=255)
                                                                                                                ttaCGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTg                                                                                                           <  1:1617916/72‑1 (MQ=255)
                                                                                                                ttaCGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTg                                                                                                           >  2:1617916/1‑72 (MQ=255)
                                                                                                                 taCGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTg                             <  4:1507179/149‑1 (MQ=255)
                                                                                                                     aaccAGAG‑GTGCGGCTTCCGCACTCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGa                                                  <  2:955469/122‑1 (MQ=38)
                                                                                                                     aaccAGAG‑GTGCGGCTTCCGCACTCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGa                                                  >  1:955469/3‑124 (MQ=38)
                                                                                                                      tCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGc                        <  3:564414/149‑1 (MQ=255)
                                                                                                                      tCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGc                        >  1:1018841/1‑149 (MQ=255)
                                                                                                                      tCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGc                        >  1:1517468/1‑149 (MQ=255)
                                                                                                                      tCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGc                        <  1:665110/149‑1 (MQ=255)
                                                                                                                        cAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGa                                                                        <  3:942896/99‑1 (MQ=255)
                                                                                                                        cAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGa                                                                        >  4:942896/1‑99 (MQ=255)
                                                                                                                         agagagTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGa                                                            <  1:530828/110‑1 (MQ=255)
                                                                                                                         agagagTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGa                                                            >  2:530828/1‑110 (MQ=255)
                                                                                                                            gagTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgt                  <  3:252659/149‑1 (MQ=255)
                                                                                                                             agTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgtt                 >  3:1487747/1‑149 (MQ=255)
                                                                                                                             agTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgtt                 <  4:1894094/149‑1 (MQ=255)
                                                                                                                             agTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgtt                 <  2:1549931/149‑1 (MQ=255)
                                                                                                                              gTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgttt                <  2:1517468/149‑1 (MQ=255)
                                                                                                                               tGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGtttgtttg               >  1:194566/1‑149 (MQ=255)
                                                                                                                                gCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTttgtttgt              <  2:76170/149‑1 (MQ=255)
                                                                                                                                 cGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTc             <  2:194566/149‑1 (MQ=255)
                                                                                                                                 cGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTc             >  3:1579536/1‑149 (MQ=255)
                                                                                                                                  ggAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCtttt                                                                >  2:812168/1‑97 (MQ=255)
                                                                                                                                  ggAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCtgtt                                                                <  1:812168/97‑1 (MQ=255)
                                                                                                                                   gAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgt                  <  4:1062035/142‑1 (MQ=255)
                                                                                                                                   gAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgt                  >  3:1062035/1‑142 (MQ=255)
                                                                                                                                      gCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGccc        <  2:1018841/149‑1 (MQ=255)
                                                                                                                                        cGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCCGCGTTTGTTTGTCGGCCGGt      >  4:255833/1‑149 (MQ=255)
                                                                                                                                         gCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTTGGCCCGTg     >  2:1532774/1‑149 (MQ=255)
                                                                                                                                         gCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCTGCCCGTg     >  1:640012/1‑149 (MQ=255)
                                                                                                                                         gCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTg     >  3:1644151/1‑149 (MQ=255)
                                                                                                                                         gCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTg     >  3:1527900/1‑149 (MQ=255)
                                                                                                                                            c‑cTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGAAACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGGCGGCCAGTAGGc  >  4:1464132/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
GTTAACAAATTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGC  >  NC_000913/1664459‑1664747

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