Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,488,624 T→C R131R (CGT→CGC aldA → aldehyde dehydrogenase A, NAD‑linked

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,488,6240TC94.4% 273.1 / 6.2 89R131R (CGT→CGCaldAaldehyde dehydrogenase A, NAD‑linked
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base C (33/51);  total (36/53)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-02

GCCAGTGAAATCAGTGCGCTGATTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGTCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCTCCCGCTCTTTTGACCGGTAATACCATCGTCATTAAACCTAGT  >  NC_000913/1488478‑1488765
                                                                                                                                                  |                                                                                                                                             
gcCAGTGAAATCAGTGCGCTGATTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCa                                                                                                                                            <  1:943134/150‑1 (MQ=255)
 ccAGTGAAATCAGTGCGCTGATTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAgg                                                                                                                                          <  1:394406/151‑1 (MQ=255)
 ccAGTGAAATCAGTGCGCTGATTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAgg                                                                                                                                          <  2:846519/151‑1 (MQ=255)
 ccAGTGAAATCAGTGCGCTGATTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAgg                                                                                                                                          <  1:16630/151‑1 (MQ=255)
   aGTGAAATCAGTGCGCTGATTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGgag                                                                                                                                        <  2:963395/151‑1 (MQ=255)
   aGTGAAATCAGTGCGCTGATTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGgag                                                                                                                                        <  2:992274/151‑1 (MQ=255)
    gTGAAATCAGTGCGCTGATTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGgaga                                                                                                                                       <  2:658102/151‑1 (MQ=255)
      gAAATCAGTGCGCTGATTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGaaa                                                                                                                                     <  1:1127531/151‑1 (MQ=255)
       aaaTCAGTGCGCTGATTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAATCGATCGCCCAGGAGaaaa                                                                                                                                    <  1:823587/151‑1 (MQ=255)
       aaaTCAGTGCGCTGATTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGaaaa                                                                                                                                    <  1:297702/151‑1 (MQ=255)
        aagcAGTGCGCTGATTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAat                                                                                                                                   <  1:1144281/148‑1 (MQ=255)
         aTCAGTGCGCTGATTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAata                                                                                                                                  <  1:1030425/151‑1 (MQ=255)
         aTCAGTGCGCTGATTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAat                                                                                                                                   >  1:7051/1‑150 (MQ=255)
           cAGTGCGCTGATTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAatat                                                                                                                                 <  2:1189353/150‑1 (MQ=255)
              tGCGCTGATTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCtt                                                                                                                             >  1:590138/1‑151 (MQ=255)
               gcgcTGATTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCttt                                                                                                                            <  2:6774/151‑1 (MQ=255)
                cgcTGATTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCtttt                                                                                                                           <  2:884536/151‑1 (MQ=255)
                   tGATTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGtt                                                                                                                        >  1:814973/1‑151 (MQ=255)
                    gATTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGTCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGtt                                                                                                                        <  1:404314/150‑1 (MQ=255)
                    gATTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGttt                                                                                                                       <  1:1002028/151‑1 (MQ=255)
                     aTTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGTCCAGGAGAAAATATTCTTTTGttt                                                                                                                       >  1:400437/1‑150 (MQ=255)
                     aTTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTa                                                                                                                      <  1:654543/151‑1 (MQ=255)
                       tgttgAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTaa                                                                                                                     <  1:1367490/150‑1 (MQ=255)
                         ttgAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACg                                                                                                                  >  2:1135997/1‑151 (MQ=255)
                            aagaagGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTgc                                                                                                               <  2:291564/151‑1 (MQ=255)
                                  ggggCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTgg                                                                                                         >  2:583880/1‑151 (MQ=255)
                                   gggCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTgg                                                                                                         >  1:692663/1‑150 (MQ=255)
                                   gggCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGgt                                                                                                        <  2:357210/151‑1 (MQ=255)
                                    ggCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGgtg                                                                                                       <  1:1229504/151‑1 (MQ=255)
                                     gCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGgtgt                                                                                                      <  1:1175646/151‑1 (MQ=255)
                                       aaGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGtgtg                                                                                                     <  1:170999/150‑1 (MQ=255)
                                          aTCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGTCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTa                                                                                                 <  2:177393/151‑1 (MQ=255)
                                            ccAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTAc                                                                                                <  2:1412622/150‑1 (MQ=255)
                                                  gctggctgAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCAtt                                                                                         >  1:582361/1‑151 (MQ=255)
                                                    tggctgAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCt                                                                                       >  2:1269730/1‑151 (MQ=255)
                                                     ggctgAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTg                                                                                      >  1:475736/1‑151 (MQ=255)
                                                        tgAAGTCGAAGCGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTTCGAGGGCGAGATTATTCAAAGCGATCGCCCCGGCGAAAATATTCTTTTGTTTAAACGTGCGCTTGGAGTGACTACCGGCATTATGCCg                                                                                   >  1:639191/1‑151 (MQ=255)
                                                          aaGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGt                                                                                  >  2:1113350/1‑150 (MQ=255)
                                                          aaGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTg                                                                                 <  1:1038750/151‑1 (MQ=255)
                                                            gTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGa                                                                               <  2:446842/151‑1 (MQ=255)
                                                             tCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGa                                                                               >  1:315132/1‑150 (MQ=255)
                                                              cGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGCAc                                                                             >  2:1366608/1‑151 (MQ=255)
                                                              cGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAAc                                                                             >  1:1355926/1‑151 (MQ=255)
                                                               gAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACt                                                                            >  2:223883/1‑151 (MQ=255)
                                                               gAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACt                                                                            <  1:474848/151‑1 (MQ=255)
                                                               gAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACt                                                                            <  2:796990/151‑1 (MQ=255)
                                                                aaGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACtt                                                                           >  1:1095442/1‑151 (MQ=255)
                                                                    ggCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCg                                                                       >  2:1205612/1‑151 (MQ=255)
                                                                    ggCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTAAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATACTGCCGTGGAACTTCCCg                                                                       <  1:1033132/151‑1 (MQ=255)
                                                                      cTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGTCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGGGTGACTACCGGCATTCTGCCGTGGAACTTCCCGtt                                                                     >  1:293479/1‑151 (MQ=255)
                                                                              gCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCt                                                               <  1:1232947/149‑1 (MQ=255)
                                                                              gCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCa                                                             <  2:300162/151‑1 (MQ=255)
                                                                                cGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTACTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCAtt                                                           <  2:1064667/151‑1 (MQ=255)
                                                                                  aCTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCACAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTg                                                          >  1:1018315/1‑150 (MQ=255)
                                                                                   cTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGcc                                                        >  1:267049/1‑151 (MQ=255)
                                                                                   cTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGAGGAACTTCCCGTTCTTCCTCATTGcc                                                        >  1:167373/1‑151 (MQ=255)
                                                                                    tatcacgtttACCTGGCGGGGGGGGCTCGGCCTCACGCGGGCGAGATGATATGAAGCGCTCGCCCAGGAGACAATATTCTTTTGGTTAACCGACCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCAGTTCTGCCTTAtggacc                                                       >  1:273560/9‑146 (MQ=255)
                                                                                      tatCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGc                                                     >  2:519570/1‑151 (MQ=255)
                                                                                       atCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGc                                                     <  2:907465/150‑1 (MQ=255)
                                                                                         cGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGTCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCaa                                                   >  1:573924/1‑150 (MQ=255)
                                                                                          gATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCaaaa                                                 <  2:73103/151‑1 (MQ=255)
                                                                                          gATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCaaa                                                  >  2:51681/1‑150 (MQ=255)
                                                                                          gATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCaaa                                                  >  2:1317656/1‑150 (MQ=255)
                                                                                          gATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCaaa                                                  >  2:296210/1‑150 (MQ=255)
                                                                                           aTTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAAt                                                <  2:441487/151‑1 (MQ=255)
                                                                                             tACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATg                                               <  1:33755/150‑1 (MQ=255)
                                                                                              aCATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGGGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATg                                               >  1:697084/1‑149 (MQ=255)
                                                                                               cATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCt                                            <  2:389680/151‑1 (MQ=255)
                                                                                                aTGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCTc                                           >  2:193814/1‑151 (MQ=255)
                                                                                                         tGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCTCCCGCTCttt                                  >  1:647340/1‑151 (MQ=255)
                                                                                                          gggCTCGGCGTTACGAGGGCGGGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTTCGCTTGGTGTGACTACCGGCATTTTGCCGTGGAACTTTCCGTTCTTCCTCATTGACCGCAAAATGGCTCCCGCTCtctt                                 <  1:1228062/151‑1 (MQ=255)
                                                                                                          gggCACGGCTTTACGAGGGCGAGATTATTCAAAGCAATCGCGCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCTCCCGCTCttg                                  <  1:529712/150‑2 (MQ=255)
                                                                                                           ggCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCTCCCGCTCTTTTg                                <  2:651005/151‑1 (MQ=255)
                                                                                                              aCGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCTCCCGCTCTTTTGAcc                             <  1:1331457/151‑1 (MQ=255)
                                                                                                                 gCGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCTCCCGCTCTTTTGACCgg                           <  2:1251/150‑1 (MQ=255)
                                                                                                                  cGTTACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCTCCCGCTCTTTTGACCGGt                          >  2:299797/1‑150 (MQ=255)
                                                                                                                     tACGAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCTCCCGCTCTTTTGACCGGTAAt                       >  1:358468/1‑150 (MQ=255)
                                                                                                                        gAGGGCGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCTCCCGCTCTTTTGACCGGTAATGCCa                   <  1:879493/151‑1 (MQ=255)
                                                                                                                             cGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCTCCCGCTCTTTTGACCGGTAATACCATCGTc              <  1:247022/151‑1 (MQ=255)
                                                                                                                             cGAGATTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCTCCCGCTCTTTTGACCGGTAATACCATCGTc              <  2:1123342/151‑1 (MQ=255)
                                                                                                                               agaTTATTCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCTCCCGCTCTTTTGACCGGTAATACCATCGTCAt            >  1:550489/1‑151 (MQ=255)
                                                                                                                                gaTTATTCAAAGCGATCGTCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCTCCCGCTCTTTTGACCGGTAATACCATCGTCAtt           <  2:400451/151‑1 (MQ=255)
                                                                                                                                  ttattCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCTCCCGCTCTTTTGACCGGTAATACCATCGTCATTa          >  1:1374019/1‑150 (MQ=255)
                                                                                                                                   tattCAAAGCGATCGCCCAGGAGAAAATAGTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCTCCCGCTCTTTTGACCGGTAATACCATCGTCATTaa         <  2:815010/150‑1 (MQ=255)
                                                                                                                                    attCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCTCCCGCTCTTTTGACCGGTAATACCATCGTCATTaaa        <  2:475758/150‑1 (MQ=255)
                                                                                                                                    attCAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCTCCCGCTCTTTTGACCGGTAATACCATCGTCATTaaa        <  2:272889/150‑1 (MQ=255)
                                                                                                                                       cAAAGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCTCCCGCTCTTTTGACCGGTAATACCATCGTCATTAAACCTa    <  1:646836/151‑1 (MQ=255)
                                                                                                                                        acaGCGATCGCCCAGGAGAAATTATTCTTTTGTTTAACCGTCCGCTTGGTGTGACTGCCGGCATTCTGCCGTGGAGCTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCTCCCGCTCTTTTGGCCGGTAATACCATCGTCATTAAACCTa    <  1:287468/148‑1 (MQ=255)
                                                                                                                                         aaGCGATCGCCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCTCCCGCTCTTTTGACCGGTAATACCATCGTCATTAAACCTAGt  >  1:518289/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
GCCAGTGAAATCAGTGCGCTGATTGTTGAAGAAGGGGGCAAGATCCAGCAGCTGGCTGAAGTCGAAGTGGCTTTTACTGCCGACTATATCGATTACATGGCGGAGTGGGCACGGCGTTACGAGGGCGAGATTATTCAAAGCGATCGTCCAGGAGAAAATATTCTTTTGTTTAAACGTGCGCTTGGTGTGACTACCGGCATTCTGCCGTGGAACTTCCCGTTCTTCCTCATTGCCCGCAAAATGGCTCCCGCTCTTTTGACCGGTAATACCATCGTCATTAAACCTAGT  >  NC_000913/1488478‑1488765

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