Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,744,085 G→T 5.0% E115* (GAA→TAA)  yiaL → DUF386 domain‑containing protein YiaL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,744,0850GT 5.0% 340.2 / 4.5 99E115* (GAA→TAA) yiaLDUF386 domain‑containing protein YiaL
Reads supporting (aligned to +/- strand):  ref base G (43/51);  new base T (3/2);  total (46/53)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 103 total aligned reads displayed.

ATCGTCGGTATATCGATATTCAGTTTCTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGCCAGAAC  >  NC_000913/3743936‑3744234
                                                                                                                                                     |                                                                                                                                                     
aTCGTCGGTATATCGATATTCAGTTTCTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAg                                                                                                                                                       <  1:290098/150‑1 (MQ=255)
 tcgtcgGTATATCGATATTCAGTTTCTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGa                                                                                                                                                      <  1:1027047/150‑1 (MQ=255)
    tcgGTATATCGATATTCAGTTTCTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAAt                                                                                                                                                   <  1:725007/150‑1 (MQ=255)
      ggTATATCGATATTCAGTTTCTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGa                                                                                                                                                 <  1:1223303/150‑1 (MQ=255)
       gTATATCGATATTCAGTTTCTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGAt                                                                                                                                                <  2:196285/150‑1 (MQ=255)
        tataTCGATATTCAGTTTCTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATa                                                                                                                                               <  2:769022/150‑1 (MQ=255)
          tatCGATATTCAGTTTCTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTCTATAAATGAttca                                                                                                                                             <  1:498605/150‑4 (MQ=255)
          tatCGATATTCAGTTTCTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATAcc                                                                                                                                             <  1:9312/150‑1 (MQ=255)
          tatCGATATTCAGTTTCTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATAcc                                                                                                                                             >  1:482751/1‑150 (MQ=255)
            tCGATATTCAGTTTCTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCgg                                                                                                                                           >  2:429507/1‑150 (MQ=255)
            tCGATATTCAGTTTCTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCgg                                                                                                                                           >  1:1110989/1‑150 (MQ=255)
              gATATTCAGTTTCTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGc                                                                                                                                         <  1:610797/150‑1 (MQ=255)
                 atTCAGTTTCTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGc                                                                                                                                      <  1:1228924/150‑1 (MQ=255)
                  ttCAGTTTCTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCt                                                                                                                                     >  1:19287/1‑150 (MQ=255)
                   tCAGTTTCTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTa                                                                                                                                    >  1:1431092/1‑150 (MQ=255)
                      gTTTCTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGTAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGc                                                                                                                                 <  1:1034015/150‑1 (MQ=255)
                      gTTTCTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGc                                                                                                                                 <  2:793578/150‑1 (MQ=255)
                         tCTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAAt                                                                                                                              <  2:659012/150‑1 (MQ=255)
                         tCTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAAt                                                                                                                              >  1:557012/1‑150 (MQ=255)
                          cTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAAtt                                                                                                                             <  1:764620/150‑1 (MQ=255)
                             gCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATtttc                                                                                                                          <  1:726972/150‑1 (MQ=255)
                              cGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACAGGGCAGCTACGAAATtttct                                                                                                                         >  1:1296480/1‑150 (MQ=255)
                               gTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATtttctt                                                                                                                        <  1:1099402/150‑1 (MQ=255)
                                tGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATtttcttt                                                                                                                       <  1:781446/150‑1 (MQ=255)
                                     cGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCa                                                                                                                  <  2:1334076/150‑1 (MQ=255)
                                     cGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCa                                                                                                                  <  2:359615/150‑1 (MQ=255)
                                      gAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAg                                                                                                                 <  2:379237/150‑1 (MQ=255)
                                         gagaAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGAt                                                                                                              >  1:1286496/1‑150 (MQ=255)
                                         gagaAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGAt                                                                                                              <  1:526978/150‑1 (MQ=255)
                                           taaaaTCGGCATTGCTATTGATACTGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGt                                                                                                            <  1:1261881/149‑1 (MQ=255)
                                                tCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCAcc                                                                                                       <  1:344569/150‑1 (MQ=255)
                                                    cgttGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATAcc                                                                                                                                             >  2:867131/3‑108 (MQ=255)
                                                    cgttGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATAcc                                                                                                                                             <  1:867131/106‑1 (MQ=255)
                                                     aTTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATTTATTATTTATATAAATGATACCTGGGAGATACTCAATTTTCTTTACGCATTATGTTCAcctagga                                                                                                  <  2:1167850/150‑6 (MQ=255)
                                                       tGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAgg                                                                                                <  1:1110456/150‑1 (MQ=255)
                                                        gCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGt                                                                                               >  2:1320157/1‑150 (MQ=255)
                                                          tATTGATACGGGAAATAATAAAGGCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTg                                                                                             <  1:1095774/150‑1 (MQ=255)
                                                           aTTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGt                                                                                            <  1:172262/150‑1 (MQ=255)
                                                           aTTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGt                                                                                            <  2:482751/150‑1 (MQ=255)
                                                              gATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGtatt                                                                                         <  1:483577/150‑1 (MQ=255)
                                                              gATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGtatt                                                                                         >  1:1261336/1‑150 (MQ=255)
                                                              gATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGAAGATACGCAATTTTCTTTCCGCAGGATGTTCACATACCAGGTTGtatt                                                                                         >  1:1443604/1‑150 (MQ=255)
                                                               aTACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGtatta                                                                                        >  2:57688/1‑150 (MQ=255)
                                                                tACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGtattat                                                                                       <  1:1477887/150‑1 (MQ=255)
                                                                 aCGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATg                                                                                      >  1:234264/1‑150 (MQ=255)
                                                                   gggAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCa                                                                                    >  2:1267768/1‑150 (MQ=255)
                                                                   gggAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCa                                                                                    <  1:1059722/150‑1 (MQ=255)
                                                                    ggAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCaa                                                                                   >  2:1017235/1‑150 (MQ=255)
                                                                        ataataAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAg                                                                               >  2:1063560/1‑150 (MQ=255)
                                                                        ataataAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAg                                                                               >  2:301395/1‑150 (MQ=255)
                                                                         taataaAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACa                                                                                >  2:714685/1‑148 (MQ=255)
                                                                         taataaAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACa                                                                                <  1:714685/148‑1 (MQ=255)
                                                                         taataaAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGc                                                                              <  2:546007/150‑1 (MQ=255)
                                                                             aaaGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCctct                                                                          >  2:48209/1‑150 (MQ=255)
                                                                              aaGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTg                                                                         <  1:1320157/150‑1 (MQ=255)
                                                                               aGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTga                                                                        <  1:70033/150‑1 (MQ=255)
                                                                                gTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTgag                                                                       >  2:214238/1‑150 (MQ=255)
                                                                                  cAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCctc                                                                           >  1:484351/1‑144 (MQ=255)
                                                                                  cAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCctc                                                                           <  2:484351/144‑1 (MQ=255)
                                                                                  cAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGAt                                                                     >  2:1073003/1‑150 (MQ=255)
                                                                                    gCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATAc                                                                   >  2:446212/1‑150 (MQ=255)
                                                                                     cGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACg                                                                  <  1:1073003/150‑1 (MQ=255)
                                                                                          cACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGaaaaa                                                             >  1:968098/1‑150 (MQ=255)
                                                                                              aTTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATtgt                                                         >  2:1150599/1‑150 (MQ=255)
                                                                                              aTTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATtgt                                                         >  1:95898/1‑150 (MQ=255)
                                                                                               ttAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATtgtg                                                        <  1:1118735/150‑1 (MQ=255)
                                                                                                  gagCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCATACAGCCTCTGAGATACTAAAAATTGTGGtt                                                     <  2:273766/150‑1 (MQ=255)
                                                                                                  gagCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGtt                                                     >  1:422406/1‑150 (MQ=255)
                                                                                                    gcagcGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTaa                                                   >  1:167715/1‑150 (MQ=255)
                                                                                                    gcagcGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTACGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTaa                                                   >  2:28088/1‑150 (MQ=255)
                                                                                                        agcAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACTAAAAATTGTGGTTAAAGTc                                               <  2:1296480/149‑1 (MQ=255)
                                                                                                         gcAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCg                                              >  2:228294/1‑150 (MQ=255)
                                                                                                         gcAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTAACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTCGTTAAAGTCg                                              >  1:1004401/1‑150 (MQ=255)
                                                                                                                 attTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACg                                      >  1:1401053/1‑150 (MQ=255)
                                                                                                                     tttATCACGACATTGAAAATGAATCATTTATAGAAATGATACAGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGATATACGAAAAATTGTGGTTAAAGTCGCATTAACGgcgc                                  <  2:1443604/150‑1 (MQ=255)
                                                                                                                       tATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTg                                <  2:1110989/150‑1 (MQ=255)
                                                                                                                       tATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTg                                <  2:578642/150‑1 (MQ=255)
                                                                                                                       tATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTg                                >  1:748882/1‑150 (MQ=255)
                                                                                                                         tCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGaa                              >  1:881987/1‑150 (MQ=255)
                                                                                                                         tCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGaa                              >  1:1146734/1‑150 (MQ=255)
                                                                                                                          cACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAAt                             >  1:509228/1‑150 (MQ=255)
                                                                                                                          cACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGATGTATTATCCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCCCTGAAt                             >  1:504123/1‑150 (MQ=255)
                                                                                                                           aCGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAAtt                            >  1:1170300/1‑150 (MQ=255)
                                                                                                                           aCGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGCTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGCCGCTGAAtt                            >  1:826365/1‑150 (MQ=255)
                                                                                                                            cGACAGTGAACATGAATCATTTATAGAAATGATACCGTTCAGCTACGCAATTTTGTTTCAGCAGGATGTTCACCGACCAGGTTTTATTATGCAAACAGCCTATGAGATACGAAAAATTGTGGTTAAAGTCGGTTTATCGGCGCTGAATTa                           <  2:1431092/150‑1 (MQ=255)
                                                                                                                                  tGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGaa                     <  2:151310/150‑1 (MQ=255)
                                                                                                                                     cccggtaTCATTTCTATAAATGATACCGGGCAGCTACGCAATTTTTTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGaataa                  >  1:941721/7‑150 (MQ=255)
                                                                                                                                     aCATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGaataa                  <  1:1452233/150‑1 (MQ=255)
                                                                                                                                      ccggtaTCATTTCTATAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGaataat                 >  2:1341482/6‑150 (MQ=255)
                                                                                                                                       cggtaTCATTTCTATAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGACGCATTAACGGCGCTGAATTAATCGAATAATc                >  2:46597/5‑150 (MQ=255)
                                                                                                                                        tGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCg               <  1:214238/150‑1 (MQ=255)
                                                                                                                                          aaTCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTc             <  1:301395/150‑1 (MQ=255)
                                                                                                                                           aTCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCt            <  2:557012/150‑1 (MQ=255)
                                                                                                                                                ttATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGcc       <  2:1401053/150‑1 (MQ=255)
                                                                                                                                                 tataGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGccc      >  2:1385210/1‑149 (MQ=255)
                                                                                                                                                  ataGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGCCAg     >  2:583515/1‑150 (MQ=255)
                                                                                                                                                   taGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGCCAGa    <  1:755595/150‑1 (MQ=255)
                                                                                                                                                     gAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGCCAGAAc  >  2:463925/1‑150 (MQ=255)
                                                                                                                                                     gAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGCCAGAAc  >  1:1307804/1‑150 (MQ=255)
                                                                                                                                                     gAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGCCAGAAc  >  2:222231/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCGTCGGTATATCGATATTCAGTTTCTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTACGCAATTTTCTTTCCGCAGGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGCCAGAAC  >  NC_000913/3743936‑3744234

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