Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,177,646 T→G 5.0% N109K (AAT→AAG lolE → lipoprotein‑releasing system transmembrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,177,6460TG 5.0% 435.7 / 12.6 119N109K (AAT→AAGlolElipoprotein‑releasing system transmembrane protein
Reads supporting (aligned to +/- strand):  ref base T (56/57);  new base G (3/3);  total (59/60)
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.60e-01

Only 100 of 120 total aligned reads displayed.

TGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGC  >  NC_000913/1177498‑1177789
                                                                                                                                                    |                                                                                                                                               
     gtggtgCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTcg                                                                                                                                            <  3:179002/149‑1 (MQ=255)
      tggtgCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTcgc                                                                                                                                           <  4:59802/149‑1 (MQ=255)
      tggtgCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTcgc                                                                                                                                           <  1:145553/149‑1 (MQ=255)
       ggtgCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTcgcg                                                                                                                                          <  3:1883853/149‑1 (MQ=255)
        gtgCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCgcgc                                                                                                                                         <  2:1884051/149‑1 (MQ=255)
                   tGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGa                                                                                                                              >  3:539915/1‑149 (MQ=255)
                   tGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAAGATTCgcgcca                                                                                                                                       >  2:810898/1‑138 (MQ=255)
                    gAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAgg                                                                                                                                 <  2:1582209/145‑1 (MQ=255)
                        tCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGc                                                                                                                         >  2:1464747/1‑149 (MQ=255)
                        tCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGc                                                                                                                         <  3:1700895/149‑1 (MQ=255)
                             gCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTaa                                                                                                                    <  1:1236268/149‑1 (MQ=255)
                              cGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAAc                                                                                                                   <  2:1169790/149‑1 (MQ=255)
                                gTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAgg                                                                                                                           >  4:1782804/1‑139 (MQ=255)
                                gTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAgg                                                                                                                           <  3:1782804/139‑1 (MQ=255)
                                  ggATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGc                                                                                                               <  1:1464747/149‑1 (MQ=255)
                                  ggATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGc                                                                                                               >  3:87516/1‑149 (MQ=255)
                                          cGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAg                                                                                                          >  2:1693361/1‑146 (MQ=255)
                                          cGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAg                                                                                                          <  1:1693361/146‑1 (MQ=255)
                                             ggACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAg                                                                                                    >  3:627604/1‑149 (MQ=255)
                                              gACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGc                                                                                                   >  2:1882628/1‑149 (MQ=255)
                                                 tAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTc                                                                                                <  3:1440023/149‑1 (MQ=255)
                                                     tGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAg                                                                                            >  1:1185820/1‑149 (MQ=255)
                                                      ggCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAgc                                                                                           <  4:714324/149‑1 (MQ=255)
                                                      ggCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAgc                                                                                           >  1:125009/1‑149 (MQ=255)
                                                         aGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCTATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCa                                                                                        >  3:809571/1‑149 (MQ=255)
                                                         aGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCa                                                                                        >  1:653608/1‑149 (MQ=255)
                                                          ggAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGaa                                                                                                       >  1:515333/1‑133 (MQ=255)
                                                          ggAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGaa                                                                                                       <  2:515333/133‑1 (MQ=255)
                                                          ggAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCAt                                                                                       >  3:1096089/1‑149 (MQ=255)
                                                            aaGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCt                                                                                               >  3:779858/1‑139 (MQ=255)
                                                            aaGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCt                                                                                               <  4:779858/139‑1 (MQ=255)
                                                            aaGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTa                                                                                     >  2:1745825/1‑149 (MQ=255)
                                                             aGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACa                                                                                                           <  3:1698253/126‑1 (MQ=255)
                                                             aGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACa                                                                                                           >  4:1698253/1‑126 (MQ=255)
                                                              gCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTAcc                                                                                   >  1:874032/1‑149 (MQ=255)
                                                                  tGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCAt                                                                                       >  3:1573120/1‑141 (MQ=255)
                                                                    gATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACa                                                                                                     >  3:40924/1‑125 (MQ=255)
                                                                    gATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACa                                                                                                     <  4:40924/125‑1 (MQ=255)
                                                                    gATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGtt                                                                             <  2:1185820/149‑1 (MQ=255)
                                                                        aCGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTc                                                                         <  1:916373/149‑1 (MQ=255)
                                                                            gCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAgg                                                                     <  2:264048/149‑1 (MQ=255)
                                                                              aaaaaGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGt                                                                              <  4:1739462/138‑1 (MQ=255)
                                                                              aaaaaGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGt                                                                              >  3:1739462/1‑138 (MQ=255)
                                                                               aaaaGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTa                                                                                     <  3:746523/130‑1 (MQ=255)
                                                                               aaaaGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTa                                                                                     >  4:746523/1‑130 (MQ=255)
                                                                                      ccGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTg                                                           <  4:539915/149‑1 (MQ=255)
                                                                                        gggTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGgc                                                         >  2:1753446/1‑149 (MQ=255)
                                                                                        gggTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGgc                                                         <  1:1753446/149‑1 (MQ=255)
                                                                                              tgccgctgccgcGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGgcg                                                        <  2:91414/144‑1 (MQ=255)
                                                                                              tgccgctgccgcGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGgcg                                                        >  1:91414/1‑144 (MQ=255)
                                                                                              tgccgctgccgcGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAAtt                                                   <  1:1745825/149‑1 (MQ=255)
                                                                                                ccgctgccgcGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAAtttt                                                 <  1:1882628/149‑1 (MQ=255)
                                                                                                   ctgccgcGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTaaa                                              <  2:653608/149‑1 (MQ=255)
                                                                                                    tgccgcGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAg                                             >  2:364135/1‑149 (MQ=255)
                                                                                                            cGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGaa                                     >  2:1395978/1‑149 (MQ=255)
                                                                                                              tataTCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGt                                                                                                                       <  3:1232575/65‑1 (MQ=255)
                                                                                                              tataTCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGt                                                                                                                       >  4:1232575/1‑65 (MQ=255)
                                                                                                               atatCAATTTCACCGGGCTGGCGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTTAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGCCGCAATTTTAAAGCGGGCGAACAg                                  >  2:274651/1‑149 (MQ=255)
                                                                                                                tatCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGc                                 >  2:1210586/1‑149 (MQ=255)
                                                                                                                 atCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCa                                >  1:1677806/1‑149 (MQ=255)
                                                                                                                  tCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAAt                                                    <  3:648117/128‑1 (MQ=255)
                                                                                                                  tCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAAt                                                    >  4:648117/1‑128 (MQ=255)
                                                                                                                  tCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCggg                                         <  1:1141856/139‑1 (MQ=255)
                                                                                                                  tCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCggg                                         >  4:1209233/1‑139 (MQ=255)
                                                                                                                  tCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCggg                                         <  3:1209233/139‑1 (MQ=255)
                                                                                                                  tCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCaa                               >  2:311828/1‑149 (MQ=255)
                                                                                                                   cAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCaaa                              <  2:874032/149‑1 (MQ=255)
                                                                                                                    aaTTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCggg                                         <  1:1558980/137‑1 (MQ=255)
                                                                                                                    aaTTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCggg                                         >  2:1558980/1‑137 (MQ=255)
                                                                                                                    aaTTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAAt                             >  3:1700484/1‑149 (MQ=255)
                                                                                                                      tttCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTa                           >  1:494430/1‑149 (MQ=255)
                                                                                                                      tttCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTa                           >  1:1450313/1‑149 (MQ=255)
                                                                                                                      tttCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTa                           <  4:1096089/149‑1 (MQ=255)
                                                                                                                        tCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGttt                                                                            >  1:1173663/1‑98 (MQ=255)
                                                                                                                        tCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGttt                                                                            <  2:1173663/98‑1 (MQ=255)
                                                                                                                         cACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCa                                <  4:151226/141‑1 (MQ=255)
                                                                                                                          aCCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTAcc                                                                                   >  3:280604/1‑89 (MQ=255)
                                                                                                                          aCCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTAcc                                                                                   <  4:280604/89‑1 (MQ=255)
                                                                                                                          aCCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTatcat                       >  1:825400/1‑149 (MQ=255)
                                                                                                                              ggCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGc                   <  2:494430/149‑1 (MQ=255)
                                                                                                                                   gtggAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTgg                                                          >  4:860931/1‑105 (MQ=255)
                                                                                                                                   gtggAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTgg                                                          <  3:860931/105‑1 (MQ=255)
                                                                                                                                     ggAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAAt                                                    <  4:1425091/109‑1 (MQ=255)
                                                                                                                                     ggAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAAt                                                    >  3:1425091/1‑109 (MQ=255)
                                                                                                                                     ggAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGc             >  4:1640278/1‑148 (MQ=255)
                                                                                                                                     ggAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGc             <  3:1640278/148‑1 (MQ=255)
                                                                                                                                       aaaGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTg          <  2:825400/149‑1 (MQ=255)
                                                                                                                                           ttgcgcGAAGATTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTatca                        <  1:1686791/129‑1 (MQ=255)
                                                                                                                                           ttgcgcGAAGATTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTatca                        >  2:1686791/3‑131 (MQ=255)
                                                                                                                                            tgcgcGAAGATTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGc        <  2:346273/145‑1 (MQ=255)
                                                                                                                                            tgcgcGAAGATTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGc        >  1:346273/2‑146 (MQ=255)
                                                                                                                                             gcgcGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTaaa                                              >  1:828989/1‑107 (MQ=255)
                                                                                                                                             gcgcGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTaaa                                              <  2:828989/107‑1 (MQ=255)
                                                                                                                                             gcgcGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGc                                        <  4:497460/113‑1 (MQ=255)
                                                                                                                                             gcgcGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGc                                        >  3:497460/1‑113 (MQ=255)
                                                                                                                                             gcgcGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGAt    <  4:1700484/149‑1 (MQ=255)
                                                                                                                                              cgcgAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATg   >  2:1739505/1‑149 (MQ=255)
                                                                                                                                               gcgAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAAtt                                                   >  2:1258176/1‑100 (MQ=255)
                                                                                                                                               gcgAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAAtt                                                   <  1:1258176/100‑1 (MQ=255)
                                                                                                                                               gcgAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATgc  <  2:125009/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGC  >  NC_000913/1177498‑1177789

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