Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,410,957 T→G 5.6% intergenic (+73/+56) dbpA → / ← ttcA ATP‑dependent RNA helicase DbpA/tRNA cytosine(32) 2‑sulfurtransferase TtcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,410,9570TG 5.6% 521.1 / 16.7 161intergenic (+73/+56)dbpA/ttcAATP‑dependent RNA helicase DbpA/tRNA cytosine(32) 2‑sulfurtransferase TtcA
Reads supporting (aligned to +/- strand):  ref base T (85/67);  new base G (5/4);  total (90/71)
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.53e-01

Only 100 of 163 total aligned reads displayed.

TTGCTCATAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTT  >  NC_000913/1410808‑1411104
                                                                                                                                                     |                                                                                                                                                   
ttGCTCATAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCt                                                                                                                                                     <  2:932099/150‑1 (MQ=255)
 tGCTCATAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCtt                                                                                                                                                    >  1:1356496/1‑150 (MQ=255)
   cTCATAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTAt                                                                                                                                                  <  1:217760/150‑1 (MQ=255)
   cTCATAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTAt                                                                                                                                                  >  2:1428985/1‑150 (MQ=255)
     cATAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATcc                                                                                                                                                >  2:773323/1‑150 (MQ=255)
       tAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCgg                                                                                                                                              <  2:793484/150‑1 (MQ=255)
        aaaGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGc                                                                                                                                             >  2:1074566/1‑150 (MQ=255)
        aaaGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGACGGGTGAAAGAGTAAACATCTTATTCGGGATTGCCGGATGAGACGCTGGCAGCGTCTTATCAGGc                                                                                                                                             >  1:574135/1‑150 (MQ=255)
         aaGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGcc                                                                                                                                            >  1:372841/1‑150 (MQ=255)
         aaGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGcc                                                                                                                                            <  1:1414490/150‑1 (MQ=255)
         aaGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGcc                                                                                                                                            <  1:1395650/150‑1 (MQ=255)
         aaGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTAATCGCGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGcc                                                                                                                                            >  2:1306192/1‑150 (MQ=255)
         aaGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGGGCAAGAGTAAACATCTTATTCGGGATTGCCGGATAAGACGCGGCCAGCGTCGCATCCGGca                                                                                                                                            <  1:1363077/150‑2 (MQ=255)
           gCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGGCGCTGGCCGCGTCTTATCCGGCCTc                                                                                                                                          >  1:866833/1‑150 (MQ=255)
            cATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTcc                                                                                                                                         <  1:1187839/150‑1 (MQ=255)
            cATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGAGACGCTGGCCGCGTCTTATCCTGCCTca                                                                                                                                         >  1:1003231/1‑149 (MQ=255)
             aTGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCa                                                                                                                                        >  1:1240334/1‑150 (MQ=255)
             aTGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCa                                                                                                                                        >  2:1051117/1‑150 (MQ=255)
               ggAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATa                                                                                                                                      <  1:1428985/150‑1 (MQ=255)
               ggAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGACTCCATa                                                                                                                                      >  2:246084/1‑150 (MQ=255)
                 aaaCAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAag                                                                                                                                    <  2:1133313/150‑1 (MQ=255)
                  aaCAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAaga                                                                                                                                   <  1:1354012/150‑1 (MQ=255)
                      gTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTac                                                                                                                               >  1:1208931/1‑149 (MQ=255)
                        tACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGcc                                                                                                                             >  1:961021/1‑150 (MQ=255)
                           aGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGa                                                                                                                          >  2:6621/1‑150 (MQ=255)
                                 ggAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCt                                                                                                                    >  2:1361511/1‑150 (MQ=255)
                                 ggAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCt                                                                                                                    >  2:231353/1‑150 (MQ=255)
                                  gAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCtt                                                                                                                   >  1:745488/1‑150 (MQ=255)
                                      aTTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGATTgcgc                                                                                                               >  2:363426/1‑150 (MQ=255)
                                      aTTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGGAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGCCTCGCCCGATACGCTTgcgc                                                                                                               >  2:1365335/1‑150 (MQ=255)
                                          aaGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATGCGCAAGCGTATCg                                                                                                           <  1:1132168/150‑1 (MQ=255)
                                          aaGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCg                                                                                                           >  1:1274317/1‑150 (MQ=255)
                                            ggAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCggg                                                                                                         >  1:880790/1‑150 (MQ=255)
                                            ggAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCggg                                                                                                         >  2:733172/1‑150 (MQ=255)
                                            ggAAAAACCTGCCGGGTGCGGTTATTAAAATAATGAATTGTTGAATTGCCGTGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGTCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATGCGCAAGCGTATCggg                                                                                                         <  2:711328/150‑1 (MQ=255)
                                             gAAAAAGGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCTACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGgc                                                                                                        <  2:681984/150‑1 (MQ=255)
                                             gAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGgc                                                                                                        <  1:965476/150‑1 (MQ=255)
                                             gAAAAACGTGCCGGGAGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTACACATCTTATTCGCGATTGCCGGATGAGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCAGATACGCTTCAGCATCGGgc                                                                                                        >  1:686756/1‑150 (MQ=255)
                                               aaaaCGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGgcgc                                                                                                      <  2:7634/150‑1 (MQ=255)
                                                 aaagTGCCGGGTGCTGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTa                                                                                                    <  2:574135/147‑1 (MQ=255)
                                                 aaCGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCt                                                                                                                                           >  1:290694/1‑111 (MQ=255)
                                                 aaCGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCt                                                                                                                                           <  2:290694/111‑1 (MQ=255)
                                                      gCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTg                                                                                               >  1:834846/1‑150 (MQ=255)
                                                       ccGGGTGTGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTgg                                                                                              <  2:774939/150‑1 (MQ=255)
                                                        cGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGt                                                                                             <  2:618479/150‑1 (MQ=255)
                                                             gCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTAttt                                                                                        >  2:701081/1‑150 (MQ=255)
                                                              cGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTc                                                                                       >  2:540440/1‑150 (MQ=255)
                                                                  tattaAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTAt                                                                                          <  2:1006689/143‑1 (MQ=255)
                                                                    ttaAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCa                                                                                 >  2:467883/1‑150 (MQ=255)
                                                                     taAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCaa                                                                                >  2:886982/1‑150 (MQ=255)
                                                                       aaaTAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAAcc                                                                              <  1:608224/150‑1 (MQ=255)
                                                                        aataatGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACca                                                                             >  2:403437/1‑150 (MQ=255)
                                                                              gAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCa                                                                       >  2:1381249/1‑150 (MQ=255)
                                                                                 aTGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGcc                                                                    >  1:1387889/1‑150 (MQ=255)
                                                                                  tGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCg                                                                   >  2:1186814/1‑150 (MQ=255)
                                                                                  tGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCCATACGCTTGCGCATCCGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCg                                                                   >  1:370596/1‑150 (MQ=255)
                                                                                    ttGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTa                                                                 >  1:901601/1‑150 (MQ=255)
                                                                                     tGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTaa                                                                >  1:1322325/1‑150 (MQ=255)
                                                                                      gAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTgg                                                                                              >  2:797526/1‑119 (MQ=255)
                                                                                      gAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTgg                                                                                              <  1:797526/119‑1 (MQ=255)
                                                                                      gAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAAc                                                               <  1:242916/150‑1 (MQ=255)
                                                                                       aaTTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACt                                                              <  1:885027/150‑1 (MQ=255)
                                                                                        aTTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTc                                                             <  1:773323/150‑1 (MQ=255)
                                                                                           gCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATc                                                          >  1:1000042/1‑150 (MQ=255)
                                                                                            ccGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATcc                                                         <  2:372841/150‑1 (MQ=255)
                                                                                             cGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATcc                                                                                                                                                >  2:700826/1‑62 (MQ=255)
                                                                                             cGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATcc                                                                                                                                                <  1:700826/62‑1 (MQ=255)
                                                                                                   cAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGAt                                                  >  1:24794/1‑150 (MQ=255)
                                                                                                     agagTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGAttt                                                <  1:678105/150‑1 (MQ=255)
                                                                                                          aaaCATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCAtc                                           >  2:967046/1‑150 (MQ=255)
                                                                                                           aaCATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCAtct                                          <  2:827093/150‑1 (MQ=255)
                                                                                                            aCATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCAtctt                                         >  2:493568/1‑150 (MQ=255)
                                                                                                              aTCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCAtcttct                                       <  1:516572/150‑1 (MQ=255)
                                                                                                               tCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCt                                                                                                <  1:1061068/92‑1 (MQ=255)
                                                                                                               tCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCt                                                                                                >  2:1061068/1‑92 (MQ=255)
                                                                                                                   aTTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTgg                                  <  1:1186814/150‑1 (MQ=255)
                                                                                                                      cGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTg                               >  2:1170390/1‑150 (MQ=255)
                                                                                                                       gggATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGATATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGc                              >  2:856192/1‑150 (MQ=255)
                                                                                                                        ggATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCTATACGCTTGCGCATCGGGCGCTATCCTGGTTTTTTCACTTCAACCACATTCAGGCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGGTGcc                             <  2:1262503/150‑1 (MQ=255)
                                                                                                                          aTTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAg                           >  2:1060985/1‑150 (MQ=255)
                                                                                                                          aTTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAg                           <  2:1356027/150‑1 (MQ=255)
                                                                                                                           ttGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGTCGCTATCCTGGTTATTTCACTTCAACCAGATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGc                          <  1:468490/150‑1 (MQ=255)
                                                                                                                           aggCCGGATAAGACGCGGCCAGCGTCGCATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTc                                                                                       <  1:652649/87‑1 (MQ=21)
                                                                                                                            tGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTATCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCa                         <  2:745488/150‑1 (MQ=255)
                                                                                                                             gCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCAc                        >  1:715687/1‑150 (MQ=255)
                                                                                                                                ggataagACGCGGCCAGCGTCGCATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGAtttt                                               <  1:86003/118‑1 (MQ=38)
                                                                                                                                ggataagACGCGGCCAGCGTCGCATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGAtttt                                               >  2:86003/7‑124 (MQ=38)
                                                                                                                                ggATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGcc                     >  1:1260311/1‑150 (MQ=255)
                                                                                                                                     cGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTg                <  2:226067/150‑1 (MQ=255)
                                                                                                                                     cGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTg                >  2:240251/1‑150 (MQ=255)
                                                                                                                                     cGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTg                <  1:1051117/150‑1 (MQ=255)
                                                                                                                                       aCGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTa              >  1:732230/1‑150 (MQ=255)
                                                                                                                                          cTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGAGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTg           >  2:1092862/1‑150 (MQ=255)
                                                                                                                                            ggCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTggg         <  2:1387889/150‑1 (MQ=255)
                                                                                                                                            ggCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTggg         <  1:540440/150‑1 (MQ=255)
                                                                                                                                             gCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGa        >  1:1284968/1‑150 (MQ=255)
                                                                                                                                             gCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGa        >  1:701495/1‑150 (MQ=255)
                                                                                                                                               cgcgTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGAtc      <  2:509682/150‑1 (MQ=255)
                                                                                                                                                   tCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCtt  <  1:608983/150‑1 (MQ=255)
                                                                                                                                                   tCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCAACTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCtt  >  2:930830/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TTGCTCATAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTT  >  NC_000913/1410808‑1411104

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