Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 497,854 T→A 5.8% intergenic (+35/‑201) adk → / → hemH adenylate kinase/ferrochelatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913497,8540TA 5.8% 379.6 / 8.5 135intergenic (+35/‑201)adk/hemHadenylate kinase/ferrochelatase
Reads supporting (aligned to +/- strand):  ref base T (62/65);  new base A (6/2);  total (68/67)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

Only 100 of 136 total aligned reads displayed.

GCTGATCGGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCT  >  NC_000913/497705‑498002
                                                                                                                                                     |                                                                                                                                                    
gCTGATCGGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCt                                                                                                                                                      <  2:1551371/150‑1 (MQ=255)
   gATCGGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCgg                                                                                                                                                   >  2:591370/1‑150 (MQ=255)
     tCGGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCgggg                                                                                                                                                 <  2:1593787/150‑1 (MQ=255)
     tCGGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCgggg                                                                                                                                                 <  2:354173/150‑1 (MQ=255)
     tCGGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCgggg                                                                                                                                                 >  2:1224067/1‑150 (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:1463951/1‑150 (MQ=255)
                 ccAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAggg                                                                                                                                     <  1:1068444/150‑1 (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:670313/1‑150 (MQ=255)
                           agaagcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAAc                                                                                                                           >  1:1038845/1‑150 (MQ=255)
                            gaagcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACa                                                                                                                          <  1:721437/150‑1 (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)
                                          cAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCa                                                                                                            <  2:176336/150‑1 (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                                                           >  2:17902/1‑150 (MQ=255)
                                                                                           gATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACaaaaa                                                           >  1:99587/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)
                                                                                                   aaaaaTCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATc                                                   <  2:1603953/150‑1 (MQ=255)
                                                                                                    aaaaTCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCt                                                  <  1:1563349/150‑1 (MQ=255)
                                                                                                      aaTCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGc                                                                                                 >  2:933417/1‑101 (MQ=255)
                                                                                                      aaTCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGc                                                                                                 <  1:933417/101‑1 (MQ=255)
                                                                                                       aTCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTgcg                                               <  2:1614384/150‑1 (MQ=255)
                                                                                                        tCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTgcgc                                              <  2:311275/150‑1 (MQ=255)
                                                                                                        tCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTgcgc                                              <  1:798815/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                                       >  1:385747/1‑150 (MQ=255)
                                                                                                               cTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCt                                       <  2:1484949/150‑1 (MQ=255)
                                                                                                               cTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACATCATTCTGCGAGATTcat                                       >  1:235606/1‑148 (MQ=255)
                                                                                                                 aaTTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCttt                                     >  1:1694887/1‑150 (MQ=255)
                                                                                                                  aTTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCtttt                                    >  1:1195793/1‑150 (MQ=255)
                                                                                                                  aTTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCtttt                                    <  1:3560/150‑1 (MQ=255)
                                                                                                                  aTTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCtttt                                    <  2:1173819/150‑1 (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)
                                                                                                                                       gTCCCCTCGCCCCCAGGGGGAGAGGGTTATGCTTAGGGGAACAGGCCCGCACAAGCAAAGTTATCAGCAATCTCAGTCCGGATATTCAATCTGCCTTTTACAAAAACCTCAATCTGCCCTAATCCTTTATCTGATTTTACCTATCACAGc               <  2:225966/150‑1 (MQ=255)
                                                                                                                                        tCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGct              <  1:1665254/150‑2 (MQ=255)
                                                                                                                                         acccTCTCCCCGAGGGGGCGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCaa             >  2:1616051/2‑150 (MQ=255)
                                                                                                                                         acccTCTCCCCGAGGGGGCGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCaa             >  2:641784/2‑150 (MQ=255)
                                                                                                                                          cccTCTCCCCGAGGGGGCGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAAt            >  2:1231467/1‑150 (MQ=255)
                                                                                                                                          cccTCTCCCCGAGGGGGCGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAAt            >  2:1225552/1‑150 (MQ=255)
                                                                                                                                          cccTCTCCCCGAGGGGGCGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAAt            >  1:1269518/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:1651744/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)
                                                                                                                                                 ccccTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGttct     >  2:878932/1‑150 (MQ=255)
                                                                                                                                                   ccTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGttcttc   >  1:1537307/1‑150 (MQ=255)
                                                                                                                                                    cTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGttcttct  >  1:472277/1‑150 (MQ=255)
                                                                                                                                                    cTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGttcttct  >  1:465813/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCTGATCGGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCT  >  NC_000913/497705‑498002

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: