Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 720,070 G→A 5.8% H131Y (CAT→TAT)  speF ← ornithine decarboxylase isozyme, inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913720,0700GA 5.8% 228.9 / 2.7 69H131Y (CAT→TAT) speFornithine decarboxylase isozyme, inducible
Reads supporting (aligned to +/- strand):  ref base G (33/32);  new base A (2/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTCGCAATTCTC  >  NC_000913/719925‑720214
                                                                                                                                                 |                                                                                                                                                
gCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGAcc                                                                                                                                               <  1:499092/149‑1 (MQ=255)
  aaTGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAg                                                                                                                                             >  3:577102/1‑149 (MQ=255)
  aaTGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAg                                                                                                                                             >  1:464721/1‑149 (MQ=255)
   aTGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAgg                                                                                                                                            <  1:94682/149‑1 (MQ=255)
    tGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATAACCAggg                                                                                                                                           >  1:122618/1‑149 (MQ=255)
    tGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATAACCAggg                                                                                                                                           >  3:85933/1‑149 (MQ=255)
                  tCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCt                                                                                                                             <  4:737215/149‑1 (MQ=255)
                      gAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTa                                                                                                                         >  4:514110/1‑149 (MQ=255)
                         tCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCt                                                                                                                      >  3:153600/1‑149 (MQ=255)
                           agcagATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCttg                                                                                                                    >  4:372059/1‑149 (MQ=255)
                            gcagATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCttga                                                                                                                   >  3:623437/1‑149 (MQ=255)
                            gcagATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCttga                                                                                                                   >  4:305580/1‑149 (MQ=255)
                            gcagATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCttga                                                                                                                   >  1:380110/1‑149 (MQ=255)
                               gATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCttgattg                                                                                                                <  2:588849/149‑1 (MQ=255)
                                     ccATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGAc                                                                                                              >  4:347224/1‑145 (MQ=255)
                                     ccATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGAc                                                                                                              <  3:347224/145‑1 (MQ=255)
                                     ccATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATaa                                                                                                          >  3:191681/1‑149 (MQ=255)
                                      cATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATAACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAAt                                                                                                         <  2:122618/149‑1 (MQ=255)
                                      cATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATAACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAAt                                                                                                         <  4:85933/149‑1 (MQ=255)
                                       aTCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATc                                                                                                        <  4:623437/149‑1 (MQ=255)
                                       aTCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATc                                                                                                        >  1:142702/1‑149 (MQ=255)
                                          gCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGAc                                                                                                     <  4:439999/149‑1 (MQ=255)
                                          gCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGAc                                                                                                     <  3:122275/149‑1 (MQ=255)
                                              cGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGt                                                                                                 <  2:380110/149‑1 (MQ=255)
                                                    cGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCgg                                                                                           >  3:309070/1‑149 (MQ=255)
                                                     gTTGCACAAGTCGGCACGGAACAGTGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGga                                                                                          >  3:639988/1‑149 (MQ=255)
                                                     gTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCACTc                                                                                                                                      >  1:237678/1‑105 (MQ=255)
                                                     gTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATc                                                                                                                                      <  2:237678/105‑1 (MQ=255)
                                                     gTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGga                                                                                          >  3:610122/1‑149 (MQ=255)
                                                      ttGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGgaa                                                                                         <  1:295333/149‑1 (MQ=255)
                                                       tGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGgaag                                                                                        <  2:641008/149‑1 (MQ=255)
                                                          acaAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGaaa                                                                                     >  1:129684/1‑149 (MQ=255)
                                                          acaAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGaaa                                                                                     <  4:153600/149‑1 (MQ=255)
                                                          acaAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGaaa                                                                                     <  1:654042/149‑1 (MQ=255)
                                                          acaAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGaaa                                                                                     >  1:677262/1‑149 (MQ=255)
                                                          acaAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGaaa                                                                                     <  2:129684/149‑1 (MQ=255)
                                                          acaAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGCCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGaaa                                                                                     >  1:138963/1‑149 (MQ=255)
                                                              gTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCTACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTg                                                                                 <  4:258914/149‑1 (MQ=255)
                                                              gTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTg                                                                                 >  2:261663/1‑149 (MQ=255)
                                                              gTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTg                                                                                 <  2:216740/149‑1 (MQ=255)
                                                               tCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTgg                                                                                >  2:445872/1‑149 (MQ=255)
                                                                              gcCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCACAGTTGAGTTTc                                                                 <  4:639988/149‑1 (MQ=255)
                                                                                cTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCAt                                                               >  1:299106/1‑149 (MQ=255)
                                                                                cTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCAt                                                               >  1:440296/1‑149 (MQ=255)
                                                                                  cACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATaa                                                             <  2:464721/149‑1 (MQ=255)
                                                                                  cACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCAGAAGAAAGGTGGGCGCAGTTGAGTTTCATaa                                                             <  2:677262/149‑1 (MQ=255)
                                                                                   aCCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGcgcg                                                                                            >  3:591221/1‑117 (MQ=255)
                                                                                   aCCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGcgcg                                                                                            <  4:591221/117‑1 (MQ=255)
                                                                                     cAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAAtgg                                                          >  3:118952/1‑149 (MQ=255)
                                                                                        aGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAAtggctg                                                       <  1:261663/149‑1 (MQ=255)
                                                                                            ttCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAg                                                   >  3:328062/1‑149 (MQ=255)
                                                                                            ttCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAg                                                   <  4:577102/149‑1 (MQ=255)
                                                                                              ccACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCg                                                 <  3:514110/149‑1 (MQ=255)
                                                                                               cacaAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCgg                                                <  4:118952/149‑1 (MQ=255)
                                                                                                    aCTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCt                                           <  3:400436/149‑1 (MQ=255)
                                                                                                        aTTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCt                                           <  3:52110/145‑1 (MQ=255)
                                                                                                        aTTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCt                                           >  4:52110/1‑145 (MQ=255)
                                                                                                             ccgccgGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGa                                  <  2:740017/149‑1 (MQ=255)
                                                                                                                   gATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATaa                            >  1:461747/1‑149 (MQ=255)
                                                                                                                       gCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAAtt                        <  1:235488/149‑1 (MQ=255)
                                                                                                                       gCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAAtt                        >  3:131624/1‑149 (MQ=255)
                                                                                                                          aCGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATtctc                     >  4:470522/1‑149 (MQ=255)
                                                                                                                              aaaaaTTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCtcgtcg                 >  2:313208/1‑149 (MQ=255)
                                                                                                                               aaaaTTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAAtggctggc                                                     >  1:289664/1‑112 (MQ=255)
                                                                                                                               aaaaTTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAAtggctggc                                                     <  2:289664/112‑1 (MQ=255)
                                                                                                                               aaaaTTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTcgc                <  4:540358/149‑1 (MQ=255)
                                                                                                                                 aaTTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGa              >  1:438922/1‑149 (MQ=255)
                                                                                                                                      gccgccCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTCGc         >  1:673259/1‑149 (MQ=255)
                                                                                                                                          cccTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTCGCAAtt     <  3:372059/149‑1 (MQ=255)
                                                                                                                                             tgatgaCCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTCGCAATtctc  <  3:305580/149‑1 (MQ=255)
                                                                                                                                             tgatgaCCAGAGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTg                                                                       >  1:94446/1‑80 (MQ=255)
                                                                                                                                             tgatgaCCAGAGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTg                                                                       <  2:94446/80‑1 (MQ=255)
                                                                                                                                                tgaCCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTCGCaa       <  2:136867/141‑1 (MQ=255)
                                                                                                                                                tgaCCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTCGCaa       >  1:136867/1‑141 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTCGCAATTCTC  >  NC_000913/719925‑720214

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