Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,990,679 T→G 7.2% S30R (AGC→CGC)  hemC ← hydroxymethylbilane synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,990,6790TG 7.2% 197.0 / 2.9 81S30R (AGC→CGC) hemChydroxymethylbilane synthase
Reads supporting (aligned to +/- strand):  ref base T (35/39);  new base G (2/5);  total (37/44)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCGATATCGGCGCGATTTTCGAGGAGCGCGACTTCCAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCA  >  NC_000913/3990532‑3990825
                                                                                                                                                   |                                                                                                                                                  
cGGCGATATCGGCGCGATTTTCGAGGAGCGCGACTTCCAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCCAGAATCACATCGCCGCGCGTCACCAGCGGGACCCGTTCACCGACCAGGCCCGGATGGCTc                                                                                                                                                   >  2:869191/1‑149 (MQ=255)
  gCGATATCGGCGCGATTTTCGAGGAGCGCGACTTCCAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCCTCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGc                                                                                                                                                 >  1:40559/1‑149 (MQ=255)
  gCGATATCGGCGCGATTTTCGAGGAGCGAGACTTCCAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCGGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGc                                                                                                                                                 <  1:585789/149‑1 (MQ=255)
    gATATCGGCGCGATTTTCGAGGAGCGCGACTTCCAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCa                                                                                                                                               <  1:903996/149‑1 (MQ=255)
    gATATCGGCGCGATTTTCGAGGAGAGCGAATTCCAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGAATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACTACCAGGCCCGGATGGCTCGCCa                                                                                                                                               >  1:1015302/1‑149 (MQ=255)
         cGGCGCGATTTTCGAGGAGCGCGACTTCCAGCTCTTTTACAAATAAGCCTTTTCCGCCTTCTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACTACCAGGCCCGGATGGCGAGCCATccgg                                                                                                                                          <  1:242561/149‑4 (MQ=255)
         cGGCGCGAGTTTCGAGGAGCGCGACTTCCAGCGCTTTTACAAATACGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCGAGCCATccgg                                                                                                                                          <  1:844259/149‑4 (MQ=255)
            cgcgATTTTCGAGGAGCGCGACTTCCAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATCGCTCGCCATCAACTTg                                                                                                                                       >  1:668749/1‑149 (MQ=255)
            cgcgATTTTCGAGGAGCGCGACTTCCAGCTCTTTTACAAATAAGACTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTg                                                                                                                                       >  1:858314/1‑149 (MQ=255)
             gcgATTTTCGAGGAGCGCGACTTCCAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGt                                                                                                                                      <  2:9921/149‑1 (MQ=255)
             gcgATTTTCGAGGAGCGCGACTTCCAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGt                                                                                                                                      <  2:119989/149‑1 (MQ=255)
               gATTATCGAGTATCGAGACTTCCAGCTCTTTTACAAATAATACTTTTCCGCCTACTATCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGGACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCt                                                                                                                                    <  2:1038243/149‑1 (MQ=255)
                 ttttCGAGGAGCGCGACTTCCAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCttt                                                                                                                                  <  1:1156658/149‑1 (MQ=255)
                 ttttCGAGGAGCGCGACTTCCAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCttt                                                                                                                                  <  2:380113/149‑1 (MQ=255)
                 ttgtCGCGGCGCCCGACTTCCCGCTCTTATACAAATAAGCCTTTGCCGCCTACTTTCGCCAGCGGCGTAGCAAGAATCACCTCGCCGCGCGTCACCATCGGTACCAGTTCTACGACCAGGCCCGGATGGCGCGCCATCAACTTGTCttt                                                                                                                                  <  2:766295/146‑1 (MQ=255)
                 ttctCGAGGAGCGCGACTTCCAGCTCCTTTACAAATAAGCCGTTTCCGCCTACATTCGCCAGCGGCGTATCAAGCATCACCTCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCttt                                                                                                                                  <  2:191562/146‑1 (MQ=255)
                  tttCGAGGAGCGCGACTTCCAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCCGCGGCGTATCAAGAATCCCATCGCCGCGCGTCACCATCGGTACCAGGTCAACGACCAGGCCCGGATGGCTCGCCATCACCTTGTcgtgg                                                                                                                                 >  2:794829/1‑145 (MQ=255)
                  tttCGAGGAGCGCGACTTCCAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTg                                                                                                                                 >  1:486863/1‑149 (MQ=255)
                   ttCGAGGAGCGCGACTTCCAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGa                                                                                                                                <  2:213279/149‑1 (MQ=255)
                   ttCGAGGAGCGCGACTTCCAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGa                                                                                                                                >  1:1155516/1‑149 (MQ=255)
                   ttCGAGGAGCGCGACTTCCAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGa                                                                                                                                <  2:176723/149‑1 (MQ=255)
                      gaggagCGCGACTTCCAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTGGTCTTTGACAt                                                                                                                             >  1:347211/1‑149 (MQ=255)
                       aggagCGCGACTTCCAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATa                                                                                                                            >  1:279094/1‑149 (MQ=255)
                           gcgcGACTTCCAGATCTTTAACAAATACGCCTTTTCCACCTACTTTCGCCAGCGGCGTATCAAGAATCACAGCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCGCGCCATCAACTTGTCTTTGACATAgtgt                                                                                                                        <  2:343411/149‑1 (MQ=255)
                              cgACTTCCAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCCGCGGCGTATCAAGAATCACATCGCCGCGCGTCACGATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCACCTTGTCTTTGACATTGTGtgcc                                                                                                                     >  2:1038179/1‑149 (MQ=255)
                                    ccAGCTCTTTTACAAATAAGCATTTTCCGCCTACTTTCGCCCGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTAGCCATCAACTTGTCTTTGACATAGTGTGCCTTCCag                                                                                                               >  1:534599/1‑149 (MQ=255)
                                      atctctTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATATGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGACagag                                                                                                             <  2:289872/147‑1 (MQ=255)
                                       gCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGt                                                                                                            <  2:908158/149‑1 (MQ=255)
                                       gCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGt                                                                                                            <  1:1087451/149‑1 (MQ=255)
                                        ctctTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTg                                                                                                           <  1:892378/149‑1 (MQ=255)
                                          ctTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGTGATAGTGTGCCTGCGAGAGTGCa                                                                                                         >  2:746646/1‑149 (MQ=255)
                                           ttttACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCTCCGCGCGTCACCATCGGTTCCAGTTAAACGCCCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGcta                                                                                                        >  1:324975/1‑147 (MQ=255)
                                            tttACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCTGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAg                                                                                                       <  1:228766/149‑1 (MQ=255)
                                             ttACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTATGACATAGTGTGCCTGCCAGAGTGCAAGt                                                                                                      >  1:714136/1‑149 (MQ=255)
                                                 aaaTAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGCCCAGGCCCGGATGGCTCGCGAGCAACTTGTCTTTGACATAGTGTGCCTGCCAGCGTGCAAGTgggg                                                                                                  >  2:832430/1‑148 (MQ=255)
                                                 aaaTAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCCGGCCCGGATGGCTCGCCATCATCTTGTCTTTGACATAGTGTGACTGCCAGAGTGCATGTGGGc                                                                                                  >  2:998462/1‑149 (MQ=255)
                                                    tAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTCCCAGTTCAACGACCAGGCCCGGATGGCTCGCCCTCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCttt                                                                                               >  1:921507/1‑149 (MQ=255)
                                                     aaGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTg                                                                                              <  1:1163906/149‑1 (MQ=255)
                                                        ccgtttCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACACTCGGTACGAGGTCAACGCCCAGGCCCGGATGGGGCGACAGGACAGTGTCTTTGACATAGTGTGCCGGCCAGAGTGCAAGTGGGCTTGGGCg                                                                                           >  2:803024/4‑149 (MQ=255)
                                                              ccGCCTACTTTCGCCAGCGGCGTATCAAGAATCACCTCGCCGCGCGTCACCATCCGTACCAGTTCAACGACCAGGCCCGGATGTCGCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGGGTGGc                                                                                     <  2:858602/149‑1 (MQ=255)
                                                              ccGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCTCCGCGCTTCACCATCTGTACAAGTTAAACGACCAGGACCGGATGGATCTCCATCAACTTGAATTTTACATAGTGTGCCTGCCAGAGTGCAAGTGGGCATTGGCGTGTGGc                                                                                     >  2:100256/1‑149 (MQ=255)
                                                              ccGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGc                                                                                     <  1:548904/149‑1 (MQ=255)
                                                                   tACTTTCTCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTc                                                                                <  2:581577/149‑1 (MQ=255)
                                                                       ttCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCTGATCAACGACCAGGCCCGGATGGCTCGCCAGCAACTTTTCTTTGACATAGTGTGCCAGCCAGAGGGCAAGTGGGCTTTGGCGTGTGGCAATTCTTaa                                                                            >  1:293593/1‑149 (MQ=255)
                                                                       ttCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTaa                                                                            >  1:1681/1‑149 (MQ=255)
                                                                       ttCGCCAGCGGCGTAACAAGAATCAAATAGCCGCGCGTCACCATCGGTACCAGTTCAACGAGCAGGCCCGGATGGCTCGCCATCAACATGTCGTTGACATAGTGTGACTGCCACCGTGCAAGTGGGCTCTGCAGTGTGGCAAtagtaaa                                                                            >  2:1045349/1‑143 (MQ=255)
                                                                          gCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCAACGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTACCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAAc                                                                         <  1:457484/149‑1 (MQ=255)
                                                                            cAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACAt                                                                       >  2:386464/1‑149 (MQ=255)
                                                                              ccggcgTATCAAGAATCACATCGCCGCGCGTCACCATCGGTCCCAGTTCCACGACCACGCCCGGAGGGCTCGCCATCAACTTGGCTTTGACATAGTGCGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTg                                                                     <  2:513879/148‑1 (MQ=255)
                                                                                ggcgTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTc                                                                   <  2:294814/149‑1 (MQ=255)
                                                                                ggcgTATCAAGAATCACATCGCCGCGCGGCACGATCTGAACCAGATCACCGACCAGGCCCGGATGGCTCGCCATCAACTAGTCTATGACATAGAGTGCCTGCCAGAGTGGAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTc                                                                   <  2:638026/149‑1 (MQ=255)
                                                                                 gagTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCt                                                                  <  1:771914/147‑1 (MQ=255)
                                                                                    tATCAAGAATCACATCGCCGCGCGTCACCAACGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGAGGGCTTTGGCGGGTGGCAATTCTAAAAACATTGTCTacc                                                               >  1:529846/1‑147 (MQ=255)
                                                                                      tCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACAt                                                             >  2:946754/1‑149 (MQ=255)
                                                                                      tCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGGTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACAt                                                             >  1:8822/1‑149 (MQ=255)
                                                                                       cAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGCCCAGGCCCGGATGGCTCGCCATCAACTTGGCGTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTATGGCGTGTGGCAATTCTTAAAACATTGTCTAACATg                                                            >  2:797370/1‑149 (MQ=255)
                                                                                       cAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATg                                                            >  1:381357/1‑149 (MQ=255)
                                                                                         aGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACGAGGCGCGGATGGCTCTCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCt                                                          <  1:775577/149‑1 (MQ=255)
                                                                                           aaTCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGGGCATGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTg                                                        >  1:1027816/1‑149 (MQ=255)
                                                                                             tCACATCGCCGCGCGTCACCATTGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTCGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGtt                                                      >  1:816804/1‑149 (MQ=255)
                                                                                             tCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGtt                                                      <  2:161278/149‑1 (MQ=255)
                                                                                              cacaTCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTa                                                     <  2:1155516/149‑1 (MQ=255)
                                                                                                caTCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTAcc                                                   >  2:337886/1‑149 (MQ=255)
                                                                                                    gccgcGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCa                                               >  1:207921/1‑149 (MQ=255)
                                                                                                      cgcgcgTCACCATCGGTACCAGTTCAACGACCAGGCCCGGAGGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCAtt                                             <  2:1084756/149‑1 (MQ=255)
                                                                                                       gcgcgTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTa                                            <  2:668749/149‑1 (MQ=255)
                                                                                                          cgTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTatca                                         <  2:486863/149‑1 (MQ=255)
                                                                                                             cACCATCGGTACCAGTTCACCGACCAGGCGCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCAGCCAGAGTGCAAGTGGGCTTTGGCGTTTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATcc                                      <  1:643021/149‑1 (MQ=255)
                                                                                                              aCCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCg                                     >  1:271296/1‑149 (MQ=255)
                                                                                                                  tgggTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGGGTGGCAATTCTTAAAACATTGTCTAACCTGCTTGTTACCGTCATTATCATCCGTGGt                                 <  1:794829/147‑1 (MQ=255)
                                                                                                                      tACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCAt                             <  2:1010141/149‑1 (MQ=255)
                                                                                                                       aCCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATc                            <  2:185945/149‑1 (MQ=255)
                                                                                                                           gTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTaa                        >  2:1056859/1‑149 (MQ=255)
                                                                                                                               aaCGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATc                    <  2:279094/149‑1 (MQ=255)
                                                                                                                                  gACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCtt                 <  1:1169089/149‑1 (MQ=255)
                                                                                                                                    ccAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGc               <  1:1056859/149‑1 (MQ=255)
                                                                                                                                        gCCCGGATGGCTCGCCATCAACTTGTCTTTGACAAAGTGTGCCGGCCAGAGTGCAAGTGGGCTTTGGCGGGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCaga           <  2:654826/149‑1 (MQ=255)
                                                                                                                                          gttgATGGCGAGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGAAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGATTGTTACCGTCATTATCATACATGTGACATCCTAACATACTTGCCAGAGt         >  2:821026/4‑149 (MQ=255)
                                                                                                                                            ggATGGCTCGCCATCAACTTGTCTTTGACTTAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTTTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCAATCCTAACATCCTTGCCAGAGTGa       <  2:1015302/149‑1 (MQ=255)
                                                                                                                                                ggCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTAGTTACCGTCATTATCATCCGTGGACCATCCTAACTACCTTGCCAGAGAGATGTc   >  2:30280/1‑149 (MQ=255)
                                                                                                                                                 gCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTAGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCa  <  1:296718/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGGCGATATCGGCGCGATTTTCGAGGAGCGCGACTTCCAGCTCTTTTACAAATAAGCCTTTTCCGCCTACTTTCGCCAGCGGCGTATCAAGAATCACATCGCCGCGCGTCACCATCGGTACCAGTTCAACGACCAGGCCCGGATGGCTCGCCATCAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCA  >  NC_000913/3990532‑3990825

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: ‑