Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,145,313 C→A 100% S683Y (TCC→TAC)  yegE → putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,145,3130CA98.1% 579.9 / 1.8 158S683Y (TCC→TAC) yegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base A (69/86);  total (71/87)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 160 total aligned reads displayed.

ACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCC  >  NC_000913/2145169‑2145458
                                                                                                                                                |                                                                                                                                                 
aCTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCat                                                                                                                                               <  2:1243484/149‑1 (MQ=255)
aCTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCat                                                                                                                                               >  3:11923/1‑149 (MQ=255)
  tGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCatga                                                                                                                                             <  3:794612/149‑1 (MQ=255)
    gcGGAGCGGCGGCTGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCatgatg                                                                                                                                           >  4:1091206/2‑149 (MQ=255)
       gAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCAc                                                                                                                                        <  4:1424871/149‑1 (MQ=255)
            gcggcAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACg                                                                                                                                   >  1:173097/1‑149 (MQ=255)
             cggcAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGc                                                                                                                                  <  2:1249991/149‑1 (MQ=255)
             cggcAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGc                                                                                                                                  <  3:356042/149‑1 (MQ=255)
                   ttACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCg                                                                                                                            >  2:881840/1‑149 (MQ=255)
                    tACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGc                                                                                                                           <  2:734124/149‑1 (MQ=255)
                    tACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGc                                                                                                                           <  3:37688/149‑1 (MQ=255)
                      cgacgTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGc                                                                                                                                  >  3:1076440/1‑140 (MQ=255)
                       gacgTGCATTACAGTATTACGCCGTTACGTACTCTGGCCGGCAGTAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCTCGCAAACTGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCAGCTCGCCaa                                                                                                                        >  4:1037866/1‑149 (MQ=255)
                       gacgTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGATCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCaa                                                                                                                        >  4:1163893/1‑149 (MQ=255)
                          gTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCatg                                                                                                                                              >  4:68798/1‑124 (MQ=255)
                          gTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAAAGCTGCACCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAAcca                                                                                                                     >  2:320072/1‑148 (MQ=255)
                               ttACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCa                                                                                                                >  4:1400578/1‑149 (MQ=255)
                               ttACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCa                                                                                                                >  4:244136/1‑149 (MQ=255)
                                     tATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTg                                                                                                          >  3:201409/1‑149 (MQ=255)
                                      aTTATGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTga                                                                                                         <  3:1433745/149‑1 (MQ=255)
                                      aTTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTga                                                                                                         <  3:1071131/149‑1 (MQ=255)
                                      aTTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTga                                                                                                         <  4:201767/149‑1 (MQ=255)
                                      aTTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTga                                                                                                         <  1:1143085/149‑1 (MQ=255)
                                      aTTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTga                                                                                                         <  2:1527631/149‑1 (MQ=255)
                                      aTTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTga                                                                                                         <  1:659570/149‑1 (MQ=255)
                                          cGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGaaa                                                                                                     >  4:670792/1‑149 (MQ=255)
                                              gTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAaacaac                                                                                                 <  2:414556/149‑1 (MQ=255)
                                              gTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAaacaac                                                                                                 <  4:11923/149‑1 (MQ=255)
                                               ttAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACt                                                                                                >  4:1458387/1‑149 (MQ=255)
                                               ttAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACt                                                                                                >  2:824161/1‑149 (MQ=255)
                                                tAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTg                                                                                               <  4:133177/149‑1 (MQ=255)
                                                    tACTCTGGACGGCAGCAATATTGGTTCGGCTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTacc                                                                                                                                                 <  1:1302098/95‑1 (MQ=255)
                                                    tACTCTGGACGGCAGCAATATTGGTTCGGCTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTacc                                                                                                                                                 >  2:1302098/1‑95 (MQ=255)
                                                      ctctGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATc                                                                                         >  2:179199/1‑149 (MQ=255)
                                                       tctGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCatgat                                                                                                                                            >  3:397998/1‑97 (MQ=255)
                                                       tctGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATcc                                                                                        >  3:52026/1‑149 (MQ=255)
                                                       tctGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATcc                                                                                        >  3:431750/1‑149 (MQ=255)
                                                       tctGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATcc                                                                                        >  2:1366990/1‑149 (MQ=255)
                                                       tctGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATcc                                                                                        >  4:1051887/1‑149 (MQ=255)
                                                          ggACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCctgc                                                                                     >  3:745544/1‑149 (MQ=255)
                                                                 agcaATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAAc                                                                              >  4:805960/1‑149 (MQ=255)
                                                                       atTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTaaa                                                                        <  3:96344/149‑1 (MQ=255)
                                                                       atTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTaaa                                                                        >  3:1276988/1‑149 (MQ=255)
                                                                         tGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACa                                                                      <  1:1366990/149‑1 (MQ=255)
                                                                          ggttcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAg                                                                     <  4:47511/149‑1 (MQ=255)
                                                                          ggttcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAg                                                                     <  1:181939/149‑1 (MQ=255)
                                                                            ttcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTa                                                                   >  3:1014049/1‑149 (MQ=255)
                                                                            ttcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAAAGGTAAACAGTa                                                                   <  3:1458387/149‑1 (MQ=255)
                                                                             tcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTac                                                                  <  4:380617/149‑1 (MQ=255)
                                                                              cggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCa                                                                                                                                                <  4:1575053/70‑1 (MQ=255)
                                                                              cggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCa                                                                                                                                                >  3:1575053/1‑70 (MQ=255)
                                                                                 ttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACAt                                                              >  1:1224722/1‑149 (MQ=255)
                                                                                     ggTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGc                                                          >  1:466604/1‑149 (MQ=255)
                                                                                     ggTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGc                                                          <  4:1489612/149‑1 (MQ=255)
                                                                                      gTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGGCGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCg                                                         >  4:446060/1‑149 (MQ=255)
                                                                                      gTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCg                                                         >  2:999419/1‑149 (MQ=255)
                                                                                        gATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGAc                                                       >  2:514090/1‑149 (MQ=255)
                                                                                         aTTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACa                                                      >  4:321797/1‑149 (MQ=255)
                                                                                           tCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCaaa                                                                               <  4:790217/122‑1 (MQ=255)
                                                                                           tCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCaaa                                                                               >  3:790217/1‑122 (MQ=255)
                                                                                            cAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACcgc                                                                                                                    >  1:494746/1‑84 (MQ=255)
                                                                                            cAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGc                                                   <  2:719158/149‑1 (MQ=255)
                                                                                            cAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGc                                                   <  3:1400578/149‑1 (MQ=255)
                                                                                             aGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGcc                                                  >  4:1032107/1‑149 (MQ=255)
                                                                                              ggACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGccc                                                 >  3:837753/1‑149 (MQ=255)
                                                                                               gACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCt                                                <  1:100809/149‑1 (MQ=255)
                                                                                               gACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCt                                                >  2:702545/1‑149 (MQ=255)
                                                                                                      ccGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGc                                                   >  1:435166/1‑139 (MQ=255)
                                                                                                      ccGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGc                                                   <  2:435166/139‑1 (MQ=255)
                                                                                                      ccGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGttt                                         >  3:161204/1‑149 (MQ=255)
                                                                                                       cGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTa                                        <  1:824161/149‑1 (MQ=255)
                                                                                                       cGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTa                                        <  4:314389/149‑1 (MQ=255)
                                                                                                       cGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTa                                        >  1:1611908/1‑149 (MQ=255)
                                                                                                       cGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTa                                        >  1:1478262/1‑149 (MQ=255)
                                                                                                        gAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTAt                                       <  2:793621/149‑1 (MQ=255)
                                                                                                         aaTCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATc                                      <  3:244136/149‑1 (MQ=255)
                                                                                                         aaTCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATc                                      <  3:1163893/149‑1 (MQ=255)
                                                                                                            cACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGAt                                   <  3:1043161/149‑1 (MQ=255)
                                                                                                            cACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGAt                                   <  2:1028436/149‑1 (MQ=255)
                                                                                                             aCGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATc                                  >  1:1523764/1‑149 (MQ=255)
                                                                                                             aCGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATc                                  <  3:321797/149‑1 (MQ=255)
                                                                                                              cGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCt                                 <  2:1458441/149‑1 (MQ=255)
                                                                                                                 aaaaTGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCa                                                            <  2:564520/119‑1 (MQ=255)
                                                                                                                  aaaTGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTaaa                                                                        >  2:2868/1‑106 (MQ=255)
                                                                                                                  aaaTGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTaaa                                                                        <  1:2868/106‑1 (MQ=255)
                                                                                                                  aaaTGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGaacatggt                             <  4:1184297/149‑8 (MQ=255)
                                                                                                                     tgctgcGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGc                          <  2:466604/149‑1 (MQ=255)
                                                                                                                       ctgcGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCtt                        <  1:1676569/149‑1 (MQ=255)
                                                                                                                         gcgcCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTa                      >  4:1228673/1‑149 (MQ=255)
                                                                                                                              aGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCg                 <  4:418540/149‑1 (MQ=255)
                                                                                                                              aGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCg                 <  3:949067/149‑1 (MQ=255)
                                                                                                                                cTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGt               <  2:1478262/149‑1 (MQ=255)
                                                                                                                                 tGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTg              >  3:271489/1‑149 (MQ=255)
                                                                                                                                 tGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTg              >  4:784118/1‑149 (MQ=255)
                                                                                                                                    gCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAAt           <  3:670792/149‑1 (MQ=255)
                                                                                                                                     cTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATg          >  4:141592/1‑149 (MQ=255)
                                                                                                                                       aCAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATCCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGAc        >  2:321568/1‑149 (MQ=255)
                                                                                                                                        cAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACa       <  3:742834/149‑1 (MQ=255)
                                                                                                                                        cAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACa       <  3:5093/149‑1 (MQ=255)
                                                                                                                                        cAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACa       <  3:1228673/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
ACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCC  >  NC_000913/2145169‑2145458

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: