Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 420,946 T→A 6.0% intergenic (+36/‑40) brnQ → / → proY branched‑chain amino acid transport system 2 carrier protein; LIV‑II transport system for Ile, Leu, and Val/proline‑specific permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913420,9460TA 6.0% 393.6 / 9.5 114intergenic (+36/‑40)brnQ/proYbranched‑chain amino acid transport system 2 carrier protein; LIV‑II transport system for Ile, Leu, and Val/proline‑specific permease
Reads supporting (aligned to +/- strand):  ref base T (50/57);  new base A (4/3);  total (54/60)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

ACCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGC  >  NC_000913/420799‑421092
                                                                                                                                                   |                                                                                                                                                  
aCCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGAt                                                                                                                                                    <  4:843723/148‑1 (MQ=255)
 ccGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGgctgcgctgc                                                                                                                                                  <  1:925917/149‑1 (MQ=255)
      ggCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGt                                                                                                                                             <  2:1074541/149‑1 (MQ=255)
       gCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTTACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCATCGCAGCCCCGTg                                                                                                                                            <  2:1005286/149‑1 (MQ=255)
       gCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGt                                                                                                                                             <  3:738820/148‑1 (MQ=255)
        ccGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCATCGCAGCCCCGTgg                                                                                                                                           <  2:347947/149‑1 (MQ=255)
         cGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGt                                                                                                                                          >  1:681347/1‑149 (MQ=255)
          gAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGAt                                                                                                                                                    <  3:582519/138‑1 (MQ=255)
          gAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGAt                                                                                                                                                    >  4:582519/1‑138 (MQ=255)
                ggTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTa                                                                                                                                    >  3:182802/1‑148 (MQ=255)
                 gTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTAtt                                                                                                                                  <  2:70108/149‑1 (MQ=255)
                 gTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTAtt                                                                                                                                  >  2:470251/1‑149 (MQ=255)
                 ctCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCGATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTAt                                                                                                                                   <  3:94671/147‑1 (MQ=255)
                  tCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTAtt                                                                                                                                  >  3:669029/1‑148 (MQ=255)
                    tGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATtgtg                                                                                                                               <  2:95786/149‑1 (MQ=255)
                         tGGTTAATGCTAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGa                                                                                                                           <  4:545550/148‑1 (MQ=255)
                         tGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGAt                                                                                                                          <  2:282044/149‑1 (MQ=255)
                           gTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATgg                                                                                                                        >  2:277401/1‑149 (MQ=255)
                            ttAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATgg                                                                                                                        >  4:105015/1‑148 (MQ=255)
                            ttAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTAATTGTGTTGATggg                                                                                                                       >  2:904437/1‑149 (MQ=255)
                             tAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATggg                                                                                                                       >  3:854095/1‑148 (MQ=255)
                               aTGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGtt                                                                                                                     >  3:470713/1‑148 (MQ=255)
                                 gCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAgg                                                                                                                  >  2:310565/1‑149 (MQ=255)
                                  ccAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTgtgt                                                                                                                              >  1:732263/1‑136 (MQ=255)
                                    aaCAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTTACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAAt                                                                                                               >  2:796884/1‑149 (MQ=255)
                                       aGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTg                                                                                                             <  4:34756/148‑1 (MQ=255)
                                       aGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTg                                                                                                             <  3:100494/148‑1 (MQ=255)
                                       aGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGa                                                                                                            <  2:779663/149‑1 (MQ=255)
                                       aGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGa                                                                                                            <  1:904437/149‑1 (MQ=255)
                                       aGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGa                                                                                                            <  2:466789/149‑1 (MQ=255)
                                          ggtgATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATgg                                                                                                         <  2:433740/149‑1 (MQ=255)
                                            tgATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGa                                                                                                        <  4:12712/148‑1 (MQ=255)
                                             gATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGaa                                                                                                       <  3:105015/148‑1 (MQ=255)
                                                  tggtTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTaa                                                                                                  <  3:379811/148‑1 (MQ=255)
                                                  tggtTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAg                                                                                                 <  1:277401/149‑1 (MQ=255)
                                                   ggtTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGa                                                                                                <  1:470251/149‑1 (MQ=255)
                                                     ttCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAAc                                                                                              >  2:600688/1‑149 (MQ=255)
                                                       cTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACa                                                                                             <  3:512149/148‑1 (MQ=255)
                                                         ggCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAg                                                                                           <  4:669029/148‑1 (MQ=255)
                                                         ggCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAg                                                                                           >  3:115556/1‑148 (MQ=255)
                                                         ggCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGc                                                                                          >  1:1001210/1‑149 (MQ=255)
                                                         ggCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGc                                                                                          >  1:554629/1‑149 (MQ=255)
                                                           ccATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGa                                                                                                            <  3:177887/129‑1 (MQ=255)
                                                           ccATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGa                                                                                                            >  4:177887/1‑129 (MQ=255)
                                                                 tCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAAtt                                                                                                              >  2:220566/1‑121 (MQ=255)
                                                                 tCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAAtt                                                                                                              <  1:220566/121‑1 (MQ=255)
                                                                 tCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGt                                                                                  >  2:987499/1‑149 (MQ=255)
                                                                  cTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTg                                                                                                             <  4:67102/121‑1 (MQ=255)
                                                                  cTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTg                                                                                                             >  3:67102/1‑121 (MQ=255)
                                                                  cTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCATCGCAGCCCCGTgg                                                                                                                                           >  1:225723/1‑91 (MQ=38)
                                                                   tGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTgg                                                                                <  1:713387/149‑1 (MQ=255)
                                                                     ggATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTggg                                                                               >  3:1011051/1‑148 (MQ=255)
                                                                      gATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGc                                                                              >  3:931233/1‑148 (MQ=255)
                                                                       aTCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTa                                                                            >  1:591238/1‑149 (MQ=255)
                                                                          gcgcgGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTa                                                                                                                    >  2:124186/1‑106 (MQ=255)
                                                                          gcgcgGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTa                                                                                                                    <  1:124186/106‑1 (MQ=255)
                                                                               gCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCg                                                                    <  1:519242/149‑1 (MQ=255)
                                                                                  ggTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGcc                                                                  >  3:115598/1‑148 (MQ=255)
                                                                                     cgtcAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCcaca                                                               >  4:10019/1‑148 (MQ=255)
                                                                                      gtcAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTAtt                                                                                                                                  >  1:1001739/1‑80 (MQ=255)
                                                                                      gtcAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTAtt                                                                                                                                  <  2:1001739/80‑1 (MQ=255)
                                                                                            tGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCtt                                                        >  3:30668/1‑148 (MQ=255)
                                                                                             gACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCttt                                                       <  4:115556/148‑1 (MQ=255)
                                                                                             gACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTa                                                      <  1:768298/149‑1 (MQ=255)
                                                                                             gACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGAACCCGCCACATACGCTTTa                                                      <  2:1001210/149‑1 (MQ=255)
                                                                                                cTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATg                                                    >  4:284091/1‑148 (MQ=255)
                                                                                                 tCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGc                                                                              >  3:117875/1‑121 (MQ=255)
                                                                                                 tCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGc                                                                              <  4:117875/121‑1 (MQ=255)
                                                                                                 tCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGc                                                  >  2:655114/1‑149 (MQ=255)
                                                                                                   cAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTg                                                                                                             <  3:773840/88‑1 (MQ=255)
                                                                                                   cAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTg                                                                                                             >  4:773840/1‑88 (MQ=255)
                                                                                                    aGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCAc                                                <  4:30668/148‑1 (MQ=255)
                                                                                                         tCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGtt                                          >  2:469685/1‑149 (MQ=255)
                                                                                                          cACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCt                                                         <  4:797122/133‑1 (MQ=255)
                                                                                                           aCTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCa                                        >  1:826027/1‑149 (MQ=255)
                                                                                                            cTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATa                                                             >  4:453197/1‑127 (MQ=255)
                                                                                                            cTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATa                                                             <  3:453197/127‑1 (MQ=255)
                                                                                                            cTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCa                                        <  4:570934/148‑1 (MQ=255)
                                                                                                                 tCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGaaa                                                                                                      <  1:820640/81‑1 (MQ=255)
                                                                                                                 tCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGaaa                                                                                                      >  2:820640/1‑81 (MQ=255)
                                                                                                                  cACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAAtt                                  >  3:104250/1‑148 (MQ=255)
                                                                                                                     tGAACATTTGTTTTAACCACGGGGCTTCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTa                                                      <  2:649889/125‑1 (MQ=255)
                                                                                                                     tGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTgg                                                                                <  1:664933/99‑1 (MQ=255)
                                                                                                                     tGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTa                                                      >  1:649889/1‑125 (MQ=255)
                                                                                                                     tGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGc                               >  4:282695/1‑148 (MQ=255)
                                                                                                                     tGAACATTTGTTTTAACCACGCGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTgg                                                                                >  2:664933/1‑99 (MQ=255)
                                                                                                                         cATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCg                           <  4:104250/148‑1 (MQ=255)
                                                                                                                                tttAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGa                                                                                                                 >  4:819289/1‑55 (MQ=255)
                                                                                                                                tttAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGa                                                                                                                 <  3:819289/55‑1 (MQ=255)
                                                                                                                                tttAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGttt                   <  2:759784/149‑1 (MQ=255)
                                                                                                                                tttAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGtt                    >  4:87927/1‑148 (MQ=255)
                                                                                                                                  tAACCACGGGGCTGCGATGCAACGTGGTTTTTCACTGTGTTGATTTGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTTCCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTCGCATCGGACTGTTTTa                 <  2:675115/149‑1 (MQ=255)
                                                                                                                                   aaCCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTAc                <  1:600688/149‑1 (MQ=255)
                                                                                                                                   aaCCACGGGGCATCGCAGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTAc                >  2:133981/1‑149 (MQ=255)
                                                                                                                                   aaCCACCGCGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTa                 <  4:235899/148‑1 (MQ=255)
                                                                                                                                     ccACGGGGCATCGCAGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACgg              <  1:409003/149‑1 (MQ=255)
                                                                                                                                      cACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGt             >  1:778376/1‑149 (MQ=255)
                                                                                                                                      cACGGGGCATCGCAGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGt             >  2:706188/1‑149 (MQ=255)
                                                                                                                                               gCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGAc    >  2:741656/1‑149 (MQ=255)
                                                                                                                                                  aTGCCCTGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGc  >  3:724734/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGC  >  NC_000913/420799‑421092

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