Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 497,855 C→G 6.3% intergenic (+36/‑200) adk → / → hemH adenylate kinase/ferrochelatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913497,8550CG 6.3% 436.7 / 12.8 135intergenic (+36/‑200)adk/hemHadenylate kinase/ferrochelatase
Reads supporting (aligned to +/- strand):  ref base C (62/64);  new base G (6/3);  total (68/67)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.36e-01

Only 100 of 136 total aligned reads displayed.

TGATCGGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTT  >  NC_000913/497707‑498003
                                                                                                                                                    |                                                                                                                                                    
 gATCGGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCgg                                                                                                                                                    >  2:591370/1‑150 (MQ=255)
  aTCGGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCggg                                                                                                                                                   <  2:1693675/150‑1 (MQ=255)
   tCGGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCgggg                                                                                                                                                  >  2:1224067/1‑150 (MQ=255)
   tCGGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCgggg                                                                                                                                                  <  2:1593787/150‑1 (MQ=255)
   tCGGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCgggg                                                                                                                                                  <  2:354173/150‑1 (MQ=255)
    cGGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGga                                                                                                                                                 >  2:1563349/1‑150 (MQ=255)
    cGGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGga                                                                                                                                                 <  2:711074/150‑1 (MQ=255)
       ctactactaCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGagag                                                                                                                                              <  2:1046053/150‑1 (MQ=255)
       ctactaCTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGagag                                                                                                                                              <  1:811913/150‑1 (MQ=255)
        tactacTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAgg                                                                                                                                             <  1:1654583/150‑1 (MQ=255)
         actactCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAggg                                                                                                                                            <  1:834173/150‑1 (MQ=255)
          ctactCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGt                                                                                                                                           >  1:1271326/1‑150 (MQ=255)
            actCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGCTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTa                                                                                                                                         <  1:1147097/150‑1 (MQ=255)
               ccAAAGTAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAggg                                                                                                                                      >  1:972459/1‑150 (MQ=255)
               ccAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAggg                                                                                                                                      <  1:1068444/150‑1 (MQ=255)
               ccAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAggg                                                                                                                                      >  1:1463951/1‑150 (MQ=255)
                cAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGt                                                                                                                                     >  2:377970/1‑150 (MQ=255)
                     aagcagcagcGGGTAATACCAAATAGGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTTATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGTGAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAggg                                                                                                                                <  2:254078/150‑1 (MQ=39)
                      agcagaagcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCTAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAgggg                                                                                                                               <  2:287413/150‑1 (MQ=255)
                       gcagaagcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGa                                                                                                                              >  1:695845/1‑150 (MQ=255)
                         agaagcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAAc                                                                                                                            >  1:1038845/1‑150 (MQ=255)
                         agaagcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAAc                                                                                                                            >  1:670313/1‑150 (MQ=255)
                          gaagcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACa                                                                                                                           <  1:721437/150‑1 (MQ=255)
                            agcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAgg                                                                                                                         <  2:218732/150‑1 (MQ=255)
                            agcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAgg                                                                                                                         >  2:1627689/1‑150 (MQ=255)
                            agcGGGTAATACCAAATACGAGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAgg                                                                                                                         <  2:855637/150‑1 (MQ=255)
                                 gTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGATAAAATCATCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGCTTAGGGTGAGGGGATCAGGCCCGc                                                                                                                    >  2:92530/1‑150 (MQ=255)
                                   aaTACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTTAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGcac                                                                                                                  <  1:1050691/150‑1 (MQ=255)
                                    aTACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGcaca                                                                                                                 <  2:1463951/150‑1 (MQ=255)
                                    aTACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGcaca                                                                                                                 <  2:1691599/150‑1 (MQ=255)
                                    aTACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGcaca                                                                                                                 <  2:972459/150‑1 (MQ=255)
                                    aTACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGcaca                                                                                                                 >  1:1603324/1‑150 (MQ=255)
                                      aCCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGCCAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAg                                                                                                               >  1:748249/1‑150 (MQ=255)
                                       ccAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGc                                                                                                              >  1:311275/1‑150 (MQ=255)
                                          aaTACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGGCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCaaa                                                                                                           <  2:756979/150‑1 (MQ=255)
                                                cgAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATc                                                                                                     <  1:591370/150‑1 (MQ=255)
                                                  aaaGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAAATCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAg                                                                                                   <  1:543501/150‑1 (MQ=255)
                                                   aaGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCATTCGGCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGc                                                                                                  <  2:375642/150‑1 (MQ=255)
                                                    aGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGAGCTGATCTGGAAAAAATCCTCGGCTAATTCGAACGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCCCAAAAGCAAACtttgccgaa                                                                                                 >  1:300906/1‑143 (MQ=38)
                                                        gACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAAtct                                                                                             >  1:159418/1‑150 (MQ=255)
                                                        gACGGCACCAAGCCGGTTGCTGAAGGTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTGGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAAtct                                                                                             <  2:51837/150‑1 (MQ=255)
                                                                 aaGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCCCGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGa                                                                                    <  2:602062/150‑1 (MQ=255)
                                                                     cGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATatt                                                                                <  1:1224067/150‑1 (MQ=255)
                                                                          gCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCa                                                                                           >  2:953948/1‑134 (MQ=255)
                                                                          gCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCa                                                                                           <  1:947850/134‑1 (MQ=255)
                                                                          gCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATattcattc                                                                           >  2:1684794/1‑150 (MQ=255)
                                                                             gAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACTAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGc                                                                        >  1:144650/1‑150 (MQ=255)
                                                                             gAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGc                                                                        <  1:1349969/150‑1 (MQ=255)
                                                                             gAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGCGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTCTCAGCAATCTCAGGCCGGATATTCATTCGGc                                                                        >  1:251880/1‑150 (MQ=255)
                                                                                   cgcgCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTa                                                                  >  2:1514245/1‑150 (MQ=255)
                                                                                   cgcgCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTa                                                                  <  2:1387175/150‑1 (MQ=255)
                                                                                      gcTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACaa                                                               <  2:901384/150‑1 (MQ=255)
                                                                                         gATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACaaaaa                                                            >  1:99587/1‑150 (MQ=255)
                                                                                         gATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACaaaaa                                                            >  2:17902/1‑150 (MQ=255)
                                                                                         gATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTCAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACaaaaa                                                            >  1:929438/1‑150 (MQ=255)
                                                                                           tCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAAcc                                                          >  2:596707/1‑150 (MQ=255)
                                                                                              ggAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCa                                                       <  2:1178209/150‑1 (MQ=255)
                                                                                               gAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCaa                                                      <  1:1627689/150‑1 (MQ=255)
                                                                                               gAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCaa                                                      >  2:1493087/1‑150 (MQ=255)
                                                                                                aaaaaaTCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAAt                                                     <  2:695845/150‑1 (MQ=255)
                                                                                                 aaaaaTCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATc                                                    <  2:1126976/150‑1 (MQ=255)
                                                                                                  aaaaTCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCt                                                   <  1:1563349/150‑1 (MQ=255)
                                                                                                    aaTCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGc                                                                                                  <  1:933417/101‑1 (MQ=255)
                                                                                                    aaTCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGc                                                                                                  >  2:933417/1‑101 (MQ=255)
                                                                                                     aTCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTgcg                                                <  2:1614384/150‑1 (MQ=255)
                                                                                                      tCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTgcgc                                               <  1:798815/150‑1 (MQ=255)
                                                                                                      tCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTgcgc                                               <  2:311275/150‑1 (MQ=255)
                                                                                                        cTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTa                                             >  2:666970/1‑150 (MQ=255)
                                                                                                          cGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTAtt                                           <  1:762623/150‑1 (MQ=255)
                                                                                                           ggCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTc                                          <  2:35395/150‑1 (MQ=255)
                                                                                                             cTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCt                                        <  2:1484949/150‑1 (MQ=255)
                                                                                                             cTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCt                                        >  1:385747/1‑150 (MQ=255)
                                                                                                             cTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACATCATTCTGCGAGATTcat                                        >  1:235606/1‑148 (MQ=255)
                                                                                                               aaTTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCttt                                      >  1:1694887/1‑150 (MQ=255)
                                                                                                                aTTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCtttt                                     <  2:1173819/150‑1 (MQ=255)
                                                                                                                aTTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCtttt                                     <  1:3560/150‑1 (MQ=255)
                                                                                                                aTTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCtttt                                     >  1:1195793/1‑150 (MQ=255)
                                                                                                                 ttCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTAATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCttttt                                    <  1:163844/150‑1 (MQ=255)
                                                                                                                   cGAAAGCGCGCACGGACAGTCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTg                                 <  2:1082684/150‑1 (MQ=255)
                                                                                                                                 gACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCACGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATCTTCCTTCGGCCTTTTACAAAAGCCTCACTCTGCGGTATTCCGTTTTCTGCGTGGACCTCTca                    >  1:274472/1‑150 (MQ=37)
                                                                                                                                   cAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAACCCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTcaca                  >  1:1361619/1‑150 (MQ=255)
                                                                                                                                     gTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGc                >  1:978359/1‑150 (MQ=255)
                                                                                                                                      tCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGct               <  1:1665254/150‑2 (MQ=255)
                                                                                                                                       acccTCTCCCCGAGGGGGCGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCaa              >  2:1616051/2‑150 (MQ=255)
                                                                                                                                       acccTCTCCCCGAGGGGGCGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCaa              >  2:641784/2‑150 (MQ=255)
                                                                                                                                        cccTCTCCCCGAGGGGGCGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAAt             >  1:1269518/1‑150 (MQ=255)
                                                                                                                                        cccTCTCCCCGAGGGGGCGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAAt             >  2:1231467/1‑150 (MQ=255)
                                                                                                                                        cccTCTCCCCGAGGGGGCGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAAt             >  2:1225552/1‑150 (MQ=255)
                                                                                                                                        cccTCTCCCCGAGGGGGCGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGAGTTGACCTCTCACAGCAAt             >  1:1262823/1‑150 (MQ=255)
                                                                                                                                        cccTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAAt             >  1:535133/1‑150 (MQ=255)
                                                                                                                                        cccTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAAt             >  1:40257/1‑150 (MQ=255)
                                                                                                                                         ccTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAAtt            >  2:833989/1‑150 (MQ=255)
                                                                                                                                          cTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTa           >  1:364011/1‑150 (MQ=255)
                                                                                                                                            cGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGt         >  1:1406952/1‑150 (MQ=255)
                                                                                                                                              cccccTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGttc       <  2:1603324/150‑1 (MQ=255)
                                                                                                                                               ccccTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGttct      >  2:878932/1‑150 (MQ=255)
                                                                                                                                                 ccTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGttcttc    >  1:1537307/1‑150 (MQ=255)
                                                                                                                                                  cTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGttcttct   >  1:465813/1‑150 (MQ=255)
                                                                                                                                                  cTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGttcttct   >  1:472277/1‑150 (MQ=255)
                                                                                                                                                   gcGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCCGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACTGCAATTAGttcttctt  <  2:236355/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGATCGGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTT  >  NC_000913/497707‑498003

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