Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 420,945 A→C 5.4% intergenic (+35/‑41) 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,9450AC 5.4% 388.2 / 11.0 129intergenic (+35/‑41)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 A (58/64);  new base C (3/4);  total (61/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.73e-01

Only 100 of 129 total aligned reads displayed.

CCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCC  >  NC_000913/420800‑421093
                                                                                                                                                 |                                                                                                                                                    
ccGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCCGCGCTCACTAAATCACTGAACATTTGTTTTAACAACGGGGCTGCGATGcc                                                                                                                                                  >  1:1695090/1‑150 (MQ=255)
  gCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCAGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGcccc                                                                                                                                                >  2:1451379/1‑150 (MQ=255)
    tGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGACATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGt                                                                                                                                              <  2:1659566/150‑1 (MQ=255)
     ggCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTg                                                                                                                                             <  2:1397240/150‑1 (MQ=255)
     ggCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCATCGCAGCCCCGTg                                                                                                                                             <  1:1008677/150‑1 (MQ=255)
     ggCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCATCGCAGCCCCGTg                                                                                                                                             <  1:1711684/150‑1 (MQ=255)
     ggCCGAACAAGGGCTGGCGGGGTTAATGCCAACAGTGGTGATGGTGGGTCTGGCCATTATCTGGGTTCGCGCGGCAGGTCGTCAGGGGTCCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTg                                                                                                                                             <  2:2583/150‑1 (MQ=255)
      gCCGACCAAGGTCTGGCGCGGTTAATGCCAACAGTGGTGAGGGTGGATCTGGCCATTGTCTGGGAGCGCGGGGCAGGTCGCCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTgg                                                                                                                                            <  2:326816/150‑1 (MQ=255)
       ccGATCAAGGTCTGGCGTGGTTTAAGCCAACAGGGGTGATGGTGGTTGTGGCCATTATCTGGGATCGCGCGGCAGGGCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGt                                                                                                                                           <  2:1067808/150‑1 (MQ=255)
       ccGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCATCGCAGCCCCGTGGt                                                                                                                                           <  2:696277/150‑1 (MQ=255)
           aCAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGttttt                                                                                                                                       <  2:1100169/150‑1 (MQ=255)
            cAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGtttttt                                                                                                                                      >  1:1154765/1‑150 (MQ=255)
                gTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATtg                                                                                                                                  <  1:968074/150‑1 (MQ=255)
                 tCTGGCGAGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATtgt                                                                                                                                 >  2:860035/1‑150 (MQ=255)
                  cTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGTCTGCGATTCCCCGTGGTTTTTTATtttg                                                                                                                                >  2:133123/1‑150 (MQ=255)
                   tGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTgtgt                                                                                                                               >  1:687128/1‑150 (MQ=255)
                      cGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGa                                                                                                                            >  2:1093875/1‑150 (MQ=255)
                        tGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATg                                                                                                                          <  2:133570/150‑1 (MQ=255)
                         ggTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATgg                                                                                                                         >  1:147133/1‑150 (MQ=255)
                          gTTAATGGCAACAGTGGTGATGGGGGGTCTGGCCATTATCTGGGATCGCGCGGCAGGGCGTCAGGTGACCTCCAGCGCTCACTAAATCCCTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATggg                                                                                                                        <  2:131603/150‑1 (MQ=255)
                             actGCAAACAGGGGTGAGGGGGGGTCTGGCCATTATCTGGGATCGCGCGGCAGGGCGGCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTa                                                                                                                     <  2:889162/148‑1 (MQ=255)
                             aaTGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTa                                                                                                                     >  1:978930/1‑150 (MQ=255)
                              aTGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAg                                                                                                                    <  1:1470350/150‑1 (MQ=255)
                               tGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAgg                                                                                                                   <  1:309468/150‑1 (MQ=255)
                                gCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGa                                                                                                                  <  2:1378552/150‑1 (MQ=255)
                                     cAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGa                                                                                                             <  2:486267/150‑1 (MQ=255)
                                       gtggtgATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCTCTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATg                                                                                                           <  1:875953/150‑1 (MQ=255)
                                       gtggtgATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAAt                                                                                                                <  1:38851/145‑1 (MQ=255)
                                         ggtgATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGa                                                                                                         >  2:205647/1‑150 (MQ=255)
                                         ggtgATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGa                                                                                                         >  2:1372046/1‑150 (MQ=255)
                                              tggtggTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTa                                                                                                    <  2:1466707/150‑1 (MQ=255)
                                                 tgggTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGATTTGATGGAAAGTAAGa                                                                                                 <  1:1048157/150‑1 (MQ=255)
                                                  ggtTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGaa                                                                                                >  1:1523407/1‑150 (MQ=255)
                                                     tCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACaa                                                                                             <  1:269943/150‑1 (MQ=255)
                                                     tCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACaa                                                                                             <  2:639476/150‑1 (MQ=255)
                                                      cTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAg                                                                                            <  2:794196/150‑1 (MQ=255)
                                                       tGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTCCGATGCCGCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAga                                                                                           >  2:1079037/1‑149 (MQ=255)
                                                              tATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCg                                                                                    <  2:1154765/150‑1 (MQ=255)
                                                              tATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCg                                                                                    >  2:459036/1‑150 (MQ=255)
                                                               aTCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGt                                                                                   >  1:792374/1‑150 (MQ=255)
                                                               aTCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGt                                                                                   >  2:569851/1‑150 (MQ=255)
                                                                tCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTg                                                                                  >  2:834738/1‑150 (MQ=255)
                                                                tCTGGGATCGCGCGGAAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACTTTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTg                                                                                  <  1:442688/150‑1 (MQ=255)
                                                                 cTGTGATCGCGCGGCTGGTCGTCAGTGGACCTCCAGCGCTCAATAAATCCATGAACATTTGTTTTAACCACTGGGCTGCGAGGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTgg                                                                                 <  2:1050289/150‑1 (MQ=255)
                                                                 cTGGGATCGCGCGGCAGGTCGTCATGTTACCTCCAGAGCTCTCTAAATCACTTAACATTTGTTTTAACCTCGTGTCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTgg                                                                                 <  1:1516594/150‑1 (MQ=255)
                                                                  tGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTggg                                                                                <  2:223195/150‑1 (MQ=255)
                                                                   gggATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGc                                                                               <  1:1358975/150‑1 (MQ=255)
                                                                    ggATCCCGCGGCAGGTCGTCTTGTGACCTCCAGCTCTCACTAAATCAATTAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCt                                                                              <  2:240947/150‑1 (MQ=255)
                                                                     gATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTa                                                                             <  2:683223/150‑1 (MQ=255)
                                                                       gcgcgcgGCAGGTCGGCAGGTGACCTCCAGCGCTCACTAATTCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAg                                                                           <  1:111164/149‑1 (MQ=255)
                                                                         gcgcgGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTa                                                                         >  1:442699/1‑150 (MQ=255)
                                                                                aGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCTATGCCCCGTGGTTTTTTATTGTGTTGATGTGTTAGTAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCca                                                                  <  2:1695090/150‑1 (MQ=255)
                                                                                  gtcgtcAGGTGACCTCCAGCGCTCTCTAACTCACTGAAGATTTGTTTTAACCACGGGGCTGCGATGCCAGGTGGTTTTTTATTGTGTGGAAGGGTAAGGAAATGATGGAAACTAAGAACAAGCTAAACCGTATGCGAAGTACACGCcaca                                                                >  2:1065527/1‑150 (MQ=255)
                                                                                  gtcgtcAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGCGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCCCcaca                                                                >  1:501374/1‑150 (MQ=255)
                                                                                    cgtcAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACACGCCACATa                                                              >  2:900705/1‑150 (MQ=255)
                                                                                      tcAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATAcc                                                            >  1:36844/1‑149 (MQ=255)
                                                                                         ggTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCtt                                                         <  1:537340/150‑1 (MQ=255)
                                                                                         ggTGACATCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCtt                                                         >  2:27241/1‑150 (MQ=255)
                                                                                            gACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTAt                                                      >  2:493723/1‑150 (MQ=255)
                                                                                            gACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTAt                                                      <  1:1451379/150‑1 (MQ=255)
                                                                                              ccTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTAATTGTGTTGATGGGGTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATgg                                                    >  2:608125/1‑150 (MQ=255)
                                                                                               cTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGc                                                   >  1:730426/1‑150 (MQ=255)
                                                                                               cTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGc                                                   <  1:114816/150‑1 (MQ=255)
                                                                                                tCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCa                                                  >  1:559820/1‑150 (MQ=255)
                                                                                                    gcgcTCACTAAATCACTGCACATTTGTTTTAACCACGGCGCTGAGATCCCCCGTGGATTTTTATTGTGTTGATGGGTTAGGAATTGATGTAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGAACttg                                              >  1:699268/1‑148 (MQ=255)
                                                                                                    gcgcTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTgg                                              <  1:1661446/150‑1 (MQ=255)
                                                                                                    gcgcTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTgg                                              <  2:1100954/150‑1 (MQ=255)
                                                                                                      ncTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCGCCACATACGCTTTATGGCACTGGGtt                                           >  2:619292/2‑150 (MQ=255)
                                                                                                      gatCACTAAATCACTGAACTTTTTTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGATAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGt                                            <  2:690260/148‑1 (MQ=255)
                                                                                                        tCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATTCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTc                                          <  1:834738/150‑1 (MQ=255)
                                                                                                         cACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCa                                         >  2:1189264/1‑150 (MQ=255)
                                                                                                         cACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCa                                         >  1:970821/1‑150 (MQ=255)
                                                                                                         cACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCa                                         >  2:364245/1‑150 (MQ=255)
                                                                                                           cTAAATCACTGTACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGc                                       <  2:687128/150‑1 (MQ=255)
                                                                                                              aaTCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAAt                                    <  1:1093875/150‑1 (MQ=255)
                                                                                                               aTCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAAtt                                   <  2:912053/150‑1 (MQ=255)
                                                                                                                tCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTg                                  <  1:569851/150‑1 (MQ=255)
                                                                                                                  aCTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGc                                <  1:1209993/150‑1 (MQ=255)
                                                                                                                  aCTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGc                                >  1:1329410/1‑150 (MQ=255)
                                                                                                                     gAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAATAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAAATGGCTcc                             >  2:248846/1‑150 (MQ=255)
                                                                                                                      aaCATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAATAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCg                            <  2:1489007/150‑1 (MQ=255)
                                                                                                                         aTTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGc                         >  1:646745/1‑150 (MQ=255)
                                                                                                                          tttGTTTTAACCACGGGGCTGCGATGCCCCGTGGTATTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGTGTTCAGCAATTGGCTCGGGGCt                        >  1:331169/1‑150 (MQ=255)
                                                                                                                             gTTTTAACCACGGGGCTGCGATGCCCCGAGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAATCGTGGGCTAAGTACCAGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGtt                     >  2:233696/1‑150 (MQ=255)
                                                                                                                              ttttAACCACGGGGCTGCGATGCCCCGTGGTTTTTTAGTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATTCGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCttgtt                    >  1:1518335/1‑146 (MQ=255)
                                                                                                                               tttAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTCCCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGtttt                   >  1:830354/1‑150 (MQ=255)
                                                                                                                               tttAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGGTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATAGGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGgttt                   >  1:1621790/1‑150 (MQ=255)
                                                                                                                                ttAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTa                  <  1:493723/150‑1 (MQ=255)
                                                                                                                                   aCCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACgg               <  1:27241/150‑1 (MQ=255)
                                                                                                                                   aCCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACgg               >  2:143993/1‑150 (MQ=255)
                                                                                                                                   aCCACGGGGCATCGCAGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACgg               <  2:1198916/150‑1 (MQ=255)
                                                                                                                                    ccACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGt              >  1:1271932/1‑150 (MQ=255)
                                                                                                                                    ccACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGt              <  2:1271932/150‑1 (MQ=255)
                                                                                                                                    ccACGGGGCATCGCAGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGt              >  1:520296/1‑150 (MQ=255)
                                                                                                                                    ccACGGGGCATCGCAGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGt              >  2:1563586/1‑150 (MQ=255)
                                                                                                                                         gggCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGc         <  1:459036/150‑1 (MQ=255)
                                                                                                                                         gggCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCGGTTTTACGGTTCGGc         <  1:1285926/150‑1 (MQ=255)
                                                                                                                                          ggCTGCGATGCCCCGAGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCAGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCa        <  1:233696/150‑1 (MQ=255)
                                                                                                                                                gATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGcc  <  2:1474056/150‑1 (MQ=255)
                                                                                                                                                gATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTAATCGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGcc  <  1:248846/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCC  >  NC_000913/420800‑421093

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