Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,874,050 C→A 5.2% C142F (TGC→TTC) ‡ dgoK ← 2‑oxo‑3‑deoxygalactonate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,874,0500CA 5.2% 307.8 / 5.7 97C142F (TGC→TTC) ‡dgoK2‑oxo‑3‑deoxygalactonate kinase
Reads supporting (aligned to +/- strand):  ref base C (40/52);  new base A (2/3);  total (42/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.96e-01

AAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAACTGACGT  >  NC_000913/3873902‑3874198
                                                                                                                                                    |                                                                                                                                                    
aaGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAGCCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGc                                                                                                                                                      >  1:556628/1‑149 (MQ=255)
 aGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTGACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCa                                                                                                                                                     >  2:1157286/1‑149 (MQ=255)
    ccAGCTGTGAAGTCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTTATTTAACAGTAACTGATGTAATTCACCGGTCATCACGGTGCGAAAATGGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATg                                                                                                                                                  <  1:382369/149‑1 (MQ=255)
    ccAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATg                                                                                                                                                  <  2:1331494/149‑1 (MQ=255)
     cAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATgg                                                                                                                                                 <  2:1199540/149‑1 (MQ=255)
      aGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATggg                                                                                                                                                <  1:1291769/149‑1 (MQ=255)
       gCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAATTGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTGCAAATGGGt                                                                                                                                               <  1:315671/149‑1 (MQ=255)
           tgAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTGCAAATGGGtgcag                                                                                                                                           <  1:980286/149‑5 (MQ=255)
           ggAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTGCAAATGGGtgcag                                                                                                                                           <  1:975156/148‑5 (MQ=255)
                gCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCAt                                                                                                                                      <  2:952730/149‑1 (MQ=255)
                gCATCGGCAGAGTTTTCCTGCGGCGGCAAACCAGCTCAAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCac                                                                                                                                      >  1:77359/1‑148 (MQ=255)
                  aTCGGCAGAGTTTTCCTGTGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGa                                                                                                                                    >  2:1179542/1‑149 (MQ=255)
                  aTCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGa                                                                                                                                    >  2:1138559/1‑149 (MQ=255)
                   tCGGCAGAGTGTAGCTGCGGCGGCAAACCTGCGCCAAGCAATGAGTGATATAACAGTAACTGATGTAATGCACCGGTCATCACGGTGCGAACATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGAc                                                                                                                                   <  1:262819/149‑1 (MQ=255)
                    cGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACa                                                                                                                                  <  1:732618/149‑1 (MQ=255)
                        agagTTTTCCTGCGGCGGCAAGCCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATaaa                                                                                                                              <  2:556628/149‑1 (MQ=255)
                        agagTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATaaa                                                                                                                              <  1:65398/149‑1 (MQ=255)
                        agagTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATaaa                                                                                                                              <  2:1414144/149‑1 (MQ=255)
                          agTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGAGGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAaga                                                                                                                            <  1:113775/149‑1 (MQ=255)
                           gTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAagag                                                                                                                           <  1:989246/149‑1 (MQ=255)
                            ttttCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAACTGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAagaga                                                                                                                          <  1:154399/149‑1 (MQ=255)
                            ttttCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAagaga                                                                                                                          >  1:47760/1‑149 (MQ=255)
                            ttttCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAagaga                                                                                                                          <  2:1292422/149‑1 (MQ=255)
                            ttttCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAagaga                                                                                                                          <  2:101485/149‑1 (MQ=255)
                               tCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCGTGACATAAAGAGagga                                                                                                                       >  1:401295/1‑149 (MQ=255)
                               tCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGagga                                                                                                                       <  1:779976/149‑1 (MQ=255)
                                 cTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGaggaag                                                                                                                     <  1:192603/149‑1 (MQ=255)
                                 cTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGaggaag                                                                                                                     >  1:49709/1‑149 (MQ=255)
                                 cTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGaggaag                                                                                                                     <  2:1199685/149‑1 (MQ=255)
                                       cggcAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCa                                                                                                               <  1:1252798/149‑1 (MQ=255)
                                        ggcAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCag                                                                                                              <  1:1244398/149‑1 (MQ=255)
                                        ggcAAACCTGCGCCAATCAATGAGAGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGATGATTTGCTGGCTATCGGCCTGCACACATTTGCAATGCGTTACGGGCATGACATAAAGAGAGGAAGGAGCCag                                                                                                              >  2:1209977/1‑149 (MQ=255)
                                         gcAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCaga                                                                                                             <  1:1224116/149‑1 (MQ=255)
                                         gcAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCaga                                                                                                             >  2:421721/1‑149 (MQ=255)
                                           aaaCCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGc                                                                                                           <  2:692058/149‑1 (MQ=255)
                                             aCCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTc                                                                                                         <  1:661484/149‑1 (MQ=255)
                                                tGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTcgcg                                                                                                      >  1:117248/1‑149 (MQ=255)
                                                 gcgcCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTcgcgc                                                                                                     >  1:188601/1‑149 (MQ=255)
                                                  cgcCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTcgcgcg                                                                                                    >  2:824108/1‑149 (MQ=255)
                                                    ccAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGcc                                                                                                  <  1:1015102/149‑1 (MQ=255)
                                                    ccAATCAATGAGGGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGcc                                                                                                  <  2:261746/149‑1 (MQ=255)
                                                          aatGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCaa                                                                                            <  1:1021461/149‑1 (MQ=255)
                                                          aatGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCaa                                                                                            >  1:37819/1‑149 (MQ=255)
                                                          aatGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCaa                                                                                            >  2:717478/1‑149 (MQ=255)
                                                           atGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAAt                                                                                           <  1:910925/149‑1 (MQ=255)
                                                            tGAGGGATTTAACAGTAAATGTTGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAAtt                                                                                          <  1:732182/149‑1 (MQ=255)
                                                                tGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTgtgt                                                                                      <  2:343891/149‑1 (MQ=255)
                                                                 gATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGtt                                                                                     >  2:1043567/1‑149 (MQ=255)
                                                                  aTTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGttt                                                                                    >  1:1085095/1‑149 (MQ=255)
                                                                  aTTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGttt                                                                                    >  2:751679/1‑149 (MQ=255)
                                                                      aaCAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGTCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTttcttc                                                                                >  2:1525466/1‑149 (MQ=255)
                                                                      aaCAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTttcttc                                                                                >  2:148148/1‑149 (MQ=255)
                                                                       aCAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTcg                                                                               >  2:1241117/1‑149 (MQ=255)
                                                                        cAGTAAATGATTTAATTCACCGTTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTcgc                                                                              <  2:1220495/149‑1 (MQ=255)
                                                                          gTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTcgccg                                                                            >  1:637187/1‑149 (MQ=255)
                                                                          gTAAATGATGTAATTCACCGGGCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTcgccg                                                                            <  1:932137/149‑1 (MQ=255)
                                                                             aaTGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCgcgc                                                                         >  2:1234112/1‑149 (MQ=255)
                                                                              atgatgTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCa                                                                        >  2:270473/1‑149 (MQ=255)
                                                                               tgatgTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGAGCGCCGATCAATTGTGTATCTTCGCCGCGCAt                                                                       >  1:779268/1‑149 (MQ=255)
                                                                                  tgTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATcac                                                                    >  2:193141/1‑149 (MQ=255)
                                                                                       ttCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATtgt                                                               >  1:321625/1‑149 (MQ=255)
                                                                                       ttCACCGGTCATCACGGTGCGAAAAGCGGTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATtgt                                                               <  2:300587/149‑1 (MQ=255)
                                                                                         cACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCGATCGGCCTGCACCCATTTGCAATGGGGTCCGGGCATGACATAACGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTgg                                                             <  2:49782/149‑1 (MQ=255)
                                                                                                 atcaCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTc                                                     >  1:771199/1‑149 (MQ=255)
                                                                                                  tcaCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCa                                                    <  1:270473/149‑1 (MQ=255)
                                                                                                   caCGGTGCGAACATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGTAAGTAGCCAGATCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCAt                                                   <  1:1388354/149‑1 (MQ=255)
                                                                                                       gTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATgaga                                               <  1:770544/149‑1 (MQ=255)
                                                                                                       gTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATgaga                                               <  1:1138559/149‑1 (MQ=255)
                                                                                                       gTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATgaga                                               <  1:181698/149‑1 (MQ=255)
                                                                                                           gAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGGTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGGTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGaca                                           >  2:103184/1‑149 (MQ=255)
                                                                                                            aaaaTCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGacac                                          >  1:660091/1‑149 (MQ=255)
                                                                                                               aTCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATa                                       >  2:658148/1‑149 (MQ=255)
                                                                                                               aTCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATa                                       >  1:842529/1‑149 (MQ=255)
                                                                                                               aTCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATa                                       <  2:304954/149‑1 (MQ=255)
                                                                                                                  gTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATga                                                 >  2:1340203/1‑136 (MQ=255)
                                                                                                                  gTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATga                                                 <  1:1340203/136‑1 (MQ=255)
                                                                                                                     gATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCgg                                 <  2:509664/149‑1 (MQ=255)
                                                                                                                     gATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCgg                                 <  2:47760/149‑1 (MQ=255)
                                                                                                                     gATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCgg                                 <  2:593123/149‑1 (MQ=255)
                                                                                                                     gATATGCTGGCTATCGGCCTGCACCCATATGCAATGGGTTCCCGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATAGTGTTTATTCGCCGCGCATCACATTGAGGTTATCGTCATGAGAGACACATAATCCgg                                 >  2:1257603/1‑149 (MQ=255)
                                                                                                                             ggCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATc                         >  2:622733/1‑149 (MQ=255)
                                                                                                                             ggCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATc                         <  1:148148/149‑1 (MQ=255)
                                                                                                                              gCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATcc                        >  1:1311821/1‑149 (MQ=255)
                                                                                                                              gCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATcc                        <  1:653422/149‑1 (MQ=255)
                                                                                                                               ctttCGGCCTGCACCCATTTTCAATGGGTTCCGGGCATGAAATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCATTTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGTAATAATCCa                       <  1:1209977/146‑1 (MQ=255)
                                                                                                                                tATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAg                      >  1:650038/1‑149 (MQ=255)
                                                                                                                                 aTCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGa                     >  2:523060/1‑149 (MQ=255)
                                                                                                                                   cGGCCTGCACCTATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGata                   <  2:1564544/149‑1 (MQ=255)
                                                                                                                                    ggCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGatat                  >  2:1480138/1‑149 (MQ=255)
                                                                                                                                     gCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATAtt                 >  1:1188512/1‑149 (MQ=255)
                                                                                                                                      ccggaaCCCATTGCAAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGACACACATAATCCGGGAATAATCCAGATAttt                >  1:971541/6‑148 (MQ=255)
                                                                                                                                      ccTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTg                >  1:1457949/1‑149 (MQ=255)
                                                                                                                                       cTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGt               <  2:49709/149‑1 (MQ=255)
                                                                                                                                       cTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGt               <  2:842529/149‑1 (MQ=255)
                                                                                                                                        ttcACCCATTTGCAATGGGGTCCGGGCATGCCATTAAGAGAGGAAGGAGCCAGGGCTCGCGCGCCGATCAATTGTGTTGCTTCGCCGCGCATCACATTGTGGGTTTCGTCAGGAGAGACACATAATCCGGGAATACTCCAGATATTGTc              <  2:889905/147‑1 (MQ=255)
                                                                                                                                        ggaaCCCATTGCAAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTc              >  1:1017252/4‑149 (MQ=255)
                                                                                                                                          cACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGc            <  1:512606/149‑1 (MQ=255)
                                                                                                                                              cATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAAc        <  1:1157286/149‑1 (MQ=255)
                                                                                                                                                    cAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAACTGACGt  <  1:1241117/149‑1 (MQ=255)
                                                                                                                                                    cAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTCTCTTCGCCGCGCATCACATTGTGGTAATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAACTGACGt  <  2:779268/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAACTGACGT  >  NC_000913/3873902‑3874198

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