Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,336,599 T→C 21.2% 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,5990TC21.2% 391.4 / 115.9 184L141P (CTT→CCT) pyrForotidine‑5'‑phosphate decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (60/85);  new base C (15/24);  total (75/109)
Fisher's exact test for biased strand distribution p-value = 8.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

Only 100 of 185 total aligned reads displayed.

CGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTACGCCGG  >  CP009273/1336453‑1336745
                                                                                                                                                  |                                                                                                                                                  
   gtgGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGTAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGc                                                                                                                                                >  1:2106624/1‑148 (MQ=255)
   gtgGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGc                                                                                                                                                <  2:1852859/148‑1 (MQ=255)
      gATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATg                                                                                                                                             <  1:275850/148‑1 (MQ=255)
      gATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATg                                                                                                                                             >  1:3368324/1‑148 (MQ=255)
         ggTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGaca                                                                                                                                          <  2:867083/148‑1 (MQ=255)
         ggTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGaca                                                                                                                                          <  1:2169544/148‑1 (MQ=255)
         ggTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGaca                                                                                                                                          <  1:3061175/148‑1 (MQ=255)
         ggTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGaca                                                                                                                                          <  1:869856/148‑1 (MQ=255)
          gTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGacac                                                                                                                                         >  2:164185/1‑148 (MQ=255)
            gAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTg                                                                                                                                       <  2:576876/148‑1 (MQ=255)
            gAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTg                                                                                                                                       <  2:576792/148‑1 (MQ=255)
             aaTGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGt                                                                                                                                      <  2:782511/148‑1 (MQ=255)
              aTGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTc                                                                                                                                     <  2:2154557/148‑1 (MQ=255)
                gTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCAc                                                                                                                                   <  2:3096494/148‑1 (MQ=255)
                     tGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCa                                                                                                                              >  1:3557585/1‑148 (MQ=255)
                       ccTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTTCAGa                                                                                                                            >  1:204960/1‑148 (MQ=255)
                       ccTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGa                                                                                                                            <  2:553991/148‑1 (MQ=255)
                          ctGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTa                                                                                                                         <  2:1341282/148‑1 (MQ=255)
                           tggtggGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTAt                                                                                                                        <  1:113554/148‑1 (MQ=255)
                             gtggGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCCCCTGCAGATTATGc                                                                                                                      >  2:46247/1‑148 (MQ=255)
                                    gcgTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGGCCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGt                                                                                                               <  1:830995/148‑1 (MQ=255)
                                    gcgTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGt                                                                                                               <  1:2033329/148‑1 (MQ=255)
                                       tATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTg                                                                                                            >  2:1292743/1‑148 (MQ=255)
                                       tATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTg                                                                                                            >  2:450150/1‑148 (MQ=255)
                                       tATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTg                                                                                                            >  1:3194840/1‑148 (MQ=255)
                                        atgatgACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTgg                                                                                                           >  2:3138226/1‑148 (MQ=255)
                                        atgatgACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTgg                                                                                                           <  2:3262757/148‑1 (MQ=255)
                                          gatgaCCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTggcg                                                                                                         >  2:3201674/1‑148 (MQ=255)
                                          gatgaCCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTggcg                                                                                                         >  1:2202771/1‑148 (MQ=255)
                                               ccGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACt                                                                                                    <  2:2232451/148‑1 (MQ=255)
                                                cGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTg                                                                                                   <  1:1230040/148‑1 (MQ=255)
                                                    gcgcGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATTCAGAACGTATGGCGGCACTGACGc                                                                                               >  1:1013097/1‑148 (MQ=255)
                                                    gcgcGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGc                                                                                               >  2:1130254/1‑148 (MQ=255)
                                                    gcgcGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGAAGAACGTCTGGCGGCACTGACGc                                                                                               >  1:1013096/1‑148 (MQ=255)
                                                     cgcgTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCa                                                                                              >  2:3706801/1‑148 (MQ=255)
                                                           aGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAAtg                                                                                        >  2:2716509/1‑148 (MQ=255)
                                                           aGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAAtg                                                                                        <  2:1519342/148‑1 (MQ=255)
                                                           aGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAAtg                                                                                        <  2:1049578/148‑1 (MQ=255)
                                                            ggCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAAtgt                                                                                       <  2:713677/148‑1 (MQ=255)
                                                            ggCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAAtgt                                                                                       <  2:2180337/148‑1 (MQ=255)
                                                            ggCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAAtgt                                                                                       <  2:1792632/148‑1 (MQ=255)
                                                             gCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAAtgtg                                                                                      >  2:2312030/1‑148 (MQ=255)
                                                                cTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGcc                                                                                   >  1:1810496/1‑148 (MQ=255)
                                                                  ggTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCtt                                                                                 >  1:564103/1‑148 (MQ=255)
                                                                     tCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAt                                                                              <  2:1385262/148‑1 (MQ=255)
                                                                     tCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAt                                                                              <  2:3287059/148‑1 (MQ=255)
                                                                        gTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAtggt                                                                           <  1:1423679/148‑1 (MQ=255)
                                                                          ttGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAtggtgt                                                                         <  1:2407592/148‑1 (MQ=255)
                                                                            ggCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAtggtgtgg                                                                       <  1:2816450/148‑1 (MQ=255)
                                                                              cAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAtggtgtggtg                                                                     <  2:2832074/148‑1 (MQ=255)
                                                                              cAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAtggtgtggtg                                                                     >  2:1949078/1‑148 (MQ=255)
                                                                                aaGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGgtgtggtgtg                                                                   >  1:1593452/1‑148 (MQ=255)
                                                                                aaGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGgtgtggtgtg                                                                   <  1:1970109/148‑1 (MQ=255)
                                                                                aaGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGgtgtggtgtg                                                                   >  1:724655/1‑148 (MQ=255)
                                                                                 aGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGgtgtgt                                                                  >  1:837634/1‑148 (MQ=255)
                                                                                  gATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGtt                                                                 <  2:3555955/148‑1 (MQ=255)
                                                                                  gATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGtt                                                                 >  1:1453515/1‑148 (MQ=255)
                                                                                       aCCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCt                                                            <  1:3622721/148‑1 (MQ=255)
                                                                                        ccGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTc                                                           <  1:119198/148‑1 (MQ=255)
                                                                                        ccGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTc                                                           >  1:1141828/1‑148 (MQ=255)
                                                                                           cTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAgg                                                        >  2:2471702/1‑148 (MQ=255)
                                                                                           cTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAgg                                                        >  2:2471927/1‑148 (MQ=255)
                                                                                            ttttGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGa                                                       >  2:2803780/1‑148 (MQ=255)
                                                                                              ttGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGAAACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAg                                                     >  1:3533504/1‑148 (MQ=255)
                                                                                                gATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCt                                                   <  1:1539596/148‑1 (MQ=255)
                                                                                                 aTTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCtg                                                  <  2:1542424/148‑1 (MQ=255)
                                                                                                  ttGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCtgt                                                 >  1:2495786/1‑148 (MQ=255)
                                                                                                   tGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGTGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCtgtg                                                >  2:1238420/1‑148 (MQ=255)
                                                                                                   tGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCtgtg                                                <  1:1501293/148‑1 (MQ=255)
                                                                                                   tGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCtgtg                                                <  2:303248/148‑1 (MQ=255)
                                                                                                     cTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTgcg                                              <  1:1976641/148‑1 (MQ=255)
                                                                                                       gtgACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCt                                            <  1:3333085/148‑1 (MQ=255)
                                                                                                       gtgACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCt                                            <  1:977682/148‑1 (MQ=255)
                                                                                                           cAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTaa                                        <  2:3634084/148‑1 (MQ=255)
                                                                                                             gtgtTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAAc                                      <  1:2889191/148‑1 (MQ=255)
                                                                                                               gtTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAg                                    >  2:3395453/1‑148 (MQ=255)
                                                                                                                ttGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAgg                                   <  2:2715472/148‑1 (MQ=255)
                                                                                                                ttGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAgg                                   <  1:1780563/148‑1 (MQ=255)
                                                                                                                ttGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAgg                                   <  1:1780421/148‑1 (MQ=255)
                                                                                                                 tGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGt                                  <  2:3048235/148‑1 (MQ=255)
                                                                                                                   aCCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCGTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTAt                                >  2:2287068/1‑148 (MQ=255)
                                                                                                                   aCCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTAt                                >  2:3157836/1‑148 (MQ=255)
                                                                                                                      aGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCg                             <  1:152824/148‑1 (MQ=255)
                                                                                                                          tGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCa                         <  2:2305645/148‑1 (MQ=255)
                                                                                                                          tGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCa                         >  2:1932071/1‑148 (MQ=255)
                                                                                                                          tGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCa                         <  2:2929647/148‑1 (MQ=255)
                                                                                                                             aaGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGa                      >  1:1425051/1‑148 (MQ=255)
                                                                                                                                    cGACCTGGTCGCTCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCaaa               <  2:2546302/148‑1 (MQ=255)
                                                                                                                                    cGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCaaa               <  2:1593452/148‑1 (MQ=255)
                                                                                                                                    cGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCaaa               <  1:2061289/148‑1 (MQ=255)
                                                                                                                                    cGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCaaa               <  2:629706/148‑1 (MQ=255)
                                                                                                                                    cGACCTGGTCGATCCTGGCAAGACACTGTCAACTGCAGATTATGCAGAACGTCAGGCGGCACTGAAGCAAAAATGTGGCCTTGATGGAGTGGTGTGTTCTTCTCAGGAAGCTGTGCGCTTTAAACAGGTATACGGTAAGGAGTTCaaa               >  2:676193/1‑148 (MQ=255)
                                                                                                                                          ggTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGtt         <  2:26438/148‑1 (MQ=255)
                                                                                                                                             cGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTACg      >  2:1987560/1‑148 (MQ=255)
                                                                                                                                             cGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTACg      >  2:1442463/1‑148 (MQ=255)
                                                                                                                                               aTCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTACGcc    >  1:2274660/1‑148 (MQ=255)
                                                                                                                                                tCTTGGCATGACACTGTCACCTTCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTACGCCg   <  2:1214441/148‑1 (MQ=255)
                                                                                                                                                tCTTGGCATGACACTGTCACCTTCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGCGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTACGCCg   <  2:204960/148‑1 (MQ=255)
                                                                                                                                                tCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTACGCCg   <  1:2505005/148‑1 (MQ=255)
                                                                                                                                                 cTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTACGCCgg  >  2:1570713/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTACGCCGG  >  CP009273/1336453‑1336745

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: