Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 753,309 G→T 5.4% A387E (GCA→GAA)  gltA ← citrate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913753,3090GT 5.4% 272.3 / 3.6 84A387E (GCA→GAA) gltAcitrate synthase
Reads supporting (aligned to +/- strand):  ref base G (34/44);  new base T (2/3);  total (36/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

TATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTACGGGTCGTTCAGCGCGATGTTTTCCAGCTCCATAGCCACTTC  >  NC_000913/753161‑753457
                                                                                                                                                    |                                                                                                                                                    
tatTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTg                                                                                                                                                      >  5:367451/1‑149 (MQ=255)
       aaCTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGc                                                                                                                                               <  6:376269/149‑1 (MQ=255)
             gCAATCAACCATTCACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGaaaa                                                                                                                                          <  5:273638/148‑1 (MQ=255)
             gCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGaaaa                                                                                                                                          <  5:115855/148‑1 (MQ=255)
             gCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAAt                                                                                                                                         <  8:67242/149‑1 (MQ=255)
                 tCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACg                                                                                                                                     <  4:173878/149‑1 (MQ=255)
                 tCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTAATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACg                                                                                                                                     >  1:57320/1‑149 (MQ=255)
                   aaCCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACgg                                                                                                                                    <  5:212182/148‑1 (MQ=255)
                     ccATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTg                                                                                                                                  >  4:292213/1‑148 (MQ=255)
                       aTTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAAc                                                                                                                               <  2:354785/149‑1 (MQ=255)
                           aCGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCGAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATgg                                                                                                                           <  8:316882/149‑1 (MQ=255)
                           aCGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATACTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATg                                                                                                                            >  4:321591/1‑148 (MQ=255)
                             gCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATgg                                                                                                                           >  4:130532/1‑147 (MQ=255)
                             gCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGaa                                                                                                                         <  4:100740/149‑1 (MQ=255)
                               ttGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTTATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGa                                                                                                                       <  2:57320/149‑1 (MQ=255)
                               ttGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACAt                                                                                                                             <  5:269193/143‑1 (MQ=255)
                               ttGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACAt                                                                                                                             >  6:269193/1‑143 (MQ=255)
                                 tatatCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTACCATTGCGAAAATGACGGTGAACATGGAAGACg                                                                                                                     <  1:92347/148‑1 (MQ=255)
                                   tatCGCTTTTAAAGTCGTGTTTTTCATTTCCTGTATACAGCTGACGCGGACGGGGAAAATACATACCGTCACACTGAATTTAGCTACAGTGGGCGATGCCGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGCCgg                                                                                                                    >  6:15604/1‑148 (MQ=255)
                                   tatCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGa                                                                                                                   <  7:279074/149‑1 (MQ=255)
                                     tCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGCCGGaa                                                                                                                  >  3:104637/1‑148 (MQ=255)
                                     tCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCAATTGCGAAAATGACGGTGAACATGGAAGACGGaa                                                                                                                  >  5:202704/1‑148 (MQ=255)
                                        cTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATAcc                                                                                                              >  6:352116/1‑149 (MQ=255)
                                            tAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCAt                                                                                                           <  3:65368/148‑1 (MQ=255)
                                              aaGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCg                                                                                                         <  1:341396/148‑1 (MQ=255)
                                                 tCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGGGATGCAGCGAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCttt                                                                                                     <  2:368552/149‑1 (MQ=255)
                                                   gcgTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCa                                                                                                   <  5:124949/149‑1 (MQ=255)
                                                   gcgTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCa                                                                                                   <  5:352116/149‑1 (MQ=255)
                                                    cgTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAg                                                                                                  <  7:232801/149‑1 (MQ=255)
                                                        tttCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGgat                                                                                               >  8:11665/1‑148 (MQ=255)
                                                        tttAATCTCCTGTTTACAGCTGACGCGGACGGGCAATATTCATACCGTCACTTTGCATTTCGCTCCAGTGGGCTATCCAGCCTACGGTACGTGCCATTGCGAAAATGACGGTGAACATTGAAGACGGAATACCCATCTCTTTCAGgat                                                                                               <  4:104637/148‑1 (MQ=255)
                                                         ttCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGgatg                                                                                              <  8:191970/148‑1 (MQ=255)
                                                             tatCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATa                                                                                           >  3:316444/1‑147 (MQ=255)
                                                              atCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAAGGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGGTTTCAGGATGATAc                                                                                          >  4:287963/1‑147 (MQ=255)
                                                              atCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATAc                                                                                          >  1:259695/1‑147 (MQ=255)
                                                              atCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATAc                                                                                          >  3:276946/1‑147 (MQ=255)
                                                               tCCTTTATACAGCTGACTCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCTGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCa                                                                                        <  6:315120/148‑1 (MQ=255)
                                                               tCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCa                                                                                        <  1:131640/148‑1 (MQ=255)
                                                                   gTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGa                                                                                                                       <  7:18066/113‑1 (MQ=255)
                                                                   gTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGa                                                                                                                       >  8:18066/1‑113 (MQ=255)
                                                                    tataCAGCTGACGCGGAAGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTa                                                                                   <  6:186238/148‑1 (MQ=255)
                                                                    gataCCGCTGACGCGTACGGTCAATCTTCATACCGTCACTGGTCCTTTCGCTCCAGGGGGCGATCCAGCCATCGGTACGGTCCATTGCGAAAATGACGGTGAACATGGAAGACGTAATACCCATCGCTGTGAGGATGATACCAGAGTa                                                                                   <  6:260232/147‑1 (MQ=255)
                                                                      taCAGCTGACGCGGATGGGCAATCTTCATACCTGCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGCAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGaa                                                                                <  4:124776/149‑1 (MQ=255)
                                                                         aGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGTCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATc                                                                             >  8:294180/1‑149 (MQ=255)
                                                                         aGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATc                                                                             >  5:19385/1‑149 (MQ=255)
                                                                              aCGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGAc                                                                          <  3:292213/147‑1 (MQ=255)
                                                                                gcgGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTc                                                                      >  2:73021/1‑149 (MQ=255)
                                                                                gcgGACGGGCAATCTTCATACAGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTc                                                                      >  1:211/1‑149 (MQ=255)
                                                                                  ggACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCgg                                                                    <  4:57399/149‑1 (MQ=255)
                                                                                     cGGGCAATGTTCATACCGTCTCTGTGCATTTCGCTCCTGTGGGCGATCGAGCAAACGGTACGTGCCATTGCGATAATGACGGTGAACATGGTTGACGGATTACCCATCGCTTTCAGGATGATTCCAGAGTAGAAATCGACGTTCGGGTa                                                                 <  4:298963/149‑1 (MQ=255)
                                                                                      gggCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTCCGGGTa                                                                 >  3:344223/1‑148 (MQ=255)
                                                                                        gCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAg                                                              >  7:229765/1‑149 (MQ=255)
                                                                                                cATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTttc                                                          <  6:19385/145‑1 (MQ=255)
                                                                                                  tACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGTCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCg                                                    <  7:294180/149‑1 (MQ=255)
                                                                                                     cGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTTCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGAt                                                  >  3:325845/1‑148 (MQ=255)
                                                                                                     cGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTTCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGAt                                                  >  3:207127/1‑148 (MQ=255)
                                                                                                     cGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATa                                                 <  3:321591/149‑1 (MQ=255)
                                                                                                      gTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATaa                                                <  7:53661/149‑1 (MQ=255)
                                                                                                       tCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATaaa                                               <  4:276946/149‑1 (MQ=255)
                                                                                                         aCTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAg                                              >  4:278735/1‑148 (MQ=255)
                                                                                                           tgtgCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGCCGGAATACCCCTCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTAc                                           >  1:45793/1‑149 (MQ=255)
                                                                                                           tgtgCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTAc                                           >  8:149107/1‑149 (MQ=255)
                                                                                                            gtgCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTACg                                          >  2:32691/1‑149 (MQ=255)
                                                                                                             tgCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTAc                                           <  2:3294/147‑1 (MQ=255)
                                                                                                             tgCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGACACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTAc                                           <  5:377709/147‑1 (MQ=255)
                                                                                                             tgCATTTCGCTCCAGTGGGCGATCCAGCCAACGGGACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTACgg                                         <  4:271639/149‑1 (MQ=255)
                                                                                                               cATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTTCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTACGGGt                                       <  4:325845/149‑1 (MQ=255)
                                                                                                               cATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTTCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTACGGGt                                       <  4:207127/149‑1 (MQ=255)
                                                                                                                 tttCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTACGGGTCg                                     >  2:97471/1‑149 (MQ=255)
                                                                                                                    cGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTACGGGTCGtt                                   <  1:97471/148‑1 (MQ=255)
                                                                                                                    cGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTACGGGTCGtt                                   <  6:202704/148‑1 (MQ=255)
                                                                                                                     gCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTACGGGTCGtt                                   >  2:357248/1‑147 (MQ=255)
                                                                                                                     gCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTACGGGTCGTTCa                                 <  1:357248/149‑1 (MQ=255)
                                                                                                                         cAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTa                                                                                   <  8:394126/95‑1 (MQ=255)
                                                                                                                         cAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTa                                                                                   >  7:394126/1‑95 (MQ=255)
                                                                                                                         cAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTACGGTGCGttacccgc                             >  6:131501/1‑143 (MQ=255)
                                                                                                                         cAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTACGGGTCGTTCAgcgc                             >  1:354605/1‑149 (MQ=255)
                                                                                                                           gTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTACGGGTCGTTCAGCGCGa                           <  2:211/149‑1 (MQ=255)
                                                                                                                            tGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGATCATGGAAGACGGAATACCCATCGCTTTCAGGATGATTCCAGAGTAGATATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTACGGGTCGTTCAGCGCGAt                          <  5:15604/149‑1 (MQ=255)
                                                                                                                            tGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGCTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCAATAAAGTACGTGTCGTTCAGCGCGAt                          >  8:371957/1‑149 (MQ=255)
                                                                                                                              ggCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTACGGGTCGTTCAGCGCGATGt                        >  3:205135/1‑149 (MQ=255)
                                                                                                                                cGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTACGGGTCGTTCAGCGCGATGttt                      <  8:149881/149‑1 (MQ=255)
                                                                                                                                  aTCCAGCCAACGGTACGTGCCATTGCGAAACTGACGGTGATCATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGTGTACCGTTTCTTGTCGATAAAGTACGGGTCGTTAAGCGCTATGTTTTc                    >  8:160824/1‑149 (MQ=255)
                                                                                                                                         cAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTACGGGTCGTTCAGCGCGATGTTTTCCAGCTcc             <  3:130532/149‑1 (MQ=255)
                                                                                                                                                  ttGCCATTTCGAAAATTACGTTGAACATCGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGATTAGAAATCGAATTTCGTGTACACTTTCTTCTCGATAAATTAAGCGTCGTTAAGCTCGATGTTATCCAGATCCATAGCCaat    >  8:171777/2‑147 (MQ=255)
                                                                                                                                                    gCCCTTGCGAAAATGACGGTGAACATGGAAGCCGGAAGACCCCTCGCGTTCAGGATGATACCAGAGTAGCAATCGACGGTCGGGTACAGTTGCTTCTCGATAAAGTACGGGTCGTTCAGCGCGATGTTTTCCAGCTCCATAGCCGCTTc  >  4:292872/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTACGGGTCGTTCAGCGCGATGTTTTCCAGCTCCATAGCCACTTC  >  NC_000913/753161‑753457

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