Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,455,369 C→T 6.5% D81N (GAT→AAT)  fbp ← fructose‑1,6‑bisphosphatase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,455,3690CT 6.5% 189.1 / 2.6 61D81N (GAT→AAT) fbpfructose‑1,6‑bisphosphatase I
Reads supporting (aligned to +/- strand):  ref base C (34/23);  new base T (2/2);  total (36/25)
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.36e-01

GAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGC  >  NC_000913/4455226‑4455515
                                                                                                                                               |                                                                                                                                                  
gaAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTg                                                                                                                                               >  2:127342/1‑149 (MQ=255)
 aaaaTGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGc                                                                                                                                              >  5:1904/1‑149 (MQ=255)
          ccGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGc                                                                                                                                     <  1:127342/149‑1 (MQ=255)
                     gTTAACATCGATGTTGGACGAGCCATCCAGTGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGtt                                                                                                                          <  1:30259/149‑1 (MQ=255)
                        aaCATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGtttt                                                                                                                        <  3:13255/148‑1 (MQ=255)
                           aTCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATCTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCAt                                                                                                                    >  5:184569/1‑149 (MQ=255)
                            tCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCa                                                                                                                             >  2:204031/1‑139 (MQ=255)
                            tCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCa                                                                                                                             <  1:204031/139‑1 (MQ=255)
                              gATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTa                                                                                                                  <  3:192948/148‑1 (MQ=255)
                                     gACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACaa                                                                                                          <  3:69663/149‑1 (MQ=255)
                                                 aGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGTCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAg                                                                                                         >  8:75783/1‑138 (MQ=255)
                                                 aGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAg                                                                                                         <  7:75783/138‑1 (MQ=255)
                                                         cATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGCACCt                                                                                      >  5:34338/1‑149 (MQ=255)
                                                          aTCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTa                                                                                                                  >  5:191943/1‑120 (MQ=255)
                                                          aTCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTa                                                                                                                  <  6:191943/120‑1 (MQ=255)
                                                           tCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCg                                                                                    >  8:40750/1‑149 (MQ=255)
                                                              gCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGccc                                                                                 <  7:69049/149‑1 (MQ=255)
                                                                accacGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTg                                                                               >  4:166430/1‑149 (MQ=255)
                                                                 ccacGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGc                                                                              <  1:32411/149‑1 (MQ=255)
                                                                  cacGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCa                                                                             >  7:13237/1‑149 (MQ=255)
                                                                  cacGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCa                                                                             >  3:192305/1‑149 (MQ=255)
                                                                     gTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGc                                                                                   <  3:119958/140‑1 (MQ=255)
                                                                     gTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGc                                                                                   >  4:119958/1‑140 (MQ=255)
                                                                      tATCTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGtt                                                                         <  6:184569/149‑1 (MQ=255)
                                                                                gtTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCg                                                               >  5:223787/1‑149 (MQ=255)
                                                                                 ttCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGc                                                              >  6:155542/1‑149 (MQ=255)
                                                                                    acaGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTg                                                            <  6:1904/148‑1 (MQ=255)
                                                                                        ccTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATTGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCAcc                                                       >  1:107746/1‑149 (MQ=255)
                                                                                        ccTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATTGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCAcc                                                       >  1:104364/1‑149 (MQ=255)
                                                                                        ccTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCAcc                                                       >  3:139668/1‑149 (MQ=255)
                                                                                        ccTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTGTCATTATCGAACAAGTCGAGTGTCTGCTGCACCTCGCCCGGCACGTGCTCAGCACCGCTGGCAcc                                                       >  4:183956/1‑149 (MQ=255)
                                                                                         cTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCAccc                                                      >  2:249160/1‑149 (MQ=255)
                                                                                          ttCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTg                                                            <  4:176416/142‑1 (MQ=255)
                                                                                          ttCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTg                                                            >  3:176416/1‑142 (MQ=255)
                                                                                          ttCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGAAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCAcc                                                       >  1:183750/1‑147 (MQ=255)
                                                                                            cAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATTGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCa                                                     <  2:104364/147‑1 (MQ=255)
                                                                                            cAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATTGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCa                                                     <  2:107746/147‑1 (MQ=255)
                                                                                             aaaGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCa                                                     >  8:18928/1‑146 (MQ=255)
                                                                                             aaaGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCa                                                     <  7:18928/146‑1 (MQ=255)
                                                                                                   gacAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGatat                                               <  3:22335/146‑1 (MQ=255)
                                                                                                   gacAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGatat                                               >  4:22335/1‑146 (MQ=255)
                                                                                                   gacAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCaa                                            >  6:182249/1‑149 (MQ=255)
                                                                                                    acAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAAc                                           <  5:155542/149‑1 (MQ=255)
                                                                                                    acAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAAc                                           <  4:192305/149‑1 (MQ=255)
                                                                                                          tcATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTcc                                     >  2:83280/1‑149 (MQ=255)
                                                                                                          tcATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTAGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTcc                                     >  3:200462/1‑149 (MQ=255)
                                                                                                          tcATCTTCTTCTTCAGAGGCAATGCCCGCAAAGATATCGCGTGCGTTCAGTGCGGCGTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCAGGCACACAGGATATCAACCAGTcc                                     >  8:201653/1‑149 (MQ=255)
                                                                                                           cATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCt                                    >  1:91274/1‑149 (MQ=255)
                                                                                                            aTCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTc                                      <  5:182249/146‑1 (MQ=255)
                                                                                                                       cAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGata                        >  5:30347/1‑149 (MQ=255)
                                                                                                                           ggCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATcgc                    <  6:223787/149‑1 (MQ=255)
                                                                                                                            gCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATcgc                    >  2:45458/1‑148 (MQ=255)
                                                                                                                                   ccGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGAt             >  1:242126/1‑148 (MQ=255)
                                                                                                                                    cGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATaa           >  8:167165/1‑149 (MQ=255)
                                                                                                                                     gCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGTACAAGTCGAGTTTCTGCTGAAACTCGCCCTGCACGTTCTCAGCAGCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGAt                 <  3:78287/142‑1 (MQ=255)
                                                                                                                                     gCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCa                                         >  1:152866/1‑118 (MQ=255)
                                                                                                                                     gCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCa                                         <  2:152866/118‑1 (MQ=255)
                                                                                                                                     gCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGAAAACGCGAt                 >  4:78287/1‑142 (MQ=255)
                                                                                                                                      cAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACCTTATCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATc         >  1:74235/1‑149 (MQ=255)
                                                                                                                                         cGATATCTCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCGGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTg      >  6:172769/1‑149 (MQ=255)
                                                                                                                                              tCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGgcgc  <  7:40750/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGCTGAACCTCGCCCTGCACGTTCTCAGCACCGCTGGCACCCAGGATATCAACCAGTCCTGCTTTGTTGATATCGCGATGGATAATCTTGGCGC  >  NC_000913/4455226‑4455515

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