Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,870,261 G→A 5.8% intergenic (+23/‑124) yeaH → / → yeaI UPF0229 family protein/putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,870,2610GA 5.8% 331.3 / 6.4 104intergenic (+23/‑124)yeaH/yeaIUPF0229 family protein/putative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (58/39);  new base A (2/4);  total (60/44)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.46e-01

Only 100 of 104 total aligned reads displayed.

AGAATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTC  >  NC_000913/1870115‑1870407
                                                                                                                                                  |                                                                                                                                                  
 gaATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGcc                                                                                                                                                  <  7:129766/148‑1 (MQ=255)
  aaTATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCt                                                                                                                                                 <  6:16739/148‑1 (MQ=255)
   atatGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAA‑‑ATCAGCCAGGCGATAATGACCTGGc                                                                                                                                              <  6:52882/148‑1 (MQ=255)
     atGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGc                                                                                                                                              >  8:284794/1‑148 (MQ=255)
      tGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACGGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCACCTGGCt                                                                                                                                             >  7:399766/1‑148 (MQ=255)
       gAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCt                                                                                                                                             >  7:116659/1‑147 (MQ=255)
              tGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGcc                                                                                                                                                  >  8:235314/1‑135 (MQ=255)
              tGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGcc                                                                                                                                                  <  7:235314/135‑1 (MQ=255)
               gCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAg                                                                                                                                    >  1:445602/1‑148 (MQ=255)
               gCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAg                                                                                                                                    >  8:168328/1‑148 (MQ=255)
               gCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGAAATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGACATTATCGCATGGCTGATTTTTAg                                                                                                                                    >  2:270179/1‑148 (MQ=255)
                cAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGc                                                                                                                                   >  4:391230/1‑148 (MQ=255)
                  aTCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCt                                                                                                                                  >  3:93560/1‑147 (MQ=255)
                   tCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGAttt                                                                                                                                        >  2:84785/1‑140 (MQ=255)
                   tCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGAttt                                                                                                                                        <  1:84785/140‑1 (MQ=255)
                   tCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTa                                                                                                                                     >  2:81001/1‑143 (MQ=255)
                   tCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTa                                                                                                                                     <  1:81001/143‑1 (MQ=255)
                    cTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTAc                                                                                                                               >  7:28307/1‑148 (MQ=255)
                       cTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTAc                                                                                                                               <  7:284794/145‑1 (MQ=255)
                         ttCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTa                                                                                                                                     >  8:318709/1‑137 (MQ=255)
                         ttCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTa                                                                                                                                     <  7:318709/137‑1 (MQ=255)
                         ttCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTa                                                                                                                           >  7:105415/1‑147 (MQ=255)
                             aCAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCAGGTGTTCCGTGAACTGATTCATAAACAAAATGCAAAAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTa                                                                                                                      >  5:387269/1‑148 (MQ=255)
                              cAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCtt                                                                                                                                 >  5:308728/1‑136 (MQ=255)
                              cAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCtt                                                                                                                                 <  6:308728/136‑1 (MQ=255)
                                  tttGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCa                                                                                                                 >  3:192922/1‑148 (MQ=255)
                                   ttGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCAt                                                                                                                >  3:208535/1‑148 (MQ=255)
                                    tGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACGGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCACCTGGCTGATTTTTAGCTTACTGTAAATTATCTCAt                                                                                                                <  8:399766/147‑1 (MQ=255)
                                       gATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCACCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATttat                                                                                                            <  1:336767/148‑1 (MQ=255)
                                            agcaCATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATtacat                                                                                                       <  4:192922/148‑1 (MQ=255)
                                               acaTCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATtacataca                                                                                                    >  5:27435/1‑148 (MQ=255)
                                                caTCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTacatacat                                                                                                   <  8:105415/148‑1 (MQ=255)
                                                caTCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTacatacat                                                                                                   <  7:701/148‑1 (MQ=255)
                                                 aTCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACAtt                                                                                                  <  6:166133/148‑1 (MQ=255)
                                                   ccGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTaa                                                                                                                          >  1:378216/1‑122 (MQ=255)
                                                   ccGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTaa                                                                                                                          <  2:378216/122‑1 (MQ=255)
                                                      cgACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCtt                                                                                             <  6:27435/148‑1 (MQ=255)
                                                           aGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAAt                                                                                        >  2:261614/1‑148 (MQ=255)
                                                               tgatATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCt                                                                                    >  6:380027/1‑148 (MQ=255)
                                                                        tCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAatat                                                                            >  8:210958/1‑147 (MQ=255)
                                                                        tCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATg                                                                           <  3:391230/148‑1 (MQ=255)
                                                                          cGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGAc                                                                         <  8:116659/148‑1 (MQ=255)
                                                                            gtgtTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAg                                                                       <  8:28307/148‑1 (MQ=255)
                                                                                tCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCtt                                                                                             >  4:97816/1‑122 (MQ=255)
                                                                                tCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCtt                                                                                             <  3:97816/122‑1 (MQ=255)
                                                                                    tGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCt                                                               >  5:330494/1‑148 (MQ=255)
                                                                                          gTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAAtt                                                                                                                       >  5:112224/1‑86 (MQ=255)
                                                                                          gTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAAtt                                                                                                                       <  6:112224/86‑1 (MQ=255)
                                                                                          gTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACa                                                          >  4:130100/1‑147 (MQ=255)
                                                                                               aTAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAAttt                                                    <  5:380027/148‑1 (MQ=255)
                                                                                                  aaCAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAg                                                  <  4:93560/147‑1 (MQ=255)
                                                                                                   aCAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGtt                                                >  3:310441/1‑148 (MQ=255)
                                                                                                    cAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTc                                               >  6:37227/1‑148 (MQ=255)
                                                                                                    cAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTc                                               >  6:114742/1‑148 (MQ=255)
                                                                                                       aaTGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTc                                               >  3:172884/1‑145 (MQ=255)
                                                                                                       aaTGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAg                                            <  5:37227/148‑1 (MQ=255)
                                                                                                             aCAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATTACAGCATAACATTTACATAATTTAGTTCCAGAAAAc                                       >  4:67635/1‑147 (MQ=255)
                                                                                                             aCAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATaa                                                                  <  3:214875/120‑1 (MQ=255)
                                                                                                             aCAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATaa                                                                  >  4:214875/1‑120 (MQ=255)
                                                                                                             aCAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTa                                                   <  4:39601/135‑1 (MQ=255)
                                                                                                             aCAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTa                                                   >  3:39601/1‑135 (MQ=255)
                                                                                                                    aaGGCTTAATCTATCAGCCATGTCATTATCTCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATTTGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCAtt                                <  1:270179/147‑1 (MQ=255)
                                                                                                                    aaGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCAtt                                <  2:445602/147‑1 (MQ=255)
                                                                                                                     aGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCg                              >  4:251776/1‑148 (MQ=255)
                                                                                                                     aGGCTAAAACTATCAGACAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCg                              >  7:130221/1‑148 (MQ=255)
                                                                                                                          aaaaCTATCAGCCAGGTCATTATCGCCTTGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAat                                                                              <  7:313188/95‑1 (MQ=255)
                                                                                                                          aaaaCTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAat                                                                              >  8:313188/1‑95 (MQ=255)
                                                                                                                          aaaaCTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTc                               <  8:211546/142‑1 (MQ=255)
                                                                                                                          aaaaCTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTc                               >  7:211546/1‑142 (MQ=255)
                                                                                                                          aaaaCTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGaaa                          >  4:306332/1‑147 (MQ=255)
                                                                                                                          aaaaCTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGaaa                          >  4:300642/1‑147 (MQ=255)
                                                                                                                          aaaaCTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCGTTCAATATGACAGCATACCCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGaaa                          >  5:362024/1‑147 (MQ=255)
                                                                                                                           aaaCTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGaaaaa                        <  4:208535/148‑1 (MQ=255)
                                                                                                                            aaCTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAAt                       >  8:64946/1‑148 (MQ=255)
                                                                                                                            aaCTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAAt                       >  1:263232/1‑148 (MQ=255)
                                                                                                                              cTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGa                     <  1:261614/148‑1 (MQ=255)
                                                                                                                                aTCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGAtt                   >  6:219730/1‑148 (MQ=255)
                                                                                                                                 tCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCa                                                                     >  3:106957/1‑97 (MQ=255)
                                                                                                                                 tCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCa                                                                     <  4:106957/97‑1 (MQ=255)
                                                                                                                                 tCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGAt                    <  7:64946/146‑1 (MQ=255)
                                                                                                                                 tCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTc                  <  3:130100/148‑1 (MQ=255)
                                                                                                                                  cAGCCAGGCGATAATGACCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCtt                                                              >  4:57219/1‑103 (MQ=38)
                                                                                                                                  cAGCCAGGCGATAATGACCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCtt                                                              <  3:57219/103‑1 (MQ=38)
                                                                                                                                       aGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCaa            <  4:310441/148‑1 (MQ=255)
                                                                                                                                        ggTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTATATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTGCATAATTTAGTTCCAGAAAACAAACAGTCGGAAAAATGATTCAGTCaa            >  2:323483/1‑147 (MQ=255)
                                                                                                                                        ggTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAg                                            <  5:309378/115‑1 (MQ=255)
                                                                                                                                        ggTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAg                                            >  6:309378/1‑115 (MQ=255)
                                                                                                                                         gTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGaaa                                         <  4:288835/117‑1 (MQ=255)
                                                                                                                                         gTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGaaa                                         >  3:288835/1‑117 (MQ=255)
                                                                                                                                         gTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAaca          >  5:161703/1‑148 (MQ=255)
                                                                                                                                          tCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATtacataca                                                                                                    >  4:38306/1‑57 (MQ=255)
                                                                                                                                          tCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATaa                                                                  <  7:32313/91‑1 (MQ=255)
                                                                                                                                          tCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATaa                                                                  >  8:32313/1‑91 (MQ=255)
                                                                                                                                          tCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAacac         >  4:275830/1‑148 (MQ=255)
                                                                                                                                          tCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAacac         >  6:219954/1‑148 (MQ=255)
                                                                                                                                          tCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAacac         >  6:25162/1‑148 (MQ=255)
                                                                                                                                          tCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTAAATAATTTAGTACCAGAAAACAATCATTCGGAAAAATGATTCAGTCAacac         >  2:33774/1‑148 (MQ=255)
                                                                                                                                               aTCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGt               >  8:215312/1‑137 (MQ=255)
                                                                                                                                               aTCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGt               <  7:215312/137‑1 (MQ=255)
                                                                                                                                                 cGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTc  >  2:75556/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AGAATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTC  >  NC_000913/1870115‑1870407

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