Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,336,599 T→C 6.5% L141P (CTT→CCT)  pyrF → orotidine‑5'‑phosphate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,336,5990TC 6.5% 517.4 / 21.6 168L141P (CTT→CCT) pyrForotidine‑5'‑phosphate decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (68/89);  new base C (3/8);  total (71/97)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 171 total aligned reads displayed.

GCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTACGCCGGG  >  CP009273/1336452‑1336746
                                                                                                                                                   |                                                                                                                                                   
     tgGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCa                                                                                                                                                <  2:1902325/148‑1 (MQ=255)
       gATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATg                                                                                                                                              <  1:2202840/148‑1 (MQ=255)
       gATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATg                                                                                                                                              >  2:1586146/1‑148 (MQ=255)
       gATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATg                                                                                                                                              >  2:1586305/1‑148 (MQ=255)
        aTGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGAGGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGa                                                                                                                                             <  1:955400/148‑1 (MQ=255)
          ggTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGaca                                                                                                                                           <  1:2348199/148‑1 (MQ=255)
               aTGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTc                                                                                                                                      <  2:1958136/148‑1 (MQ=255)
                 gTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCAc                                                                                                                                    <  1:2289637/148‑1 (MQ=255)
                   tCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCt                                                                                                                                  >  2:838918/1‑148 (MQ=255)
                         ctctGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGAt                                                                                                                            <  1:2428864/148‑1 (MQ=255)
                          tctGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGAtt                                                                                                                           <  2:83171/148‑1 (MQ=255)
                                ggggCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAg                                                                                                                     <  2:1096032/148‑1 (MQ=255)
                                 gggCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGAATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGa                                                                                                                    <  1:664388/148‑1 (MQ=255)
                                      cgTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAAAGTc                                                                                                               >  2:696404/1‑148 (MQ=255)
                                           gatgaCCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTggcg                                                                                                          >  1:877465/1‑148 (MQ=255)
                                            atgaCCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTggcgg                                                                                                         >  2:2901313/1‑148 (MQ=255)
                                                   cAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGAc                                                                                                  >  2:2085805/1‑148 (MQ=255)
                                                   cAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGAc                                                                                                  >  1:3072126/1‑148 (MQ=255)
                                                   cAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGAc                                                                                                  >  1:1029963/1‑148 (MQ=255)
                                                      cgcgTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCa                                                                                               >  2:2844877/1‑148 (MQ=255)
                                                         gTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCaaaa                                                                                            <  1:111896/148‑1 (MQ=255)
                                                         gTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCaaaa                                                                                            <  2:1604281/148‑1 (MQ=255)
                                                            aGGCACTGGGTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAAtg                                                                                         <  1:1805710/148‑1 (MQ=255)
                                                              gCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAAtgtg                                                                                       >  2:262168/1‑148 (MQ=255)
                                                              gCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAAtgtg                                                                                       <  1:1976355/148‑1 (MQ=255)
                                                                aCTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGc                                                                                     <  2:188776/148‑1 (MQ=255)
                                                                 cTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGCTTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGcc                                                                                    >  2:1120073/1‑148 (MQ=255)
                                                                 cTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGcc                                                                                    <  1:3021592/148‑1 (MQ=255)
                                                                 cTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGcc                                                                                    >  2:308891/1‑148 (MQ=255)
                                                                  tGGTACCGTTTGGCAAAGATGCACCGCTTTTGATGGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCt                                                                                   <  1:865956/148‑1 (MQ=255)
                                                                  tGGATCCGTTTGGCAAAGATGCACCGCTTTTGAATGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTAGATCTTCGCATGACAATCTCACCTGCAGATTATTCAGAAACTCTGTCGGCACTGACCCAAAATTGTTGcca                                                                                   >  2:2409959/1‑147 (MQ=255)
                                                                   ggTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCtt                                                                                  >  1:833666/1‑148 (MQ=255)
                                                                    gTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCCGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTg                                                                                 <  2:2033564/148‑1 (MQ=255)
                                                                    gTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCCGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTg                                                                                 <  2:2033566/148‑1 (MQ=255)
                                                                     ttCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGa                                                                                <  2:938244/148‑1 (MQ=255)
                                                                     ttCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGa                                                                                <  2:362345/148‑1 (MQ=255)
                                                                     ttCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGa                                                                                >  1:2908661/1‑148 (MQ=255)
                                                                       ccGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAtg                                                                              <  2:2437140/148‑1 (MQ=255)
                                                                         gTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAtggt                                                                            <  1:472745/148‑1 (MQ=255)
                                                                          tttGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAtggtg                                                                           >  1:1027054/1‑148 (MQ=255)
                                                                            tGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAtggtgtg                                                                         <  1:838996/148‑1 (MQ=255)
                                                                               cAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAtggtgtggtg                                                                      <  2:3272881/148‑1 (MQ=255)
                                                                                aaaGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCGTGATggtgtggtgt                                                                     >  1:967609/1‑148 (MQ=255)
                                                                                 aaGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGgtgtggtgtg                                                                    >  2:3290053/1‑148 (MQ=255)
                                                                                  aGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGgtgtgt                                                                   >  2:2228294/1‑148 (MQ=255)
                                                                                  aGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGgtgtgt                                                                   <  2:335671/148‑1 (MQ=255)
                                                                                  aGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGgtgtgt                                                                   <  2:1545814/148‑1 (MQ=255)
                                                                                   gATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGtt                                                                  <  1:3274116/148‑1 (MQ=255)
                                                                                   gATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGtt                                                                  <  1:3198507/148‑1 (MQ=255)
                                                                                   gATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGtt                                                                  <  2:384332/148‑1 (MQ=255)
                                                                                   gATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGtt                                                                  <  2:565266/148‑1 (MQ=255)
                                                                                    aTGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTc                                                                 >  1:2715967/1‑148 (MQ=255)
                                                                                    aTGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTc                                                                 <  2:611708/148‑1 (MQ=255)
                                                                                      gCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTg                                                               <  2:2467872/148‑1 (MQ=255)
                                                                                          cGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCa                                                           >  1:2870633/1‑148 (MQ=255)
                                                                                             ttttGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGa                                                        <  2:888410/148‑1 (MQ=255)
                                                                                               ttGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAg                                                      <  2:1254484/148‑1 (MQ=255)
                                                                                               ttGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAg                                                      >  2:3303511/1‑148 (MQ=255)
                                                                                                  aTTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCtg                                                   <  2:703601/148‑1 (MQ=255)
                                                                                                  aTTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCtg                                                   <  1:2844863/148‑1 (MQ=255)
                                                                                                  aTTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCtg                                                   <  1:1660012/148‑1 (MQ=255)
                                                                                                   ttGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCtgt                                                  <  2:49920/148‑1 (MQ=255)
                                                                                                   ttGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCtgt                                                  >  1:2805345/1‑148 (MQ=255)
                                                                                                    tGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTCGATGGTGTGGTGTGTTCTGCTCAGGAAGCtgtg                                                 >  2:2779880/1‑148 (MQ=255)
                                                                                                      cTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCGTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTgcg                                               >  1:1524935/1‑148 (MQ=255)
                                                                                                       tgtgACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTgcgc                                              <  2:224792/148‑1 (MQ=255)
                                                                                                        gtgACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCt                                             >  1:1745618/1‑148 (MQ=255)
                                                                                                         tgACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCtt                                            >  2:1196801/1‑148 (MQ=255)
                                                                                                            cAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTaa                                         <  2:2672153/148‑1 (MQ=255)
                                                                                                             aGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTaaa                                        <  1:545738/148‑1 (MQ=255)
                                                                                                               tgtTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACa                                      >  2:2127320/1‑148 (MQ=255)
                                                                                                                gtTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAg                                     <  1:320504/148‑1 (MQ=255)
                                                                                                                 ttGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAgg                                    <  1:2901313/148‑1 (MQ=255)
                                                                                                                 ttGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAgg                                    <  1:1217395/148‑1 (MQ=255)
                                                                                                                    aCCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTAt                                 <  2:3229805/148‑1 (MQ=255)
                                                                                                                     ccAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTAtt                                <  1:2517025/148‑1 (MQ=255)
                                                                                                                     ccAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTAtt                                <  2:1119314/148‑1 (MQ=255)
                                                                                                                      cAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTc                               >  2:1405139/1‑148 (MQ=255)
                                                                                                                       aGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGGAATCg                              >  2:122543/1‑148 (MQ=255)
                                                                                                                          aTGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTc                           <  1:2561017/148‑1 (MQ=255)
                                                                                                                            ggAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAg                         <  1:826521/148‑1 (MQ=255)
                                                                                                                              aaGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGa                       >  1:2893942/1‑148 (MQ=255)
                                                                                                                               aGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAg                      >  2:443491/1‑148 (MQ=255)
                                                                                                                                gCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGt                     >  2:629797/1‑148 (MQ=255)
                                                                                                                                 ccAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGtt                    <  1:3155072/148‑1 (MQ=255)
                                                                                                                                   aGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCa                  >  2:2593084/1‑148 (MQ=255)
                                                                                                                                     cGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATTGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCaaa                <  1:700882/148‑1 (MQ=255)
                                                                                                                                     cGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCaaa                >  2:1670092/1‑148 (MQ=255)
                                                                                                                                        ccTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTg             >  1:2981538/1‑148 (MQ=255)
                                                                                                                                          tGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGt           <  2:213748/148‑1 (MQ=255)
                                                                                                                                            gTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAGTGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTa         >  1:1172410/1‑148 (MQ=255)
                                                                                                                                             tCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTAc        >  1:2568673/1‑148 (MQ=255)
                                                                                                                                             tCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTAc        >  1:2983032/1‑148 (MQ=255)
                                                                                                                                             tCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTAc        >  2:1490181/1‑148 (MQ=255)
                                                                                                                                             tCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTAc        >  2:136940/1‑148 (MQ=255)
                                                                                                                                             tCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGAGGAACAGACGCAAAAATGTGGCCTAGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAAAAAGTATTAGGTCAGGAGATCAAACTGTTTAc        >  2:1490328/1‑148 (MQ=255)
                                                                                                                                               gATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTACGc      <  1:1586146/148‑1 (MQ=255)
                                                                                                                                               gATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTACGc      <  1:1586305/148‑1 (MQ=255)
                                                                                                                                                 tCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTACGCCg    >  2:1029580/1‑148 (MQ=255)
                                                                                                                                                   ttGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTTGTGTGTTCTGCTCAGGAAGCTGTGCGATTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTACGCCggg  <  1:21321/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTACGCCGGG  >  CP009273/1336452‑1336746

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