Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,043,629 A→G I185T (ATT→ACT)  yihG ← putative endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,043,6290AG99.4% 1121.3 / ‑3.8 355I185T (ATT→ACT) yihGputative endonuclease
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (149/204);  total (150/205)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 390 total aligned reads displayed.

GCCGCCTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAATGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTAATAAATAA  >  NC_000913/4043539‑4043718
                                                                                          |                                                                                         
gccgccTTTGGTGGCAATAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtg                                                                                         <  2:1498681/93‑1 (MQ=255)
gccgccTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtg                                                                                         <  2:7899250/93‑1 (MQ=255)
  cgccTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtggt                                                                                       <  2:3217233/93‑1 (MQ=255)
  cgccTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtggt                                                                                       <  1:3612471/93‑1 (MQ=255)
    ccTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtggtgg                                                                                     <  2:3938442/93‑1 (MQ=255)
    ccTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtggtgg                                                                                     >  1:4327512/1‑93 (MQ=255)
      tttGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAgtggtggtg                                                                                   >  1:5258702/1‑93 (MQ=255)
       ttGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtggtggtgg                                                                                  <  1:2623017/93‑1 (MQ=255)
        tggtggCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTggg                                                                                 <  1:5152805/93‑1 (MQ=255)
           tggCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATg                                                                              >  1:4636594/1‑93 (MQ=255)
           tggCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATg                                                                              >  1:3109960/1‑93 (MQ=255)
           tggCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATg                                                                              >  2:6617443/1‑93 (MQ=255)
              caacaaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAg                                                                           >  1:5081139/1‑93 (MQ=255)
               aacaaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAgg                                                                          >  2:7785567/1‑93 (MQ=255)
               aacaaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAgg                                                                          >  1:5406795/1‑93 (MQ=255)
               aacaaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAgg                                                                          >  2:2297016/1‑93 (MQ=255)
                acaaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGc                                                                         <  1:881435/93‑1 (MQ=255)
                acaaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGc                                                                         <  2:6162568/93‑1 (MQ=255)
                 caaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCg                                                                        >  2:6762535/1‑93 (MQ=255)
                 caaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCg                                                                        >  2:5449828/1‑93 (MQ=255)
                  aaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGa                                                                       <  1:5478026/93‑1 (MQ=255)
                  aaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGa                                                                       <  1:1956319/93‑1 (MQ=255)
                   aGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGaa                                                                      <  1:718293/93‑1 (MQ=255)
                     ttttGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGaaat                                                                    >  1:6579169/1‑92 (MQ=255)
                     ttttGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAAc                                                                    <  1:5819671/93‑1 (MQ=255)
                        tGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACttt                                                                 <  1:6022961/93‑1 (MQ=255)
                        tGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACttt                                                                 <  2:1680446/93‑1 (MQ=255)
                         gAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACtttt                                                                <  1:1291817/93‑1 (MQ=255)
                           aaaGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTca                                                              <  2:2962418/93‑1 (MQ=255)
                           aaaGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTca                                                              <  1:682809/93‑1 (MQ=255)
                             agtggagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTcaca                                                            <  2:3496073/93‑1 (MQ=255)
                              gtggagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACaa                                                           <  2:7191329/93‑1 (MQ=255)
                              gtggagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACaa                                                           >  1:1014399/1‑93 (MQ=255)
                                ggagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGa                                                         >  1:6291570/1‑93 (MQ=255)
                                ggagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGa                                                         <  1:3739292/93‑1 (MQ=255)
                                 gagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGaa                                                        <  2:6901643/93‑1 (MQ=255)
                                 gagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGaa                                                        <  1:6865254/93‑1 (MQ=255)
                                 gagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGaa                                                        <  1:6953004/93‑1 (MQ=255)
                                 gagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACAAAATTGACAGTGGTGGTAGGATGCAGGCGAAACTTTTCACAAGaa                                                        <  2:2067373/93‑1 (MQ=255)
                                  agtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAAc                                                       <  1:1881282/93‑1 (MQ=255)
                                   gtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACg                                                      <  1:1611597/93‑1 (MQ=255)
                                   gtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACg                                                      >  2:533722/1‑93 (MQ=255)
                                    tggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGc                                                     <  1:7681657/93‑1 (MQ=255)
                                    tggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGc                                                     <  1:4717130/93‑1 (MQ=255)
                                        tttGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAg                                                 <  1:7886313/93‑1 (MQ=255)
                                         ttGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGt                                                <  2:4146038/93‑1 (MQ=255)
                                         ttGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGt                                                <  2:5068918/93‑1 (MQ=255)
                                          tGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTg                                               <  1:7507787/93‑1 (MQ=255)
                                          tGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTg                                               >  2:7660502/1‑93 (MQ=255)
                                           gCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTgg                                              >  1:7312976/1‑93 (MQ=255)
                                            gttgtGTTTTTTCTGCGGGGAACCGGAACCCTTAAAGAAATTGAACAGGGGGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGt                                             <  2:7711584/89‑1 (MQ=255)
                                            cTGAGGTTTTTCTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGtt                                            <  2:284444/93‑1 (MQ=255)
                                             tGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGtt                                            <  1:211563/93‑1 (MQ=255)
                                             tGATGTTTTTCTTGCGTGAAACGGCAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGtt                                            <  1:2900610/93‑1 (MQ=255)
                                               aTGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTTGTTTc                                          >  2:6994168/1‑93 (MQ=255)
                                               aTGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGGGGTTTc                                          >  2:4916022/1‑93 (MQ=255)
                                                 gTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCaa                                        >  2:2648557/1‑93 (MQ=255)
                                                   ttttCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACa                                      <  2:4210330/93‑1 (MQ=255)
                                                    tttCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACAt                                     <  1:1793238/93‑1 (MQ=255)
                                                     ttCTTTCCTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATc                                    <  1:4271630/93‑1 (MQ=255)
                                                       cTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCtt                                  <  1:3861468/93‑1 (MQ=255)
                                                       cTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCtt                                  <  1:4989979/93‑1 (MQ=255)
                                                        ttGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCttt                                 >  1:3417172/1‑93 (MQ=255)
                                                        ttGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCttt                                 >  2:3293544/1‑93 (MQ=255)
                                                         tGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTa                                <  2:7418527/93‑1 (MQ=255)
                                                         tGCGCGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTa                                <  2:2896104/93‑1 (MQ=255)
                                                          gCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTAc                               >  2:5722813/1‑93 (MQ=255)
                                                          gCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTAc                               <  1:2374838/93‑1 (MQ=255)
                                                           cGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTAcc                              <  1:321298/93‑1 (MQ=255)
                                                            gTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCa                             <  2:803836/93‑1 (MQ=255)
                                                            gTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCa                             <  1:5730526/93‑1 (MQ=255)
                                                              gAAACGGGAGCCTTCAACGAAATTAACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAGCGCCGAGCGGTTTCAAAATCTTTACCACg                           >  2:7586054/1‑93 (MQ=255)
                                                                aaCGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTc                         >  2:4613920/1‑93 (MQ=255)
                                                                 aCGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCg                        >  1:6701336/1‑93 (MQ=255)
                                                                    ggAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTc                     <  2:531766/93‑1 (MQ=255)
                                                                     gAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTcc                    <  2:3328757/93‑1 (MQ=255)
                                                                     gAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTcc                    <  2:2390344/93‑1 (MQ=255)
                                                                      aGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCg                   <  1:2754557/93‑1 (MQ=255)
                                                                       gCCTTCAACGAAATTGACAGTGGTGGTGGGACGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCgg                  >  1:611795/1‑93 (MQ=255)
                                                                        ccTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGa                 >  2:6278017/1‑93 (MQ=255)
                                                                          ttCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGa                 <  2:3883170/91‑1 (MQ=255)
                                                                          ttCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATg               >  1:4896795/1‑93 (MQ=255)
                                                                          ttCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATg               >  2:1922650/1‑93 (MQ=255)
                                                                          ttCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATg               <  2:6557746/93‑1 (MQ=255)
                                                                          ttCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATg               >  2:1132335/1‑93 (MQ=255)
                                                                          ttCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATg               >  1:6463000/1‑93 (MQ=255)
                                                                           tCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGa              >  1:2077943/1‑93 (MQ=255)
                                                                           tCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGa              <  1:458842/93‑1 (MQ=255)
                                                                             aaCGAAATTGACAGTGGTGGGGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACg            >  1:2536331/1‑93 (MQ=255)
                                                                               cGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGta          <  2:431602/93‑1 (MQ=255)
                                                                                gAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGtaa         <  1:4005569/93‑1 (MQ=255)
                                                                                gAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGtaa         <  2:3120614/93‑1 (MQ=255)
                                                                                gAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGtaa         >  2:4081856/1‑93 (MQ=255)
                                                                                 aaaTTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGtaat        >  2:1867904/1‑93 (MQ=255)
                                                                                 aaaTTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGtaat        <  2:2002906/93‑1 (MQ=255)
                                                                                  aagtGACAGTGGGGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGtaata       <  2:5598247/90‑1 (MQ=255)
                                                                                     tGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTaataaat    <  2:3712737/93‑1 (MQ=255)
                                                                                     tGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTaataaat    <  2:4179670/93‑1 (MQ=255)
                                                                                      gACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTaataaata   <  2:3596724/93‑1 (MQ=255)
                                                                                       aCAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTAataaataa  >  2:1366419/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCCGCCTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAATGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTAATAAATAA  >  NC_000913/4043539‑4043718

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