Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,177,646 T→G 7.5% 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 7.5% 304.8 / 17.0 93N109K (AAT→AAGlolElipoprotein‑releasing system transmembrane protein
Reads supporting (aligned to +/- strand):  ref base T (44/42);  new base G (5/2);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01

TGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCG  >  NC_000913/1177498‑1177790
                                                                                                                                                    |                                                                                                                                                
tGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAAt                                                                                                                                                  >  2:934827/1‑149 (MQ=255)
    ggtggtGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTc                                                                                                                                              <  1:641108/149‑1 (MQ=255)
      tggtgCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTcgc                                                                                                                                            <  2:2618/149‑1 (MQ=255)
      tggtgCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTcgc                                                                                                                                            <  2:814257/149‑1 (MQ=255)
          gCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCaa                                                                                                                                        >  3:275669/1‑149 (MQ=255)
          gCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAAGATTCgcgcca                                                                                                                                        <  1:991956/149‑3 (MQ=255)
           ccGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAAGATTCgcgccac                                                                                                                                       <  4:790447/149‑4 (MQ=255)
                tGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAgg                                                                                                                                  <  4:383405/149‑1 (MQ=255)
                cggTGAAATCGAGGCGGGGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAgg                                                                                                                                  <  2:957994/148‑1 (MQ=255)
                       aTCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAggg                                                                                                                           >  1:741495/1‑149 (MQ=255)
                        tCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGc                                                                                                                          <  4:487516/149‑1 (MQ=255)
                          gAGGCGGTGGATCAGCCGTGGACTAACTTGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCCCCGGGCTGGTGGACAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGt                                                                                                                        >  3:924375/1‑149 (MQ=255)
                           aGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGtt                                                                                                                       <  4:503389/149‑1 (MQ=255)
                             gCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTaa                                                                                                                     <  3:424316/149‑1 (MQ=255)
                                gTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAAccc                                                                                                                  <  1:330640/149‑1 (MQ=255)
                                gTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAAccc                                                                                                                  >  2:844524/1‑149 (MQ=255)
                                     tCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAg                                                                                                                                   <  1:97415/127‑1 (MQ=255)
                                     tCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCCGGTGGAAGGTGGCGCGAATCTTCGCGCAATCCAg                                                                                                                                   >  2:97415/1‑127 (MQ=255)
                                     tCAGCCGTGGACTAACTAGCAGGAAGCACTGGATCACATGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCGTCGCGCAATCCAGGTGAAGGGCGTTAACCCACAAc                                                                                                             >  1:831979/1‑149 (MQ=255)
                                      cAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACa                                                                                                            <  3:25576/149‑1 (MQ=255)
                                      cAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTCGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACa                                                                                                            <  3:921791/149‑1 (MQ=255)
                                       aGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAg                                                                                                           >  1:25422/1‑149 (MQ=255)
                                       aGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAg                                                                                                           >  3:47853/1‑149 (MQ=255)
                                       aGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAg                                                                                                           <  1:1064336/149‑1 (MQ=255)
                                          cGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGaa                                                                                                        <  1:9970/149‑1 (MQ=255)
                                           gTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAAc                                                                                                       <  4:324475/149‑1 (MQ=255)
                                                      ggCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAgc                                                                                            <  4:545148/149‑1 (MQ=255)
                                                          ggAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAg                                                                                                                             <  2:931416/112‑1 (MQ=255)
                                                          ggAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAg                                                                                                                             >  1:931416/1‑112 (MQ=255)
                                                          ggAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCAt                                                                                        >  3:1038736/1‑149 (MQ=255)
                                                          ggAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCAt                                                                                        >  3:759545/1‑149 (MQ=255)
                                                             aGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTAc                                                                                     <  4:309455/149‑1 (MQ=255)
                                                              gCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTAcc                                                                                    <  2:741495/149‑1 (MQ=255)
                                                              gCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTAcc                                                                                    >  4:479748/1‑149 (MQ=255)
                                                              gCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTAcc                                                                                    <  4:47853/149‑1 (MQ=255)
                                                                aCTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCt                                                                                  >  2:249705/1‑149 (MQ=255)
                                                                 cTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTc                                                                                 >  4:112161/1‑149 (MQ=255)
                                                                 cTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTc                                                                                 >  3:696032/1‑149 (MQ=255)
                                                                 cTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTc                                                                                 >  3:899950/1‑149 (MQ=255)
                                                                   ggATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGt                                                                               >  1:792777/1‑149 (MQ=255)
                                                                    gATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGTGTCTGAGCGCATTACCCTCGtt                                                                              >  4:67724/1‑149 (MQ=255)
                                                                     aTCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGttt                                                                             <  1:488406/149‑1 (MQ=255)
                                                                      tCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAgg                                                                                                                            >  1:160999/1‑101 (MQ=255)
                                                                      tCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAgg                                                                                                                            <  2:160999/101‑1 (MQ=255)
                                                                      tCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTg                                                                            >  4:535672/1‑149 (MQ=255)
                                                                      tCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTg                                                                            <  2:116645/149‑1 (MQ=255)
                                                                         cGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTcc                                                                         <  2:91428/149‑1 (MQ=255)
                                                                          gTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCa                                                                        >  1:746437/1‑149 (MQ=255)
                                                                           tGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAg                                                                       >  1:37129/1‑149 (MQ=255)
                                                                           tGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAg                                                                       >  3:570242/1‑149 (MQ=255)
                                                                                aaaGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGttt                                                                             <  3:323905/138‑1 (MQ=255)
                                                                                aaaGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGttt                                                                             >  4:323905/1‑138 (MQ=255)
                                                                                   gTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGTGCGAATCt                                                                                                                                                <  4:637663/68‑1 (MQ=255)
                                                                                   gTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGTGCGAATCt                                                                                                                                                >  3:637663/1‑68 (MQ=255)
                                                                                        gggTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGgc                                                          >  3:352342/1‑149 (MQ=255)
                                                                                         ggTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGgcg                                                         <  4:696032/149‑1 (MQ=255)
                                                                                             ttGCCGCTGTCGCGCTGCACATCAATTTCACCGGCCTGGTGGAAAGCGGCGCGAATGTTCGCGCAATCCAAGTGAAGGGCGTTAACCTGCAAAACGAACAGGGTCTGAGCGCATTACCATCGTTTGTCCAGGGCGATGCGTGGCGCCAt                                                     <  4:347544/149‑1 (MQ=255)
                                                                                             ttGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAAt                                                     <  4:437561/149‑1 (MQ=255)
                                                                                                       cgcgCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCgg                                           <  2:25422/149‑1 (MQ=255)
                                                                                                        gcgcCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCggg                                          <  1:36922/149‑1 (MQ=255)
                                                                                                          gcCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCg                                        >  3:1027632/1‑149 (MQ=255)
                                                                                                            cGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGaa                                      <  1:795533/149‑1 (MQ=255)
                                                                                                             gTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGTGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAAc                                     <  3:67724/149‑1 (MQ=255)
                                                                                                               atatCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCAt                                                                                        >  4:972281/1‑96 (MQ=255)
                                                                                                               atatCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCAt                                                                                        <  3:972281/96‑1 (MQ=255)
                                                                                                                 atCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCa                                 >  3:260292/1‑149 (MQ=255)
                                                                                                                  tCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCaa                                >  3:218490/1‑149 (MQ=255)
                                                                                                                  tCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCaa                                >  3:800610/1‑149 (MQ=255)
                                                                                                                  tCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCaa                                >  4:456885/1‑149 (MQ=255)
                                                                                                                  tCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCaa                                >  4:755633/1‑149 (MQ=255)
                                                                                                                   cAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCaaa                               <  3:609771/149‑1 (MQ=255)
                                                                                                                    aaTTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAAt                              >  3:459539/1‑149 (MQ=255)
                                                                                                                    aaTTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAAt                              <  4:218490/149‑1 (MQ=255)
                                                                                                                    aaTTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAAt                              >  2:59363/1‑149 (MQ=255)
                                                                                                                     aTTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAAtt                             >  1:1001937/1‑149 (MQ=255)
                                                                                                                      tttCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAAt                                                     <  2:861918/124‑1 (MQ=255)
                                                                                                                      tttCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAAt                                                     >  1:861918/1‑124 (MQ=255)
                                                                                                                      tttCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTa                            >  4:438531/1‑149 (MQ=255)
                                                                                                                         cACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTatca                         >  3:660248/1‑149 (MQ=255)
                                                                                                                           ccGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTatcatc                       <  2:391896/149‑1 (MQ=255)
                                                                                                                           ccGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTatcatc                       <  4:260292/149‑1 (MQ=255)
                                                                                                                            cGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCg                      >  1:1055240/1‑149 (MQ=255)
                                                                                                                            cGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCg                      >  4:345671/1‑149 (MQ=255)
                                                                                                                            cGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCg                      <  2:354599/149‑1 (MQ=255)
                                                                                                                                cTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCaa                  >  3:480635/1‑149 (MQ=255)
                                                                                                                                     ggAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCg             >  4:988468/1‑149 (MQ=255)
                                                                                                                                       aaaGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGtgtg           >  2:796902/1‑149 (MQ=255)
                                                                                                                                        ggattgcgcGAAGATTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTgg          >  3:248454/6‑149 (MQ=255)
                                                                                                                                        aaGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTgg          <  1:844569/149‑1 (MQ=255)
                                                                                                                                           ttgcgcGAAGATTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCgg       >  3:191411/3‑149 (MQ=255)
                                                                                                                                           ttgcgcGAAGATTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCgg       >  1:713897/3‑149 (MQ=255)
                                                                                                                                           ttgcgcGAAGATTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCgg       >  1:439164/3‑149 (MQ=255)
                                                                                                                                           ttgcgcGAAGATTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCAGCAAAGGCGTGGCgg       >  3:438291/3‑149 (MQ=255)
                                                                                                                                                cgAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATgcg  <  3:327987/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCG  >  NC_000913/1177498‑1177790

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