Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 926,168 C→T 9.4% intergenic (‑197/+57) serW ← / ← infA tRNA‑Ser/translation initiation factor IF‑1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913926,1680CT 9.4% 321.8 / 20.6 106intergenic (‑197/+57)serW/infAtRNA‑Ser/translation initiation factor IF‑1
Reads supporting (aligned to +/- strand):  ref base C (52/44);  new base T (6/4);  total (58/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 = 1.00e-01

Only 100 of 106 total aligned reads displayed.

AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTAATGATAATTATTGTTTGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGAACGGTACCTTGCATTTCAATATTGTCTTCTTTGGCC  >  NC_000913/925905‑926441
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                 
aaCCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTAATGATAATTATTGTTTGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAAcc                                                                                                                                                                                                                                                                         <  2:531860/274‑1 (MQ=255)
                          acTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTAATGATAATTATTGTTTGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGc                                                                                                                                                                                                                                              >  2:464661/1‑275 (MQ=255)
                                      tGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTAATGATAATTATTGTTTGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCgg                                                                                                                                                                                                                                   >  4:69935/1‑274 (MQ=255)
                                      tGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTAATGATAATTATTGTTTGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCgg                                                                                                                                                                                                                                   >  2:449127/1‑274 (MQ=255)
                                                                             gCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTAATGATAATTATTGTTTGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGc                                                                                                                                                                                           <  4:449126/275‑1 (MQ=255)
                                                                                                                                                aTTTTAATGATAATTATTGTTTGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTgcg                                                                                                                        <  2:70008/275‑1 (MQ=255)
                                                                                                                                                                   tttGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGAGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTAGGGGGTCAGTTCAACAGTCACTTGTACCCCCGTCAGGTGGGGGAGGNAGGTTNNGGGCAgttttaaggtgaggtgg                                                                                                     >  3:9654/1‑262 (MQ=255)
                                                                                                                                                                     tGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGAtgtgtg                                                                                                    >  2:147638/1‑274 (MQ=255)
                                                                                                                                                                     tGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGAGACTACGGAAAACAATGCGGCCTTTGCTCAGGTCGGACGGGGTCAGTTCAAAAGTCACTTTTCCGCCCGTCCGGATGCGGCTTTAATTTTTGCGCATTTTACCGgagcttgttg                                                                                                    >  4:80760/1‑267 (MQ=255)
                                                                                                                                                                         cATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTaa                                                                                                                                                                                                                                                                            <  5:888964/102‑1 (MQ=255)
                                                                                                                                                                         cATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGtt                                                                                                                                                                                                                                                          <  7:520949/120‑1 (MQ=255)
                                                                                                                                                                         cATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGtt                                                                                                                                                                                                                                                          >  8:520949/1‑120 (MQ=255)
                                                                                                                                                                              tGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttcttt                                                                                                                                                                                                                                                      <  8:216545/119‑1 (MQ=255)
                                                                                                                                                                              tGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttcttt                                                                                                                                                                                                                                                      >  7:216545/1‑119 (MQ=255)
                                                                                                                                                                               gTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTaaa                                                                                                                                                                                                                               <  8:575161/141‑1 (MQ=255)
                                                                                                                                                                               gTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTaaa                                                                                                                                                                                                                               >  7:575161/1‑141 (MQ=255)
                                                                                                                                                                                    cAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTaa                                                                                                                                                                                                                                                                            <  7:815076/91‑1 (MQ=255)
                                                                                                                                                                                    cAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTaa                                                                                                                                                                                                                                                                            >  8:815076/1‑91 (MQ=255)
                                                                                                                                                                                     aaGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCg                                                                                   <  4:147637/275‑1 (MQ=255)
                                                                                                                                                                                          gtgCAGTATATCTACATCGAGACAAGTTACGGCCTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGCCAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTGTTCGCACGTCAGGCTTCGGATGTAGTTTTTCAGCATTTTACCGTACATGTTTGAAGTACCCCCGCTACCGTTTTc                                                                              >  4:358526/1‑275 (MQ=255)
                                                                                                                                                                                               gTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAAct                                                                         <  2:358527/275‑1 (MQ=255)
                                                                                                                                                                                                     cTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTAcgc                                                                   >  2:78854/1‑275 (MQ=255)
                                                                                                                                                                                                     cTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTAAGGTCGTACGGGGTCAGGTATACAGTAACTTTGTCGCCCGTCAGGCTGCGGCTGTAGTTTTTGCGGATTTTACCGGAGCTGTGTGCAGTAACCACGGCACACGTGTCTAACTCTatta                                                                   >  4:536778/1‑272 (MQ=255)
                                                                                                                                                                                                       aCATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttc                                                                                                                                                                                                                                                         <  7:928450/91‑1 (MQ=255)
                                                                                                                                                                                                        cATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGa                                                                >  2:300143/1‑275 (MQ=255)
                                                                                                                                                                                                        cATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGa                                                                >  2:476118/1‑275 (MQ=255)
                                                                                                                                                                                                          tCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGccc                                                                                                                                                                                                                                            >  2:276017/1‑101 (MQ=255)
                                                                                                                                                                                                          tCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCa                                                                                                                                                                                                                                           <  8:915469/102‑1 (MQ=255)
                                                                                                                                                                                                          tCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTCCGGGGTCAGTTAAACAGTCACTTTGTCGCCCGTCCGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGTTTTGTGAAGTAACCAGGTGCCCGTTTTCTAAATCTACGCGGaa                                                               >  4:3156/1‑274 (MQ=255)
                                                                                                                                                                                                          tCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGaa                                                               >  2:356571/1‑274 (MQ=255)
                                                                                                                                                                                                          tCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTAAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGaa                                                               >  2:172205/1‑274 (MQ=255)
                                                                                                                                                                                                            gagaCAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACa                                                             <  2:69936/274‑1 (MQ=255)
                                                                                                                                                                                                              gaaaaGTTACGGACTTATACTGCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTTGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGt                                                                                                                                                                           <  2:98198/159‑1 (MQ=255)
                                                                                                                                                                                                              gaCAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCa                                                                                                                                                                                                                                        <  4:276016/101‑1 (MQ=255)
                                                                                                                                                                                                              gaCAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCa                                                                                                                                                                                                                        >  1:80419/1‑117 (MQ=255)
                                                                                                                                                                                                               ccAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATAACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTATGCGGAACa                                                             <  2:536779/270‑1 (MQ=255)
                                                                                                                                                                                                                cAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATc                                                                                                                                                                                                                         <  6:992982/114‑1 (MQ=255)
                                                                                                                                                                                                                cAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATc                                                                                                                                                                                                                         >  5:992982/1‑114 (MQ=255)
                                                                                                                                                                                                                  aGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGa                                                                                                                                                                                                                    <  3:80419/117‑1 (MQ=255)
                                                                                                                                                                                                                             tataCTTCCAAAGTACTTCATACACATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCTCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTTCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGttt                                           >  4:77577/1‑275 (MQ=255)
                                                                                                                                                                                                                              ataCTTCCAAAGTACTTCATACATATAACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTc                                          <  1:9654/275‑1 (MQ=255)
                                                                                                                                                                                                                                   tCCAAAGTACTTCATACATAACACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGaa                                     <  2:3156/275‑1 (MQ=255)
                                                                                                                                                                                                                                    ccAAAGTACTTCATACATATCACAAGATAAAAAGGCCGGTTaa                                                                                                                                                                                                                                                                            <  6:676590/43‑1 (MQ=255)
                                                                                                                                                                                                                                    ccAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTaa                                                                                                                                                                                                                                                                            >  5:676590/1‑43 (MQ=255)
                                                                                                                                                                                                                                    ccAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCgg                                                                 >  1:7787/1‑246 (MQ=255)
                                                                                                                                                                                                                                    ccAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCGCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGCTGTAGTTTTTGCTCATTTTACCGGAGCTGGGTGCAGTCACCACGTGACCGTTTTCTAACTCTCAGCGGAACATGGTATTAGGCAACGTTTCAAGAAc                                    >  2:381377/1‑275 (MQ=255)
                                                                                                                                                                                                                                        aGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCt                                                                                                                                                                                          <  2:121103/121‑1 (MQ=255)
                                                                                                                                                                                                                                        aGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGt                                                                                                                                                               <  7:343852/148‑1 (MQ=255)
                                                                                                                                                                                                                                        aGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGt                                                                                                                                                               >  8:343852/1‑148 (MQ=255)
                                                                                                                                                                                                                                            cTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGAACGGTACCtt                            >  2:97797/1‑275 (MQ=255)
                                                                                                                                                                                                                                             ttaatacataTAAAAAAATAAAAAGGCCGGTTACACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCGTGCGGTANAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGAACGGTACCTTg                           <  3:127437/272‑1 (MQ=255)
                                                                                                                                                                                                                                                 tacataTCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGa                                                                                                                                                                                                             <  7:902519/93‑1 (MQ=255)
                                                                                                                                                                                                                                                 tacataTCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGa                                                                                                                                                                                                             >  8:902519/1‑93 (MQ=255)
                                                                                                                                                                                                                                                  acataTCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACg                                              >  4:424150/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    atatCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGAACGGTACCTTGCATTTCa                    >  4:391950/1‑275 (MQ=255)
                                                                                                                                                                                                                                                      atCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCa                                                                                                                                                                                                                        >  1:122744/1‑77 (MQ=255)
                                                                                                                                                                                                                                                      atCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGAACGTTAc                               >  1:127437/1‑262 (MQ=255)
                                                                                                                                                                                                                                                      atCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGAACGGTACCt                             >  4:205775/1‑264 (MQ=255)
                                                                                                                                                                                                                                                      atCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGAACGGTACCTTGCATTTCaa                   >  2:364549/1‑274 (MQ=255)
                                                                                                                                                                                                                                                        cacaAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAAc                                                                                                                                                                  >  7:327175/1‑129 (MQ=255)
                                                                                                                                                                                                                                                        cacaAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAAc                                                                                                                                                                  <  8:327175/129‑1 (MQ=255)
                                                                                                                                                                                                                                                          caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCgg                                                                                                                                                                                                                                   >  7:159990/1‑62 (MQ=255)
                                                                                                                                                                                                                                                          caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCgg                                                                                                                                                                                                                                   <  8:159990/62‑1 (MQ=255)
                                                                                                                                                                                                                                                          caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCgg                                                                                                                                                                                                                                   >  7:106917/1‑62 (MQ=255)
                                                                                                                                                                                                                                                          caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCgg                                                                                                                                                                                                                                   <  8:106917/62‑1 (MQ=255)
                                                                                                                                                                                                                                                          caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCATCGa                                                                                                                                                                                                                    <  3:122744/77‑1 (MQ=255)
                                                                                                                                                                                                                                                          caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCttt                                                                                                                                                                                             <  6:270672/100‑1 (MQ=255)
                                                                                                                                                                                                                                                          caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCttt                                                                                                                                                                                             >  5:270672/1‑100 (MQ=255)
                                                                                                                                                                                                                                                          caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGAACGGTACCTTGCa                         <  2:205776/264‑1 (MQ=255)
                                                                                                                                                                                                                                                           aaaaTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAgg                                                                                                                                            >  5:775260/1‑148 (MQ=255)
                                                                                                                                                                                                                                                            aaaTAAAAAGGCCGGTTAAACCGACCTTTTACTCGt                                                                                                                                                                                                                                                           >  8:740667/1‑36 (MQ=255)
                                                                                                                                                                                                                                                            aaaTAAAAAGGCCGGTTAAACCGACCTTTTACTCGt                                                                                                                                                                                                                                                           <  7:740667/36‑1 (MQ=255)
                                                                                                                                                                                                                                                            aaaTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACggg                                                                                                                                                                             <  8:353003/114‑1 (MQ=255)
                                                                                                                                                                                                                                                            aaaTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACggg                                                                                                                                                                             >  7:353003/1‑114 (MQ=255)
                                                                                                                                                                                                                                                             aaTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCgg                                                                 >  2:507921/1‑221 (MQ=255)
                                                                                                                                                                                                                                                              aTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATg                                                                                                                                         >  7:528917/1‑148 (MQ=255)
                                                                                                                                                                                                                                                              aTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATg                                                                                                                                         <  6:775260/148‑1 (MQ=255)
                                                                                                                                                                                                                                                               tAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATg                                                                                                                                                                                                     >  7:1095631/1‑87 (MQ=255)
                                                                                                                                                                                                                                                               tAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATg                                                                                                                                                                                                     <  8:1095631/87‑1 (MQ=255)
                                                                                                                                                                                                                                                               tAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGt                                                                                                                                                       >  8:531771/1‑133 (MQ=255)
                                                                                                                                                                                                                                                               tAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGt                                                                                                                                                       <  7:531771/133‑1 (MQ=255)
                                                                                                                                                                                                                                                                taaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCATTTCAACAGTCACTTTg                                                                                                                                                        >  5:82041/2‑131 (MQ=255)
                                                                                                                                                                                                                                                                taaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTg                                                                                                                                                        <  6:82041/130‑1 (MQ=255)
                                                                                                                                                                                                                                                                aaaaaGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACgg                                                                                                                                                                              <  7:1090976/109‑1 (MQ=255)
                                                                                                                                                                                                                                                                aaaaaGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACgg                                                                                                                                                                              >  8:1090976/1‑109 (MQ=255)
                                                                                                                                                                                                                                                                 aaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGc                                                                                                                                                                                                                                     <  6:1105657/53‑1 (MQ=21)
                                                                                                                                                                                                                                                                 aaaaGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACa                                                             <  4:507920/221‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aaaGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTAc                                                                                                                >  4:526637/1‑169 (MQ=255)
                                                                                                                                                                                                                                                                   aaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGc                                                                                                                                                                                                 <  5:1046275/87‑1 (MQ=39)
                                                                                                                                                                                                                                                                   aaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGc                                                                                                                                                                                                 <  5:744247/87‑1 (MQ=39)
                                                                                                                                                                                                                                                                   aaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGc                                                                                                                                                                                                 >  6:744247/1‑87 (MQ=39)
                                                                                                                                                                                                                                                                   aaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGc                                                                                                                                                                                                 >  6:1046275/1‑87 (MQ=39)
                                                                                                                                                                                                                                                                    aGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAg                                                                                                                                                                                                           >  7:291480/1‑76 (MQ=38)
                                                                                                                                                                                                                                                                    aGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAAc                                                                                                                                                                  >  5:190810/1‑117 (MQ=255)
                                                                                                                                                                                                                                                                    aGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGAACGGTACCTTGCATTTCAATATTGTCTTCTTTgg    <  2:434308/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACaa                                                                                                                                                                                                                           <  5:58648/58‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACaa                                                                                                                                                                                                                           >  6:58648/1‑58 (MQ=255)
                                                                                                                                                                                                                                                                      gCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGga                                                                                                            <  2:526638/169‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGt                                                                                                <  1:96380/181‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGAACGGTACCTTGCATTTCAATATTGTCTTCTTTGGcc  <  4:464660/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                 
AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTAATGATAATTATTGTTTGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGAACGGTACCTTGCATTTCAATATTGTCTTCTTTGGCC  >  NC_000913/925905‑926441

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