Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,660,611 G→T 12.6% Y307* (TAC→TAA) ‡ iscS ← cysteine desulfurase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,660,6110GT12.6% 331.1 / 29.6 111Y307* (TAC→TAA) ‡iscScysteine desulfurase
Reads supporting (aligned to +/- strand):  ref base G (45/52);  new base T (3/11);  total (48/63)
Fisher's exact test for biased strand distribution p-value = 9.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.18e-01

Only 100 of 113 total aligned reads displayed.

AAAGAGAAACGGATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGC  >  NC_000913/2660462‑2660757
                                                                                                                                                     |                                                                                                                                                  
aaaGAGAAACGGATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCAGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACg                                                                                                                                                    <  1:1459493/150‑1 (MQ=255)
     gaAACGGATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGtt                                                                                                                                               >  1:1589330/1‑150 (MQ=255)
         cGGATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGTGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAg                                                                                                                                           <  1:724923/150‑1 (MQ=255)
          ggATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGATACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAgg                                                                                                                                          <  1:175459/150‑2 (MQ=255)
          ggATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAgg                                                                                                                                          <  1:618591/150‑2 (MQ=255)
          ggATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAgg                                                                                                                                          <  2:1522618/150‑2 (MQ=255)
          ggATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGc                                                                                                                                          >  1:27916/1‑150 (MQ=255)
          ggATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGc                                                                                                                                          >  2:907599/1‑150 (MQ=255)
           gATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAggt                                                                                                                                         <  1:8243/150‑3 (MQ=255)
           gATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCt                                                                                                                                         >  1:871330/1‑150 (MQ=255)
            aTAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAGGTg                                                                                                                                        <  1:1286586/150‑1 (MQ=255)
            aTAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAGGTg                                                                                                                                        <  1:841733/150‑1 (MQ=255)
            aTAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTg                                                                                                                                        <  2:592690/150‑1 (MQ=255)
             tAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAGGTGa                                                                                                                                       <  2:1610556/150‑1 (MQ=255)
              agagCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAGGTGag                                                                                                                                      <  2:857494/150‑2 (MQ=255)
              agagCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGAc                                                                                                                                      <  2:1559105/150‑1 (MQ=255)
                 gCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGACTTGCTGACGTACAGGCGGAACCTGAAGAAACTCCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGtt                                                                                                                                   >  1:931747/1‑150 (MQ=255)
                  cTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCAGCTTCAACTACGTTGAAGGTGagtcgc                                                                                                                                  <  2:251258/150‑6 (MQ=255)
                  cTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTg                                                                                                                                  >  1:336932/1‑150 (MQ=255)
                  cTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTg                                                                                                                                  <  2:905463/150‑1 (MQ=255)
                   tATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTga                                                                                                                                 <  2:27916/150‑1 (MQ=255)
                   aaTGTCCCAGCTCGTCGTTCACCCCCAGCGCGCGCAGCACGTAGGCCGGCTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGGCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTga                                                                                                                                 <  2:940567/149‑1 (MQ=255)
                    aTGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTgag                                                                                                                                >  1:1387873/1‑150 (MQ=255)
                      gtgCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGaa                                                                                                                              >  2:1487947/1‑150 (MQ=255)
                       tgCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGa                                                                                                                                       <  1:1049902/140‑1 (MQ=255)
                       tgCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGa                                                                                                                                       >  2:1049902/1‑140 (MQ=255)
                       tgCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAAt                                                                                                                             >  1:830983/1‑150 (MQ=255)
                       tgCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAAt                                                                                                                             >  1:886536/1‑150 (MQ=255)
                         ccAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGt                                                                                                                           >  2:111368/1‑150 (MQ=255)
                          cAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGtt                                                                                                                          >  2:1596036/1‑150 (MQ=255)
                          cAGCACGTCGTTCACCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTCACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGtt                                                                                                                          >  2:890835/1‑150 (MQ=255)
                           aGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTc                                                                                                                         >  2:596794/1‑150 (MQ=255)
                              tcgtcgTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGgc                                                                                                                      <  1:111368/150‑1 (MQ=255)
                              tcgtcgTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGgc                                                                                                                      >  2:907871/1‑150 (MQ=255)
                                   gtTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCAcc                                                                                                                 <  1:1428218/150‑1 (MQ=255)
                                        gCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCgg                                                                                                                       >  2:684253/1‑139 (MQ=255)
                                        gCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCgg                                                                                                                       <  1:684253/139‑1 (MQ=255)
                                        gCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGtt                                                                                                            <  2:1449988/150‑1 (MQ=255)
                                           ccAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCa                                                                                                         <  2:1201566/150‑1 (MQ=255)
                                            cAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAg                                                                                                        <  1:252872/150‑1 (MQ=255)
                                             aGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACATACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAgg                                                                                                       <  1:720201/150‑1 (MQ=255)
                                               cgcgcgcAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTc                                                                                                     >  2:906184/1‑150 (MQ=255)
                                                     cagcaCGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGtt                                                                                               >  2:1275643/1‑150 (MQ=255)
                                                      agcaCGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTc                                                                                              <  2:527207/150‑1 (MQ=255)
                                                        caCGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAg                                                                                            <  1:1196906/150‑1 (MQ=255)
                                                         aCGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAgg                                                                                           <  1:1464150/150‑1 (MQ=255)
                                                          cGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGt                                                                                          <  1:1218867/150‑1 (MQ=255)
                                                           gTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTa                                                                                         >  1:937863/1‑150 (MQ=255)
                                                             aGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTaaa                                                                                       <  1:1416596/150‑1 (MQ=255)
                                                             aGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTaaa                                                                                       <  1:1411027/150‑1 (MQ=255)
                                                              ggACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTGGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAAc                                                                                      <  1:1488704/150‑1 (MQ=255)
                                                               gACGGTTCGAGGCTTGCTGACGTACAGGCGGAACGTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAAct                                                                                     <  1:800259/150‑1 (MQ=255)
                                                               gACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAAct                                                                                     <  1:1472996/150‑1 (MQ=255)
                                                                  ggTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAActtct                                                                                  <  1:596794/150‑1 (MQ=255)
                                                                    ttCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACttcttc                                                                                >  2:35814/1‑150 (MQ=255)
                                                                    ttCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACttcttc                                                                                <  1:1596036/150‑1 (MQ=255)
                                                                     tCGAGGCTTGCTGACGTACATGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCg                                                                               <  1:1619362/150‑1 (MQ=255)
                                                                      cGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGGGAACCGTGTTCCAGTTCACCGTTCAGGTAAACTTCTTCGa                                                                              <  2:336932/150‑1 (MQ=255)
                                                                      cGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGa                                                                              <  1:683516/150‑1 (MQ=255)
                                                                        aGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGata                                                                            <  2:830983/150‑1 (MQ=255)
                                                                         ggCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGatat                                                                           >  1:343748/1‑150 (MQ=255)
                                                                           cTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCt                                                                         >  1:1061092/1‑150 (MQ=255)
                                                                                   cGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGcc                                                                 >  1:188893/1‑150 (MQ=255)
                                                                                        aGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTcc                                                            <  1:572848/150‑1 (MQ=255)
                                                                                         ggCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCca                                                           >  2:1294485/1‑150 (MQ=255)
                                                                                           cGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCcaca                                                         <  1:405069/150‑1 (MQ=255)
                                                                                                ccTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACgg                                                    >  2:856952/1‑150 (MQ=255)
                                                                                                  tGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCttt                                                                       <  1:4580/129‑1 (MQ=255)
                                                                                                  tGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCttt                                                                       >  2:4580/1‑129 (MQ=255)
                                                                                                    aagaaACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTAc                                                <  2:188893/150‑1 (MQ=255)
                                                                                                     agaaACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACg                                               >  2:292277/1‑150 (MQ=255)
                                                                                                     agaaACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACg                                               >  2:1179612/1‑150 (MQ=255)
                                                                                                       aaaCTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCa                                             <  1:1294485/150‑1 (MQ=255)
                                                                                                       aaaCTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCa                                             <  1:906184/150‑1 (MQ=255)
                                                                                                        aaCTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAg                                            <  2:369856/150‑1 (MQ=255)
                                                                                                         aCTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAgg                                           <  1:33046/150‑1 (MQ=255)
                                                                                                          cTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGc                                          >  2:593341/1‑150 (MQ=255)
                                                                                                           tGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGcc                                         >  1:510242/1‑150 (MQ=255)
                                                                                                               cggTCTTTCAGCGCCATGATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGGCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCgcgc                                     <  2:475224/149‑1 (MQ=255)
                                                                                                                        aGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTcc                            >  2:33123/1‑150 (MQ=255)
                                                                                                                         gcgcCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCa                           >  2:1555659/1‑150 (MQ=255)
                                                                                                                          cgcCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGGTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCAt                          <  1:907871/150‑1 (MQ=255)
                                                                                                                           gcCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCAtc                         >  2:30593/1‑150 (MQ=255)
                                                                                                                            ccATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCAtct                        <  2:1387873/150‑1 (MQ=255)
                                                                                                                            ccATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCAtct                        <  1:43886/150‑1 (MQ=255)
                                                                                                                            ccATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCAtct                        <  1:37836/150‑1 (MQ=255)
                                                                                                                                 aTCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTc                   <  2:614997/150‑1 (MQ=255)
                                                                                                                                    aGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGcc                <  2:1453653/150‑1 (MQ=255)
                                                                                                                                    aGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGcc                <  2:1201017/150‑1 (MQ=255)
                                                                                                                                      cGACTCACCTTCAACGTAGTTGAAGCTGCCGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCaa              >  1:1614712/1‑149 (MQ=255)
                                                                                                                                       acgtCAGCTTCAACTACGTTGAAGGTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCAtc             >  1:172013/4‑150 (MQ=255)
                                                                                                                                        cgtCAGCTTCAACTACGTTGAAGGTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCAtct            >  1:543458/3‑150 (MQ=255)
                                                                                                                                         cTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCAtctc           >  2:1522916/1‑150 (MQ=255)
                                                                                                                                           cACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCtctt         >  2:1554893/1‑150 (MQ=255)
                                                                                                                                            agcTTCAACTACGTTGAAGGTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCtcttc        >  1:140594/3‑150 (MQ=255)
                                                                                                                                             ccTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCtcttct       <  1:593341/150‑1 (MQ=255)
                                                                                                                                              cTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTcttctt      >  1:260911/1‑150 (MQ=255)
                                                                                                                                                 cAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtg   <  2:937863/150‑1 (MQ=255)
                                                                                                                                                  aaCGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtgc  >  1:1278285/1‑150 (MQ=255)
                                                                                                                                                  aaCGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtgc  <  2:440837/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AAAGAGAAACGGATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGC  >  NC_000913/2660462‑2660757

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: ‑