Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,145,313 C→A 53.5% 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,3130CA53.5% ‑0.9 / 147.2 86S683Y (TCC→TAC) yegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (20/20);  new base A (13/33);  total (33/53)
Fisher's exact test for biased strand distribution p-value = 4.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

TTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGG  >  NC_000913/2145165‑2145461
                                                                                                                                                    |                                                                                                                                                    
ttGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTc                                                                                                                                                      >  3:466827/1‑149 (MQ=255)
   cACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTccc                                                                                                                                                    <  2:139154/148‑1 (MQ=255)
   cACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCa                                                                                                                                                   <  4:767696/149‑1 (MQ=255)
    aCTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCat                                                                                                                                                  >  3:948575/1‑149 (MQ=255)
    aCTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCa                                                                                                                                                   >  2:253236/1‑148 (MQ=255)
     cTGTCGGAGCGGCGGCAGTTATCACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGAGATATTGGTTCGGATCTGGTGATTCAGGATGTTACCGAATCACTCAAAATGAGTCGCCAACTGAGCTACAGTTCCTCCCatg                                                                                                                                                 >  4:68409/1‑149 (MQ=255)
      tGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCatg                                                                                                                                                 <  1:253236/148‑1 (MQ=255)
      tGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCatg                                                                                                                                                 <  2:1979/148‑1 (MQ=255)
           gAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCAc                                                                                                                                           <  4:466827/149‑1 (MQ=255)
                 cggcAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACg                                                                                                                                      >  2:841212/1‑148 (MQ=255)
                    cAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGAGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCAtc                                                                                                                                  <  3:114164/149‑1 (MQ=255)
                      gTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCAtct                                                                                                                                 <  1:605124/148‑1 (MQ=255)
                      gTTACGACGTGCATTACAGGATTACACCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCAtctc                                                                                                                                <  4:449345/149‑1 (MQ=255)
                       ttACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCAtctc                                                                                                                                <  2:586446/148‑1 (MQ=255)
                        tACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCg                                                                                                                               >  1:99385/1‑148 (MQ=255)
                        tACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCg                                                                                                                               >  1:847773/1‑148 (MQ=255)
                             cgTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAACATTGGTTCAGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAAcc                                                                                                                         <  3:823711/149‑1 (MQ=255)
                                  aTTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGcc                                                                                                                    >  4:883056/1‑149 (MQ=255)
                                    tACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCa                                                                                                                   >  2:747888/1‑148 (MQ=255)
                                    tACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCTAACCGCGCCAg                                                                                                                  <  3:762927/149‑1 (MQ=255)
                                          aTTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTg                                                                                                             <  2:167688/148‑1 (MQ=255)
                                          aTTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTg                                                                                                             <  2:929207/148‑1 (MQ=255)
                                          aTTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTg                                                                                                             <  2:380001/148‑1 (MQ=255)
                                                   ttAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACt                                                                                                   >  4:265625/1‑149 (MQ=255)
                                                   ttAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAaacaac                                                                                                    >  1:321790/1‑148 (MQ=255)
                                                   ttAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACt                                                                                                   <  3:335958/149‑1 (MQ=255)
                                                    tAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACt                                                                                                   >  2:213518/1‑148 (MQ=255)
                                                     aaGTACTCTGGACGGCAGCAATATTTGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGc                                                                                                 <  3:673994/149‑1 (MQ=255)
                                                     aaGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTg                                                                                                  <  2:454028/148‑1 (MQ=255)
                                                     aaGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGc                                                                                                 <  3:588803/149‑1 (MQ=255)
                                                     aaGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGc                                                                                                 <  3:645009/149‑1 (MQ=255)
                                                      aGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGc                                                                                                 >  2:587074/1‑148 (MQ=255)
                                                           tctGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATcc                                                                                           >  3:132718/1‑149 (MQ=255)
                                                               gACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCctgct                                                                                       <  3:883056/149‑1 (MQ=255)
                                                               gACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCctgc                                                                                        <  2:534404/148‑1 (MQ=255)
                                                                     agcaATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAAc                                                                                 >  4:839992/1‑149 (MQ=255)
                                                                      gcaATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAAc                                                                                 >  2:3455/1‑148 (MQ=255)
                                                                       caATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACgg                                                                               <  4:241718/149‑1 (MQ=255)
                                                                        aaTATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACgg                                                                               <  2:399579/148‑1 (MQ=255)
                                                                         atatTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGt                                                                              <  1:756072/148‑1 (MQ=255)
                                                                         atatTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTa                                                                             >  3:606595/1‑149 (MQ=255)
                                                                           atTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTaaa                                                                           <  4:923511/149‑1 (MQ=255)
                                                                           atTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATc                                                                                            <  1:115437/132‑1 (MQ=255)
                                                                           atTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATc                                                                                            >  2:115437/1‑132 (MQ=255)
                                                                               gttcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGt                                                                       >  3:746836/1‑149 (MQ=255)
                                                                                ttcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGt                                                                       <  2:952657/148‑1 (MQ=255)
                                                                                ttcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGt                                                                       <  2:154703/148‑1 (MQ=255)
                                                                                ttcggttcTGGTAATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGt                                                                       >  1:454842/1‑148 (MQ=255)
                                                                                 tcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTac                                                                     <  4:137879/149‑1 (MQ=255)
                                                                                 tcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTa                                                                      >  1:46682/1‑148 (MQ=255)
                                                                                         ggTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAg                                                              <  1:328820/148‑1 (MQ=255)
                                                                                          gTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGGTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGc                                                             >  2:461136/1‑148 (MQ=255)
                                                                                             aTTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGAc                                                          >  2:392757/1‑148 (MQ=255)
                                                                                               tCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATg                                                       >  4:868386/1‑149 (MQ=255)
                                                                                                cAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATg                                                       <  2:99385/148‑1 (MQ=255)
                                                                                                     cGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTaa                                                                            >  2:367715/1‑122 (MQ=255)
                                                                                                     cGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTaa                                                                            <  1:367715/122‑1 (MQ=255)
                                                                                                      gTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTgg                                                 <  2:465641/148‑1 (MQ=255)
                                                                                                          ccGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGttt                                            >  3:893061/1‑149 (MQ=255)
                                                                                                           cGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTa                                           >  3:31473/1‑149 (MQ=255)
                                                                                                           cGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTa                                           >  3:57525/1‑149 (MQ=255)
                                                                                                            gAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTAt                                          <  4:90999/149‑1 (MQ=255)
                                                                                                            gAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTa                                           <  1:454720/148‑1 (MQ=255)
                                                                                                             aaTCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATc                                         <  4:150531/149‑1 (MQ=255)
                                                                                                             aaTCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCCCCCATGATGCACTGCCGCATCTCGCCAACCGCGCCAGTTTTGAGAACCAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATc                                         <  3:68409/149‑1 (MQ=255)
                                                                                                                cACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGa                                       >  1:609490/1‑148 (MQ=255)
                                                                                                                cACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGa                                       <  1:138480/148‑1 (MQ=255)
                                                                                                                  cGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATc                                     >  2:369010/1‑148 (MQ=255)
                                                                                                                  cGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTg                                                                                                             <  1:726813/76‑1 (MQ=255)
                                                                                                                  cGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTg                                                                                                             >  2:726813/1‑76 (MQ=255)
                                                                                                                  cGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCt                                    <  3:624273/149‑1 (MQ=255)
                                                                                                                   gCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCtt                                   >  3:854007/1‑149 (MQ=255)
                                                                                                                   gCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCtt                                   <  4:653748/149‑1 (MQ=255)
                                                                                                                          gctgcGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGc                             <  2:800855/148‑1 (MQ=255)
                                                                                                                           ctgcGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCt                            >  1:924642/1‑148 (MQ=255)
                                                                                                                           ctgcGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCt                            <  2:674189/148‑1 (MQ=255)
                                                                                                                                ccAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTaaa                       <  2:361895/148‑1 (MQ=255)
                                                                                                                                  aGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGc                     <  2:698963/148‑1 (MQ=255)
                                                                                                                                   gCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCgg                   <  3:712283/149‑1 (MQ=255)
                                                                                                                                   gCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCg                    <  2:378481/148‑1 (MQ=255)
                                                                                                                                    cTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGt                  <  4:72835/149‑1 (MQ=255)
                                                                                                                                      gAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTg                 <  2:609490/148‑1 (MQ=255)
                                                                                                                                       aGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGa                <  1:960203/148‑1 (MQ=255)
                                                                                                                                        gCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGaa               <  1:369010/148‑1 (MQ=255)
                                                                                                                                        gCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAAt              <  4:132718/149‑1 (MQ=255)
                                                                                                                                            cAGCGCCTACCATGATGCACTGACGCATCTCGCCACCCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACAACAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACa          <  4:208248/149‑1 (MQ=255)
                                                                                                                                            cAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGAc           <  1:527034/148‑1 (MQ=255)
                                                                                                                                                gcCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAgcg       >  1:595060/1‑148 (MQ=255)
                                                                                                                                                    cccATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCggg  >  4:815100/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGG  >  NC_000913/2145165‑2145461

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