Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,174,504 A→T 5.4% L221Q (CTG→CAG)  ycfT ← inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,174,5040AT 5.4% 725.4 / 22.4 205L221Q (CTG→CAG) ycfTinner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (104/89);  new base T (5/6);  total (109/96)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-01

Only 100 of 211 total aligned reads displayed.

GTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAGAATGCGATGGGTGGTGTAGATAGCAATGGTGTTGGAGCCAATCACATTCAGCAGGCTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGC  >  NC_000913/1174232‑1174748
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
                  gcgACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAGAATGCGATGGGTGGTGTAGATAGCAATGGTGTTGGAGCCAATCACAGTCAGCAGGCTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCAc                                                                                                                                                                                                                                   <  4:534692/274‑1 (MQ=255)
                                                gCAGAGTTAAGCTGAATATTTCAACCAGAATGCGATGGGTGGTGTAGATAGCAATGGTGTTGGAGCCAATCACATTCAGCAGGCTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGccc                                                                                                                                                                                                    >  1:14901/1‑274 (MQ=255)
                                                              aaTATTTCAACCAGAATGCGATGGGTGGTGTAGATAGCAATGGTGTTGGAGCCAATCACATTCAGCAGGCTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACg                                                                                                                                                                                                                                       >  2:347982/1‑226 (MQ=255)
                                                                      aaCCATAATGCGATGGGTTGTGCAGATAGCAATGGTGTTGGAGCCAATCACATTCAGCAGGCTGGTGTAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCTGGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGtaat                                                                                                                                                                              <  2:510683/275‑1 (MQ=255)
                                                                                                   cAATGGTGTTGGAGCCAATCACATTCAGCAGGCTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCAtt                                                                                                                                                 >  2:369045/1‑275 (MQ=255)
                                                                                                    aaTGGTGTTGGAGCCAATCACATTCAGCAGGCTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCAtt                                                                                                                                                 <  2:440454/274‑1 (MQ=255)
                                                                                                                                  gCTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCTGCTGATggc                                                                                                                                                                                                                                                <  7:612349/149‑1 (MQ=255)
                                                                                                                                   cTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGAAAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGATGTAATAAGGCAAATTGCTGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAACTTAACCg                                                                                                                  >  2:595853/1‑274 (MQ=255)
                                                                                                                                            cgcATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTCACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTTTAATAAGGCCAATTTCGGATCACACTGTTCATTCCCCCCCACGGCGTGGGCACAAAATTAACCGCCACACCca                                                                                                         >  4:488036/1‑274 (MQ=255)
                                                                                                                                               aTACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCa                                                                                                                                                                                                                                    >  5:782378/1‑148 (MQ=255)
                                                                                                                                                tACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCTGCTGATggcttc                                                                                                                                                                                                                                             >  7:429009/1‑135 (MQ=255)
                                                                                                                                                tACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCTGCTGATggcttc                                                                                                                                                                                                                                             <  8:429009/138‑4 (MQ=255)
                                                                                                                                                tACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCg                                                                                                                  <  2:488038/261‑1 (MQ=255)
                                                                                                                                                 aCCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCa                                                                                                                                                                                                                                    <  6:693963/146‑1 (MQ=255)
                                                                                                                                                  ccGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACctc                                                                                                                                                                                                                                <  8:608549/149‑1 (MQ=255)
                                                                                                                                                    gAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACctct                                                                                                                                                                                                                               >  5:387234/1‑148 (MQ=255)
                                                                                                                                                    gAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCtt                                                                                                                                                                                                                              >  8:717755/1‑149 (MQ=255)
                                                                                                                                                     aaTCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCa                                                                                                                                                                                               <  2:581370/179‑1 (MQ=255)
                                                                                                                                                      accGTTGCGCGTATGGAAAAAACAGCTTCATGATAACCCAAACGGATACCAGCGACACAAGAAGCGAGATAGTAAACAACCAGGCACCGACCGCAAGAACGGTAAGAAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTacaaac                                                                                              <  4:369043/273‑1 (MQ=255)
                                                                                                                                                      aTCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTa                                                                                                                                                                                                                            >  7:306110/1‑149 (MQ=255)
                                                                                                                                                       tCGTTGCTAGTATTAATAAAACAGCTTCAGGATCACCACAAGCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTacaaac                                                                                              <  2:556247/274‑1 (MQ=255)
                                                                                                                                                        cGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTaa                                                                                                                                                                                                                           <  6:2555/148‑1 (MQ=255)
                                                                                                                                                           tGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAacac                                                                                                                                                                                                                        >  6:336760/1‑148 (MQ=255)
                                                                                                                                                              tCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGt                                                                                                                                                                                                                    <  8:76331/149‑1 (MQ=255)
                                                                                                                                                              tCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGt                                                                                                                                                                                                                    >  7:13131/1‑149 (MQ=255)
                                                                                                                                                              tCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGAGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCGAGAACGGTAAGTAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGt                                                                                                                                                                                                                    >  8:555800/1‑149 (MQ=255)
                                                                                                                                                               cGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGAGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCACCGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATTCGCCAAGGCTGTAATAAGGCAAATTCCGGATCACACTGTTGATTCCGCAACACGGCGTGGGAACGAAATGAACAGCCCCACTCGGCAGTACAAAACAGGCGAAt                                                                                     >  1:228642/1‑275 (MQ=255)
                                                                                                                                                                 tATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCa                                                                                                                                                                                                                  >  5:436667/1‑148 (MQ=255)
                                                                                                                                                                       tAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTAttgt                                                                                                                                                                                                           <  7:494300/149‑1 (MQ=255)
                                                                                                                                                                       tAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTAttg                                                                                                                                                                                                            >  6:790224/1‑148 (MQ=255)
                                                                                                                                                                        aaaaCAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCttt                                                                                                                                                                                                                             >  2:654278/1‑130 (MQ=255)
                                                                                                                                                                        aaaaCAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTAttgt                                                                                                                                                                                                           <  6:169127/148‑1 (MQ=255)
                                                                                                                                                                          aaCAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCTGCTGATggc                                                                                                                                                                                                                                                >  5:315308/1‑109 (MQ=255)
                                                                                                                                                                            cAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAaca                                                                                                                                                                                                                         <  4:654276/130‑1 (MQ=255)
                                                                                                                                                                             aGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTgcgc                                                                                                                                                                                                      >  5:225121/1‑148 (MQ=255)
                                                                                                                                                                               cTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGGGGGAACGAAATTAACCGCCACACCCAGAAGTACAAACAAGGACAATAAGGGCaagccccacc                                                                     >  4:422908/1‑266 (MQ=255)
                                                                                                                                                                                    aaaagcaagaaaTCGATAACCGCGACAACAGAAGCGAAATAGCAAACAACCAGGCACCGACCGCCAGAAAGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCAAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCt                                                                 <  4:205888/264‑1 (MQ=255)
                                                                                                                                                                                     gaaccacacgaacGAAACCAGCGACAAANGNAGGAAAATAGTAAACAACCAGACCCCGACCGACAGAACGGTAAGCNAAGAAGCCCTCAACAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCNGGCGCCAAGGCTGTAATGAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaa                                                               <  3:81307/263‑1 (MQ=255)
                                                                                                                                                                                         accacAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGcc                                                                                                                                                                                         >  7:231908/1‑149 (MQ=255)
                                                                                                                                                                                          ccacAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCTGCTGATggcttc                                                                                                                                                                                                                                             <  5:860459/96‑4 (MQ=255)
                                                                                                                                                                                          ccacAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGcc                                                                                                                                                                                         >  6:547149/1‑148 (MQ=255)
                                                                                                                                                                                           cacaATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCaa                                                                                                                                                                                       >  7:260275/1‑149 (MQ=255)
                                                                                                                                                                                             caATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGcc                                                                                                                                                                                         <  4:698566/145‑1 (MQ=255)
                                                                                                                                                                                                 cGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTg                                                                                                                                                                                  >  5:685249/1‑148 (MQ=255)
                                                                                                                                                                                                     aCCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACg                                                                                                                                                                                                                                       >  3:52483/1‑91 (MQ=255)
                                                                                                                                                                                                       cAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAg                                                                                                                                                                           <  7:713761/149‑1 (MQ=255)
                                                                                                                                                                                                       cAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAg                                                                                                                                                                           >  8:393692/1‑149 (MQ=255)
                                                                                                                                                                                                       cAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAg                                                                                                                                                                           >  8:686029/1‑149 (MQ=255)
                                                                                                                                                                                                        aGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAg                                                                                                                                                                           >  5:873807/1‑148 (MQ=255)
                                                                                                                                                                                                        aGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGAc                                                                     >  2:205889/1‑250 (MQ=255)
                                                                                                                                                                                                         gCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTCAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCTCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGCCGGCTAAAAATTTTACATACCACGaaa                                           >  1:183249/1‑275 (MQ=255)
                                                                                                                                                                                                         gCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCAc                                                                                                                                                                                                                                   <  1:52484/91‑1 (MQ=255)
                                                                                                                                                                                                            acaacaGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACg                                                                                                                                                                                                                                       >  2:36066/1‑84 (MQ=255)
                                                                                                                                                                                                                cagcagCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCAc                                                                                                                                                                                                                                   <  4:36065/84‑1 (MQ=255)
                                                                                                                                                                                                                  gcagcGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATcaca                                                                                                                                                          >  2:718518/1‑155 (MQ=255)
                                                                                                                                                                                                                   cagcGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCgg                                                                                                                                                                <  6:873807/148‑1 (MQ=255)
                                                                                                                                                                                                                        agaTAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATcaca                                                                                                                                                          <  8:280872/149‑1 (MQ=255)
                                                                                                                                                                                                                         gaTAGTAAACACCCAGGCACCGACCGCCAGAACTGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCa                            <  1:21265/274‑1 (MQ=255)
                                                                                                                                                                                                                          aTAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGTGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACt                                                                                                                                                        <  7:555800/149‑1 (MQ=255)
                                                                                                                                                                                                                           tAGTAAACAACCAGGCACCGACAGCCCGAACGGCACGCAAAGAAGCCACCAGCAGATGACGCAGCAACGGCCCCACACTAGCACTGATCATTATTCTTGCGCCACACCATGCTCCAACGCTGTAATAAGGCAAATTGCGGATCACACCGTCATTCCCCACGCCGGCGTGGTGACGATATTGACCGACACACTCAGCAGCACAAACAACGCGAACTGTCACAGCCACCGCCGGCCAAAAATTTAACTTCCCCAGAAATAGACAATTAACGcccac                         >  4:84254/1‑270 (MQ=255)
                                                                                                                                                                                                                            aGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTg                                                                                                                                                       >  5:287005/1‑148 (MQ=255)
                                                                                                                                                                                                                              tAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAgg                                                                                                                                                                                     >  7:122716/1‑116 (MQ=255)
                                                                                                                                                                                                                              tAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCaa                                                                                                                                                                       >  6:331298/1‑130 (MQ=255)
                                                                                                                                                                                                                              tAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCaa                                                                                                                                                                       <  5:331298/130‑1 (MQ=255)
                                                                                                                                                                                                                                aacaacCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCAt                                                                                                                                                  <  7:763754/149‑1 (MQ=255)
                                                                                                                                                                                                                                 acaacCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCAtt                                                                                                                                                 <  8:120568/149‑1 (MQ=255)
                                                                                                                                                                                                                                 acaacCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCAt                                                                                                                                                  >  5:639021/1‑148 (MQ=255)
                                                                                                                                                                                                                                  caacCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTc                                                                                                                                                <  8:646560/149‑1 (MQ=255)
                                                                                                                                                                                                                                   aacCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTc                                                                                                                                                >  5:819261/1‑148 (MQ=255)
                                                                                                                                                                                                                                     ccAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACCGGTCATTATTGTTGCGCCACACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGTAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGAATACAGATACcaca               >  1:138839/1‑275 (MQ=255)
                                                                                                                                                                                                                                     ccAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATACGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGTAATTCACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCCGACGGCTAAAAATTATACATACAACGAAATAGACAATTAAATAATAAAGATAcacca               >  4:584969/1‑271 (MQ=255)
                                                                                                                                                                                                                                     ccAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTcccc                                                                                                                                             >  4:172561/1‑149 (MQ=255)
                                                                                                                                                                                                                                     ccAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCCCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACca                 <  2:693147/273‑1 (MQ=255)
                                                                                                                                                                                                                                         gCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCccacca                                                                                                                                         <  2:172562/149‑1 (MQ=255)
                                                                                                                                                                                                                                         acACCGAACCGCAGAACTGAGAGAAAAGAAGCCAGAATCAGATGGCGGCGGAAGGGAAACTCTCTAACACATGTCATTATTGTTGGGCCAAACCATGCGCCAAGGCGGTAAAAAGGCAAATTGCGGATCACCCTGTTCATTACCCACCAGGGCGTGGGAACGAAGTGAAACTCCACACTCAGCAGTACAAACAAGTCGAATAGTGGCAGCGCCAGACGGCTAAAAAGTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGc            <  4:292745/273‑1 (MQ=255)
                                                                                                                                                                                                                                               accgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCCCAACGGCCCCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCCCACTGTTCATTCCCCACCACGGCGCGGGAACGACATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACgg                                                                   >  1:78357/1‑213 (MQ=255)
                                                                                                                                                                                                                                                ccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGt                                                                                                                                   <  6:756496/148‑1 (MQ=255)
                                                                                                                                                                                                                                                 cgCCAGAACGGTAAGCAAAGAAGCCATCAGCGGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGCCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGc  <  2:84255/275‑1 (MQ=255)
                                                                                                                                                                                                                                                 cgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTgg                                                                                                                                 <  8:127050/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 cgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTgg                                                                                                                                 <  7:857016/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 cgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTg                                                                                                                                  <  6:374958/148‑1 (MQ=255)
                                                                                                                                                                                                                                                 cgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTg                                                                                                                                  <  5:336760/148‑1 (MQ=255)
                                                                                                                                                                                                                                                 cgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTcag                                                                                                        <  4:562047/174‑1 (MQ=255)
                                                                                                                                                                                                                                                   ccAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACCAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATAACCAGAAATAGAAAATTAACGCATACAGATACCAAAAGCt           >  2:1016/1‑265 (MQ=255)
                                                                                                                                                                                                                                                    cAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGaa                                                                                                                              <  8:231908/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    cAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGaa                                                                                                                              <  8:260275/149‑1 (MQ=255)
                                                                                                                                                                                                                                                         cGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGaaa                                                                                                                         <  7:686029/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           gTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAAt                                                                                                                        >  5:251189/1‑148 (MQ=255)
                                                                                                                                                                                                                                                             aaGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTaa                                                                                                                     >  7:658707/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  aagaagCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGc                                                                                                                 <  6:436173/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                      cgCCATCTGCTGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGAt                                                                                                                                                              >  7:183050/2‑99 (MQ=255)
                                                                                                                                                                                                                                                                      cgCCATCTGCTGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGAt                                                                                                                                                              <  8:183050/98‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gCCATCTGCTGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCacac                                                                                                            >  6:487224/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                          aTCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTggg                                                                                                                                >  8:197657/1‑125 (MQ=255)
                                                                                                                                                                                                                                                                          aTCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTggg                                                                                                                                <  7:197657/125‑1 (MQ=255)
                                                                                                                                                                                                                                                                          aTCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAAcc                                                                                                                   >  6:122191/1‑138 (MQ=255)
                                                                                                                                                                                                                                                                            cagcagATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTcagc                                                                                                       >  6:597009/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                             agcagATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTcagcag                                                                                                     <  8:455771/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gcagATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTcagcag                                                                                                     <  5:547149/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gcagATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTcagcag                                                                                                     >  6:820632/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
GTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAGAATGCGATGGGTGGTGTAGATAGCAATGGTGTTGGAGCCAATCACATTCAGCAGGCTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGC  >  NC_000913/1174232‑1174748

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